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Preface 


This book deals with the challenge of exploiting ambient vibrational energy, which 
can be used to power small and low-power electronic devices, e.g., wireless sensor 
nodes. Generally, particularly for low voltage amplitudes, low-loss rectification is 
required to achieve high conversion efficiency. In the special case of piezoelectric 
energy harvesting, pulsed charge extraction has the potential to extract more power 
compared to a full-bridge rectifier. Therefore, a fully autonomous CMOS inte- 
grated interface circuit for piezoelectric generators which fulfills these require- 
ments is presented. This book covers three main aspects of the integrated interface 
circuit: 

First of all, the book explains in detail the different circuit blocks on transistor 
level and highlights techniques to reduce the power consumption. Hence, only a 
very small fraction of the power delivered by the generator is wasted, which is 
extremely important in order to achieve a high overall harvesting efficiency, 
especially in case the piezoelectric generator outputs little power. 

Second, the book analyzes the various loss mechanisms within the CMOS chip, 
such as conduction losses, switching losses, etc. Therefore, a mathematical method 
of approximating the conduction losses is presented, which reduces calculation 
effort and gives deep insight into the loss dependency on different parameters. A 
detailed breakdown of the actual chip losses identifies the most dominant loss 
mechanisms and gives ideas how these losses can be further reduced. 

Third, since the performance of the CMOS chip strongly depends on the used 
power source, lot of effort is spent on investigating the interaction between the 
interface circuit and the piezoelectric generator. For accurate simulations, a model 
which takes into account this electromechanical feedback is used. A CMOS chip 
has been fabricated and tested under laboratory conditions in combination with one 
custom-made and one commercially available piezoelectric generator. By com- 
paring measurement and simulation results, the used model could be verified. 

The presented CMOS chip has been shown to be fully autonomous and self- 
powered down to a piezoelectric output power in the range of 10 uW. It enables 
cold-startup and enhances the extracted power compared to the commonly known 
diode rectifiers by up to 127 %, depending on the excitation conditions. For low 
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excitations, due to the boosting effect, the chip harvests power where a diode 
rectifier would harvest nothing. The chip operates properly for piezoelectric 
voltage amplitudes in the range of 1.3—20 V and for excitation frequencies from 50 
Hz to 2 kHz. 

Due to these key properties enabling universal usage, other CMOS designers 
working in the field of energy harvesting will be encouraged to use some of the 
shown structures for their own implementations. The book highlights the design 
process from scratch to the final chip. Therefore, it gives the designer a compre- 
hensive guide of how to 


e setup an appropriate harvester model to get realistic simulation results, 

e design the integrated circuits for low power operation, 

e setup a laboratory measurement environment in order to extensively charac- 
terize the chip in combination with the real harvester, 

e and finally interpret the simulation and measurement results in order to improve 
the chip performance. 


Since the dimensions of all devices (transistors, resistors etc.) are given, readers 


and other designers can easily re-use the presented circuit concepts. 


Villingen-Schwenningen, Germany Thorsten Hehn 
Freiburg, Germany Yiannos Manoli 
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Chapter 1 
Introduction 


In the last few years, energy harvesters have decreased in size, and at the same time 
increased their output power. These properties make energy harvesting attractive for 
powering wireless sensor nodes (WSNs) which are originally battery powered, with 
the goal of prolongating their lifetime. This chapter starts with a rough review of the 
main energy harvesting mechanisms, covering the thermoelectric, radio frequency 
(RF) and vibration based principles. For each mechanism, some representative appli- 
cation examples are shown, followed by possible applications of these conversion 
mechanisms in WSNs are presented. After that, a detailed review of the state of the 
art in interface circuits for vibration based energy harvesters is given. Therefore, the 
interface circuits are separated in two categories, one of them focusing on efficient 
AC/DC conversion, and the others incorporating impedance matching methods. 

In order to classify the major achievements of the proposed interface circuit, a list 
of aspects which should be fulfilled when designing general interface circuitry for 
energy harvesters is given. Based on this list, the major achievements of this work 
are then discussed. This chapter ends up with a brief description of the organization 
of the book. 


1.1 Energy Harvesting Principles 


Energy harvesters! convert ambient energy into usable electrical energy. There 
are many possible conversion mechanisms, which are briefly presented in the 
following [49]. 


! Usually, the term energy harvester or energy scavenger is used for a microscale device converting 
small amounts of ambient energy into electrical energy. In the broader sense, such a device is of 
course a generator, and sometimes is referred to as a converter or a transducer as well. In this book, 
the terms energy harvester and generator are used as synonyms. 
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2 1 Introduction 





Fig. 1.1 Illustration of a sensor cube application integrating photovoltaic harvesters as packaging 
[34] (reproduced with kind permission of Valer Pop) 


The probably most popular ambient energy source is sunlight which is converted 
into electrical energy by means of photovoltaic cells. Outdoors, they provide a power 
density of 100 mW /cm?, which reduces to 10 — 100 pW /cm? for indoor application. 
Their efficiencies range from 5—30 96, depending on the material used. Hence, power 
densities provided by photovoltaic cells usually are in the order of 10 mW outdoors 
and 10 uW indoors. Figure 1.1 shows a sensor cube powered by photovoltaic cells, 
providing 20 uW in indoor illuminating conditions. 

Temperature gradients can also be converted into electrical energy, by exploiting 
the Seebeck effect: In case two ends of a piece of thermoelectric material are kept 
at a different temperature, i.e. a temperature gradient is present, an electrical voltage 
appears between the ends. A thermocouple made of two different materials and a 
metallic interconnect is the simplest thermoelectric generator. If several thermocou- 
ples are connected thermally in parallel and electrically in series, a thermogener- 
ator with useable output voltage and power levels can be realized. The company 
Micropelt produces small-scale thermoelectric generators with a surface area of 
14 mm^ and a thickness of 1mm, producing a matched power of 1.5mW at 
AT = 10K [29]. On one hand, these small generators allow the implementation 
of small WSNs which have a long lifetime and are reliable due to the absence of 
moving parts. But on the other hand, it is difficult to achieve a significant temperature 
gradient over the small thickness, requiring large heat sinks, increasing total system 
size. A photograph of Micropelt generator chips is depicted in Fig. 1.2. 

Another possible source for energy harvesting are radio frequency (RF) waves 
available everywhere due to public telecommunication services like Global System 
for Mobile Communications (GSM) and Wireless Local Area Network (WLAN). 
This conversion mechanism exhibits relatively low power densities. For distances of 
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Fig. 1.2 Micropelt thermoelectric generator chip MPG-D751/MPG-D651 [29] (reproduced with 
kind permission of Micropelt) 





Fig. 1.3 RF harvesting system with a small loop antenna (reproduced with kind permission of 
Tolgay Ungan) 


25 — 100 m from a GSM base station, a power of 0.1 — 1 mW /m? may be achieved 
from single frequencies. WLAN exhibits one magnitude lower power densities [48]. 
Alternatively, a dedicated RF source can be positioned in front of a WSN in order 
to directly power this device at a small distance. Using this method, a transmission 
power of 200 pW at a distance of 2 m has been reported [46]. Figure 1.3 shows the 
corresponding RF system. 

As a last conversion mechanism, harvesting energy from vibrations or motions is 
discussed in the following. For that purpose, mainly three conversion mechanisms 
exist: The capacitive [18, 52], the inductive [5, 41] and the piezoelectric [15, 20, 38] 
principle. The capacitive principle exploits the relative movement of two capaci- 
tor plates with respect to each other. This movement causes charges stored in the 
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Fig. 1.4 Capacitive microgenerator chips, partly packaged in CLCC packages (reproduced with 
kind permission of Daniel Hoffmann) 





Fig. 1.5 SEM picture of a single triangular-curved piezoelectric cantilever which is an optimum 
for a uniform load (reproduced with kind permission of Ingo Kühne) 


capacitor plates to be moved through an external load, producing electrical power. 
Figure 1.4 shows an example of a capacitive microgenerator packaged in a ceramic 
lead chip carrier (CLCC) package. For the inductive conversion principle, a perma- 
nent magnet moves relatively to a coil, inducing a voltage across the coil. Lastly, 
mechanical strain in a layer made of piezoelectric material produces a charge sepa- 
ration in the piezoelectric material which turns into a voltage between both surfaces 
of the layer. Figure 1.5 depicts a SEM image of a triangular-curved piezoelectric 
cantilever. 
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Each of the three main kinetic conversion principles (capacitive, inductive, piezo- 
electric) has its own advantages and disadvantages, which are summarized in the 
following [4, 40]: 


e Compared to the other two principles, the capacitive conversion principle can be 
realized most easily in micro-electro-mechanical systems (MEMS) technology, 
and there exists a high level of corresponding process know-how. But unfortu- 
nately, capactive harvesters require an initial polarizing voltage or charge. This 
can be achieved by special materials called electrets which can provide the 
initial polarization and these can maintain their charge over several years [3]. 
Electrostatic harvesters suffer from a high output impedance, limiting the achiev- 
able output current. In contrast, the output voltage is usually very high (— 100 V), 
posing a challenge in the interface circuit design. Parasitic capacitances within the 
structure can sometimes lead to reduced generator efficiency and cause capacitor 
electrodes shorting or sticking. Electrostatic generators exhibit by far the lowest 
energy storage density [40]. 

e The inductive conversion principle offers a well-established technique, with a 
variety of spring/mass configurations that can be used with various types of mate- 
rial. High-performance bulk magnets and multi-turn, macro-scale coils are readily 
available. A relatively high output current comes at the expense of a low output 
voltage («1 V), requiring highly efficient rectifiers and boost converters to provide 
a voltage level suitable to supply common electronics. Problems exist in designing 
MEMS devices, due to the poor properties of planar magnets, the limited number 
of turns of planar coils, and the restricted vibration amplitude. The practical maxi- 
mum energy density of inductive generators is much higher compared to capacitive 
generators, but lower than what is achieved by piezoelectric generators [40]. 

e The piezoelectric conversion principle offers the simplest approach, because there 
is no need for having complex geometries and numerous extra components. Vibra- 
tions are directly converted into electricity through the electroded piezoelectric 
material. Piezoelectric materials can be simply deposited using thin- and thick- 
film, hence it is well suited for MEMS processes. Piezoelectric harvesters are 
capable of producing relatively high output voltages but only at low electrical cur- 
rents. Due to the fact that piezoelectric materials are strained directly, the piezo- 
electric properties limit overall performance and lifetime. The commonly used 
material lead circonate titanate (PZT) is very brittle and hence prone to crack if it 
is over-stressed. The piezoelectric conversion principle is discussed more detailed 
in Chap. 2. According to Roundy et al. [40], of the three kinetic conversion mecha- 
nisms, piezoelectric generators offer the highest practical maximum energy storage 
density. 


According to [4, 39, 40], vibration-based energy harvesting is a viable means of 
obtaining the small quantities of energy necessary to power WSNs. The three main 
techniques of harvesting energy from ambient vibrations have been shown to be 
capable of generating output power levels in the range of a few microwatts to several 
hundred microwatts. 


6 ] Introduction 


This book focuses on an integrated interface circuit for piezoelectric harvesters. 
Hence, only the piezoelectric conversion principle is considered from Chap. 2 on. 


1.2 Examples of Wireless Sensor Nodes 


Over the last few years, wearable WSNs have become increasingly attractive in 
a wide range of applications [49, 50]. A general WSN is composed of an energy 
management which manages the converted energy from the ambient, the sensors, a 
signal processing unit, and an RF radio. Figure 1.6 shows the corresponding block 
diagramm of a WSN. The RF radio transmission usually consumes much more 
power than the measurement system. In most cases, the duty cycle of the WSN is 
low, meaning that there are large time intervals between the measurements, in order 
to save power. 

Body area networks (BANs) provide users with information about biological or 
physiological signals for medical, lifestyle, assisted living, sports, or entertainment 
purposes. An example device that measures the electromyogram (EMG) on human 
body is shown in Fig. 1.7. Powered by photovoltaic cells, the EMG device measures 
the “muscle activity, fatigue or power, with possible applications in sports, lifestyle, 
healthcare and rehabilitation" [34]. The EMG signals are sampled at 200 Hz, and the 
gathered data is transmitted wirelessly to a base station. 
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Fig. 1.6 Block diagram of a WSN 
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Fig. 1.7 Wireless EMG BAN application [34] (reproduced with kind permission of Valer Pop) 


In the automotive sector, WSNs are used in safety and control applications, e.g. 
in tire pressure monitoring systems (TPMSs). Generally, the TPMS can be mounted 
onto the tire as an attachable patch, or it could be integrated into the rubber itself 
as a part of the vulcanizing process. In the United States, the usage of TPMSs is 
regulated by the Department of Transportation (DOT). In the Transportation Recall 
Enhancement, Accountability, and Documentation (TREAD) Act of 2000 [47], the 
DOT has established "[. . .] anew Federal Motor Vehicle Safety Standard that requires 
the installation of tire pressure monitoring systems (TPMSs) that warn the driver 
when a tire is significantly under-inflated. The standard applies to passenger cars, 
trucks, multipurpose passenger vehicles, and buses with a gross vehicle weight rating 
of 10,000 pounds or less, except those vehicles with dual wheels on an axle. [. . .]" 
As of November 1, 2006, “[...] a vehicle's TPMS must warn the driver when the 
pressure in any single tire or in each tire in any combination of tires, up to a total of 
four tires, has fallen to 25 9o or more below the vehicle manufacturer's recommended 
cold inflation pressure for the tires, or a minimum level of pressure specified in the 
standard, whichever pressure is higher. [. ..]" 

Smart buildings can help reducing overall power consumption by detecting the 
presence or absence of human beings and turn on or off lighting, heating, ventilation, 
or air-conditioning (HVAC) accordingly. An example for making buildings smart is 
the wireless switch produced by EnOcean, as shown in Fig. 1.8. Inductive energy 
conversion is achieved by pushing a magnet through a coil. The power generated 
is more than 100 J per push, which is sufficient to send an RF signal to the base 
station which turns the light on or off. 

Predictive maintenance can increase the safety of infrastructures like bridges, 
where strain can be monitored by wireless sensors in order to detect over-stress 
[14]. This offers the possibility to react before a bridge collapses. In order to extend 
lifetime, low-power radios have to be used. Predictive maintenance is also an issue 
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Fig. 1.8 EnOcean's self-powered wireless light switch PTM200 [12] (reproduced with kind per- 
mission of EnOcean) 


in public transportation vehicles like airplanes or trains. For example, aircrafts are 
subjected to various stresses during their flights, especially mechanical stresses dur- 
ing take-offs and landings. Health and usage monitoring systems developed e.g. by 
EADS France Innovation Works improve maintenance and can consequently reduce 
costs [37]. 

Another important application is the identification or tracking functionality 
enabled by radio frequency identification (RFID) tags. RFID tags may be active 
or passive. Passive RFID tags are powered by RF energy supplied by the base station 
or reader, while active RFID tags are partially or completely powered by a battery. 
In addition to identification or tracking, RFID tags can also be used for sensing. For 
example, an active RFID tag which is able to sense and log temperature and humid- 
ity in a intercontinental fresh fish logistic chain is reported in [1]. Using a 25 mAh 
lithium battery, a lifetime of 29 days is achieved, taking a measurement every 6s. 

Advances in low-power electronics have reduced the power consumption of signal 
processing units in WSNs to below 100 uW [27, 35]. In that regime, the lifetime 
of batteries is limited to a few months [13], requiring frequent battery replacement 
which is unpractical or sometimes even impossible. It is not expected that a sig- 
nificant improvement of the power density of batteries can be achieved in the near 
future [13, 32, 49]. Since energy harvesters in practical usable forms are capable of 
providing an output power in the range 10 — 100 uW [6], it seems viable that energy 
harvesting as described in Sect. 1.1 can significantly extend the lifetime of WSNs, 
in theory towards unlimited lifetime. 
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In the following, it is motivated why interface circuits are essential for energy har- 
vesters in general. After that, a review of the research conducted in the past years 
in this field is given. Therefore, two groups of interface circuits are investigated: 
The interface circuits focusing on efficient AC/DC conversion on one hand, and the 
interface circuits based on impedance matching on the other hand. 


1.3.1 About the Importance of Interface Circuits 


Generally, any electronic circuit requires a direct current (DC) to operate properly. 
There are approaches to make signal processing circuits directly operable by an alter- 
nating current (AC) [2], but these belong to absolute fundamental research and have 
no practical relevance up to now. Hence, for the case of vibration based harvesters 
providing an AC voltage, it is directly evident that a rectifier is required to convert 
AC to DC. 

The second target of interface circuits is to optimize the output power which can 
be extracted from the harvester, 1.e. to provide some kind of impedance matching. 
For harvesters with DC output characteristics, e.g. photovoltaic and thermoelectric 
harvesters, there exists an optimal output voltage at which the power transferred to the 
load is optimal. This optimal voltage is called the maximum power point (MPP). The 
MPP depends on the open circuit voltage of the harvester, which on its part depends 
on the energy source, i.e. the sunlight power in case of the photovoltaic cell or the 
temperature gradient in case of the thermoelectric generator. Hence, the interface 
circuit has to perform a maximum power point tracking (MPPT) [11, 19] in order 
to maintain the MPP under any circumstances. For vibration based harvesters, there 
also exists a MPP in case a rectifier first converts the AC voltage into DC [26, 31]. 
Additionally, the electrical damping presented by the load circuit greatly affects the 
output power. Thus, the interface circuit should set the optimal damping in order to 
extract most of the available power out of the harvester. 

As a third goal, it is desirable that the interface circuit decouples the load from the 
harvester. Often, the load electronics (e.g. WSNs) need to be powered periodically 
only for a short time interval, so the energy has to be stored in an intermediate buffer, 
e.g. a (super-) capacitor or a rechargeable battery. This enables the load circuit to 
take small amounts of energy out of the buffer whenever it is needed, independently 
on the current output power produced by the harvester. During time intervals with 
no load activity, the harvester can recharge the buffer. 
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1.3.2 Literature Overview 


Some years ago, the research on energy harvesting focused on the evaluation and 
the development of generator structures itself, improving output power, but often 
neglecting the fact that almost no application can be powered directly by the harvester 
output. Recently, the design of efficient interface circuits for any type of energy 
harvester has attracted increased attention, and a vast amount of publications have 
emerged. Since this work focuses on the design of an interface circuit for piezoelectric 
generators which are AC type, only interface circuits made for AC type vibration 
based harvesters are considered in the following. 

In general, the interface circuits reported in the last 10 years can be sepa- 
rated in two groups. The first group focuses on the efficient AC/DC conversion of 
vibration-based harvesters, by reducing losses in the rectifier and designing 
control circuitry consuming little power. Hence, a majority of the published work 
is based on integrated complementary metal-oxide semiconductor (CMOS) design, 
which is very well suited for meeting these requirements. Many authors focus on 
overcoming the voltage drop of diodes originally used for very basic rectifiers. The 
interface circuits being part of the first group can often be used for any vibration 
based harvester, i.e. no matter if capacitive, inductive, or piezoelectric, because load 
impedance matching issues aiming at a special type of harvester is rarely an issue. 

The second group concentrates on extracting the optimal power out of the 
harvester. Whereas some interface circuits provide impedance matching using a 
MPPT scheme, others exploit the capacitive output characteristics of piezoelectric 
generators by means of a nonlinear power processing technique. Due to the higher 
complexity compared to the first group, some early publications of the second group 
introducing the basic operation principle were implemented discretely and hence 
required an external power supply. However, in the meanwhile, custom integrated 
CMOS design allows self-powered standalone operation. 


1.3.2.1 Interface Circuits with the Focus on Efficient AC/DC Conversion 


All of the publications focusing on efficient AC/DC conversion summarized in 
Table 1.1 try to overcome the voltage drop over the diodes usually being part of 
rectifiers by replacing them by metal oxide-semiconductor field effect transistor 
(MOSFET) switches. The challenge is to make the involved control circuitry ultra 
low power, in order to transfer as much power as possible to the output. 

Active full-wave rectifiers are full-wave recitifiers where the passive diodes are 
substituted by actively driven MOSFETs. In Le et al. [23], compare several half-wave 
and full-wave passive and active rectifiers in terms of efficiency and output power. 
It turns out that despite the fact that the active rectifier requires quiescent current for 
operation, it has a higher peak efficiency than the passive rectifier (0.56 compared 
to 0.66). The power management system presented in [8] by Colomer-Farrarons et 
al. combines another active full-wave rectifier with a power conditioning circuit. 
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By keeping the voltage of the buffer capacitor within a fixed voltage window 1.2 — 
1.4 V, a very low power WSN can be powered. In Peters et al. [33], published another 
active rectifier with special focus on low voltage operation. By means of bulk-input 
comparators, input voltages down to 0.38 V are accepted, verified by an inductive 
microgenerator. In Maurath et al. [28], use a common-gate stage in order to achieve 
a wide range of tolerated input voltages 0.48 — 3.3 V. In combination with a MPPT 
system, energy harvesting with inductive transducers is possible at buffer voltage 
levels above actual transducer output voltages. Due to high losses in the integrated 
capacitances, the MPP tracker suffers from comparably low efficiency of below 0.48. 

The active rectification principle can also be applied to voltage doublers. In 
Dallago et al. [9], present an active voltage doubler with a peak efficiency of 0.92 by 
keeping the current consumption of the involved active elements to 500 nA. Another 
voltage doubler topology providing a positive and a negative output voltage has 
been investigated by Cheng et al. in [7]. Although this active voltage doubler has 
been implemented as discrete-type electronics, it shows relatively high efficiency of 
0.85 and tolerates very low input voltage amplitudes down to 5 mV. Furthermore, 
the operation principle has been tested with irregular input voltages. However, an 
external voltage supply is needed to power the involved comparators, limiting the 
practicability in standalone energy harvesting applications. 


1.3.2.2 Interface Circuits with the Focus on Impedance Matching 


The interface circuits listed in Table 1.2 can be roughly categorized in two groups. 
The first group consists of the circuits performing a MPP tracking. In most cases, with 
varying harvester open circuit voltage, the duty cycle of DC/DC step-down convert- 
ers is adjusted such that the harvested power remains maximum. Although the given 
circuits have been designed for piezoelectric harvesters, the corresponding method 
can also be applied to other vibration based or DC-type harvesters. The second group 
performs nonlinear extraction techniques, actively actuating the piezoelectric device 
at dedicated instants by means of one or more switches. This can be done either 
with or without an additional inductor acting as a temporary energy storage. The 
techniques making use of an inductor are called Synchronized Switch Harvesting on 
Inductors (SSHI) and Synchronous Electric Charge Extraction (SECE), as theoreti- 
cally discussed in [24]. In case a simple switch is placed in parallel to the piezoelectric 
generator, a Bias-flip Rectifier [36] is established. Lastly, a quite new concept which 
supersedes a rectifier is presented by Kwon and Rincón-Mora in [21]. 

Ottman et al. [30, 31] were among the first to apply the MPPT scheme to piezo- 
electric harvesters using a DC/DC step-down converter. Although in [31], a dis- 
crete implementation is presented consuming much power, due to the 30 mW power 
extracted from the harvester, a stand-alone operation is claimed. The circuit is har- 
vesting power only for open circuit voltages of higher than 20 V, and for an open 
circuit voltage of 50 V, an efficiency of 0.7 is achieved. The CMOS realization pre- 
sented in [26] achieves a peak efficiency of 0.956 using the vibration tracker, and 
due to the startup capability, it is supposed to permit standalone operation. 
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The implementation of the non-linear processing techniques SSHI and SECE has 
been tackled several times in the past. The goal is always to derive the timing for 
the switch control signals from the piezoelectric voltage signal, in order to avoid the 
use of power consuming displacement sensors. Thereby, due to the fact that SSHI 
changes the shape of the piezoelectric voltage signal, the self-powered and stand- 
alone implementation of SSHI is more complex compared to SECE. Due to the power 
consumption of the control circuitry, the discrete-type self-powered SSHI imple- 
mentation presented in [25] achieves a performance which is significantly below the 
theoretical performance. Due to this, it overrides the performance of the standard 
rectifier only if the excitation level is high enough, i.e. if the output power is tending 
towards 1 mW. In contrast, Ramadass and Chandrakasan [36] present a self-powered 
CMOS integrated SSHI implementation which achieves a peak efficiency as high as 
0.87 for a harvested power level of 60 uW. However, a negative external voltage 
signal is required for proper operation, and the timing is not adaptive on varying 
excitation conditions, making it difficult to use this implementation in real world 
applications. 

In [51], a CMOS integrated implementation of SECE is presented. Providing that 
the buffer capacitor is precharged to a certain voltage level, the chip ensures self- 
powered operation at an acceptable efficiency of 0.7 and for a wide range of excitation 
frequencies. However, the timing of the switch control signals is fixed to one volt- 
age level and cannot be changed externally. In addition, the piezoelectric voltage is 
supposed to be lower than buffer voltage level, avoiding a voltage down-conversion 
and thus greatly limiting practical use. As is the case for the discrete-type realization 
of SSHI [25], the discrete-type implementation of SECE presented in [43] requires 
high output power levels in the order of milliwatts in order to achieve acceptable 
efficiencies. The third nonlinear extraction technique performing the piezoelectric 
voltage flipping is presented by Ramadass and Chandrakasan in [36] and by Sun et al. 
in [42]. By replacing the passive rectifier employed in [36] by an active counterpart, 
Sun’s implementation achieves a significantly higher efficiency than Ramadass’ chip. 
However, the timing in Sun's implementation is adopted from Ramadass' implemen- 
tation, so the problem of the invariable timing is the same. 

Lastly, a quite new nonlinear AC/DC conversion concept which supersedes a 
rectifier is presented by Kwon and Rincón-Mora in [21, 22]. Kwon makes use of 
a step-up DC/DC converter topology directly connected to the piezoelectric gen- 
erator and controls the involved switches such that the energy is directed into the 
buffer capacitor, similar to the SECE technique. The absence of a rectifier permits the 
minimum tolerated piezoelectric voltage to be as low as 0.35 V, making this imple- 
mentation on one hand appropriate for low input voltage and power levels. But on 
the other hand, in order to process negative piezoelectric voltages, an external —2 V 
voltage signal has to be applied separately, representing a significant limitation in 
practical use. 

Except for one chip presented in [26], a startup capability is missing in all of the 
impedance matchin implementations described above. This means that the buffer 
capacitor has to be precharged initially in order to make the circuits operable. In real 
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world applications where a continuous excitation of the harvester is not given, this 
would cause the circuit to fail if the buffer capacitor or the battery voltage drops 
below the minimum required supply voltage level. 


1.4 Goal of this Work and Major Achievements 


As discussed in Sect. 1.1, ofthe three kinetic transduction mechanisms, piezoelectric 
generators exhibit the highest practical maximum energy storage density. Moreover, 
they can be easily fabricated without requiring complex structures or numerous extra 
components. These points make the piezoelectric conversion mechanisms attractive 
for developing corresponding interface circuitry. 

The goal of the work presented in this book is to enhance the practicability of a 
CMOS integrated interface circuit for piezoelectric harvesters in real applications. In 
the opinion of the author, a practical interface circuit implementation should provide 
the following properties: 


1. Harvested power independency: The harvested power should be constant with 
regard to the load, 1.e. the buffer voltage. 

2. Power gain: Incorporating all circuit losses, the interface circuit should outper- 
form standard diode rectifiers. 

3. Operation independency: Proper operation should be independent of the excita- 
tion frequency, the piezoelectric voltage, and the buffer voltage. 

4. Autonomy: Every signal required for proper operation should be derived on-chip 
from either the piezoelectric voltage or the buffer voltage. This means that no 
external voltage reference signals should be needed. 

5. Self-powering/high efficiency: The control circuitry should consume as low power 
as possible, in order not to waste valuable power generated by the harvester. Any 
other losses, e.g. conduction losses or switching losses, should be also kept as 
low as possible. Thus, the control circuitry can be directly supplied by the energy 
stored in the buffer capacitor, and a high conversion efficiency can be achieved. 

6. Startup: The circuit should be capable of starting up from an uncharged buffer 
capacitor. 


By applying nonlinear extraction techniques, the output power of piezoelectric gen- 
erators can be greatly enhanced compared to what is achievable using full-wave 
rectifiers, even if there would be no losses within the rectifiers. Among the nonlinear 
extraction techniques, the SECE scheme is the only one providing a harvested power 
which is (ideally) independent of the load, i.e. the buffer voltage or the load resis- 
tance. Under certain conditions, i.e. for generators with a high electromechanical 
coupling, or for resonant excitation, the SSHIscheme can give more harvested power 
compared to SECE [16, 44, 45]. However, similar to rectifiers, the harvested power 
strongly depends on the load with SSHI.To conclude, SECE offers the highest level 
of practicability. Moreover, SECE is the technique with the highest "effectiveness", 
indicating the product of the harvesting efficiency(1.e. the harvested power compared 
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to the absolute maximum power) and the chip efficiency (i.e. the harvested power 
with regard to the input power) [10, ch. 6.3.1]. These are the reasons for concentrating 
on this scheme throughout this work. 

As discussed in Sect. 1.3.2.2, the state-of-the-art interface circuits performing 
SECE as presented in [51] and [43] are lacking of some of the requirements listed 
before, strongly limiting practicability. Hence, the main goal of this work is to present 
a pulsed synchronous charge extractor (PSCE) chip which fulfills all of the given 
requirements, as reported for the first time in [17]. The major achievements with 
respect to these requirements are as follows, whereas the performance values are 
taken from the PSCE chip v2 (cf. Sect. 6.4): 


e Due to the use of SECE, requirements 1 and 2 are fulfilled. As discussed in 
Sect. 6.4.4, due to the inevitable chip losses, the harvested power is not really 
constant over varying buffer voltages, but the dependence is much lower compared 
to standard rectifiers. The PSCE chip outperforms the standard rectifier under any 
excitation conditions by at least 20 % and up to 127 96, depending on the excitation 
conditions. 

e The operation independence as claimed in requirement 3 is analyzed in Sect. 6.4.5. 
Proper operation is achieved over a wide range of frequencies (50—2,000 Hz) and 
inductances of the external coil (100 uH — 1 H). Due to a capacitively-coupled 
peak detector and the use of high-voltage MOSFETS, piezoelectric voltages from 
1.3 to 20 V can be processed. Due to the use of a supply independent biasing 
as described in Sect. 5.7, timing is kept constant for buffer voltages in the range 
1.4-5 V. 

e Since all the signals required for proper operation are generated on-chip 
(see Chap. 5), autonomy is given according to requirement 4. 

e The control circuitry is designed to consume low power by activating the switch 

control on a per-need basis (see also Sect. 5.8.2). Based on simulations, the power 

consumption is as low as 2.5, W for a buffer voltage Vbuf = 2.5 V, as ana- 
lyzed in Sect. 6.4.3.2. In addition, a modified switching technique as described in 

Chap. 4 improves the efficiency by reducing the conduction losses. According to 

Sect. 6.4.4.1, the measured efficiency gain is in the range 0.06—0.16, depending on 

the used inductor and the buffer voltage, resulting in efficiency values up to 0.85. 

Hence, according to requirement 5, self-powering and high efficiency is achieved. 

Regarding requirement 6, the PSCE chip is capable of performing a startup from 

an uncharged buffer capacitor, as shown in Sect. 6.4.6. In case the buffer voltage is 

below the minimum supply voltage, the buffer capacitor is being charged passively 
via a bypass (see Sect. 5.5.2), until the minimum supply voltage of 1.4 V is reached. 


In addition to the chip realization, in Chap. 4, a method performing an approximation 
of the energy loss is presented. This allows an extensive analysis of the performance 
of the different switching techniques, without requiring long simulation times. 
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1.5 Organization of the Book 


The book is organized as follows. In Chap. 2, after the theoretical aspects of the 
piezoelectric effect have been briefly, it is described how this effect can be utilized 
for piezoelectric generators. Therefore, different converter configurations are inves- 
tigated. As a last point, the mechanical and electrical modeling of kinetic harvesters 
in general and piezoelectric harvesters in particular are analyzed. 

In Chap. 3, the resistor load, the full-wave rectifier load, and the SECE load are 
discussed analytically with regard to their output power performance. Finally, these 
different load circuits are compared with regard to the electromechanical coupling 
factor of the piezoelectric generator, stating which load circuit performs best under 
which conditions. 

Chapter 4 theoretically investigates the proposed PSCE circuit. The different 
switching techniques are discussed in detail and evaluated with regard to efficiency 
and especially to the efficiency gain of the improved techniques compared to the 
standard technique. 

In Chap. 5, the implementation of the PSCE chips v1 and v2 is described. First, 
the circuit blocks are illustrated on transistor level by means of schematics, and 
their function is explained using calculations and simulations. Afterwards, the chip 
implementations are shown on system level, and finally, a respective chip migrograph 
is given. 

At the beginning of Chap. 6, the vibration setup is shown. Then, the piezoelectric 
generators used as a power source for the measurement results are characterized, and 
their parameters required for the simulation model are extracted. After the demon- 
strator board using the PSCE chip v1 has been presented, the performance of the 
PSCE chip v2 is discussed in detail by means of measurements and simulations. 

Chapter 7 draws conclusions and gives an outlook on future work. 
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Chapter 2 
Piezoelectricity and Energy 
Harvester Modelling 


The direct piezoelectric effect is well-suited for harnessing environmental vibrations 
and to convert them into usable electrical energy. Coated on a cantilever beam, piezo- 
electric ceramics allows simple energy harvesting without additional mechanical 
structures. In addition, relatively small external excitations result in output voltages 
of several volts which do not require low-efficient startup mechanisms as known from 
thermoelectric generators. In order to give the reader a feeling about what is impor- 
tant for developing interface circuitry for piezoelectric energy harvesters, the basics 
of piezoelectricity and its usage in energy harevsting are presented in this chapter. 
After the physical conversion principle has been introduced roughly, its application 
in cantilever beam harvesters is explained. The last two sections are about the model- 
ing of kinetic, vibration based energy harvesters in general and piezoelectric energy 
harvesters in particular. Electrical equivalent circuits of the energy harvesters are 
essential for the simulation of entire energy harvesting systems composed of the 
mechanical harvester structure and the electrical interface circuitry on one hand, and 
on the other hand to understand their behavior by means of analytical calculations. 


2.1 Theoretical Background of the Piezoelectric Effect 


The theoretical treatment of the piezoelectric effect presented in this section is inten- 
tionally kept very superficial, because this book's focus is on the interface circuitry 
for piezoelectric energy harvesters, not the harvesters themselves. Nevertheless, the 
reader gets a basic introduction to this topic in order to be able to understand the 
challenges linked with energy extraction out of piezoelectric energy harvesters. A 
more thorough analysis of piezoelectricity can be found in [7]. 
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Fig. 2.1 Intrinsic piezoelectric effect in lead zirconate, showing the a crystallite above and below 
the Curie temperature Tc, where the charged zirconium or titanium ion moves relative to the center 
position 


2.1.1 History and Conversion Principle 


Certain solid materials like crystals and ceramics with no inversion symmetry [5] 
generate charges under an externally applied force. This so-called direct piezoelectric 
effect was initially discovered by the brothers Pierre Curie and Jacques Curie in the 
year 1880 and demonstrated using tourmaline, quartz, topaz, cane sugar and Rochelle 
salt (sodium potassium tartrate tetrahydrate). One year later, in 1881, Gabriel 
Lippmann found out that piezoelectric materials exhibit the reverse piezoelectric 
effect, meaning a deformation of the material when an electric field is applied [12]. 
Since then, intensive research has led to many applications like microbalances and 
scanning microscopy as examples for the direct piezoelectric effect, and motors and 
ultra-exact positioning as examples for the reverse piezoelectric effect. 

Lead zirconate titanate (PZT) is the most commonly used piezoelectric material 
at the moment. Since the generators used throughout this work are made of PZT, the 
explanation of the physical basics is limited to this material. The elements of PZT (Pb, 
O, Ti/Zr) crystallize as perovskite crystals with the general formula PbZrı-—xTix0O3. 
It has been shown that x ~ 0.5 maximizes the piezoelectric effect [22]. The piezo- 
electric effect is explained by means of Fig. 2.1. For temperatures below the Curie 
temperature (T < Tc), the charge center of the PZT crystallite moves away from 
the geometric center, leading to asymmetric electric dipoles with the characteris- 
tic rhomboedric/tetragonal structure. When the temperature is elevated above the 
Curie temperature (T > Tc), the crystallites are de-poled and the ceramic loses its 
piezoelectric characteristic. 

Groups of dipoles with the same alignment are called Weiss domains. Due to the 
random distribution of Weiss domains with different alignment inside the ceramics 
as depicted in Fig. 2.2a, there is no macroscopic piezoelectric effect. Similar to the 
magnetizing of a permanent magnet, the ceramics can be polarized by applying a 
strong electric field, resulting in an arrangement of the Weiss domains as shown in 
Fig. 2.2b. The ceramic is now ready to use and exhibits piezoelectric properties, i.e. 
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(a) (b) 





Fig. 2.2 Polarization of polycrystalline piezoelectric ceramic causes a the as-fired random domain 
polarity to align to b a net positive polarity, which c relaxes or ages over time 


charges are generated at the top and bottom electrodes, inducing a current when a 
load is connected, or changing its shape when a voltage is applied. Fig. 2.2c shows 
a ceramic where the Weiss domains are slightly misaligned compared to the case 
shown in Fig. 2.2b due to exceeding the mechanical, thermal or electrical limits, or 
due to aging. This results in a reduction of the piezoelectric effect. 


2.1.2 Constitutive Equations 


The piezoelectric effect can be described using the coupled equations [2, 13] 


=s".T+d'-E QA) 
+e’ -E, (2.2) 


where S is the mechanical strain, 7 is the mechanical stress, D is the electrical 
displacement (charge density), E is the electric field, s^ is the compliance under 
a zero or constant electrical field (indicated by the superscript E) and e^ is the 
dielectric permittivity under a zero or constant stress (indicated by the superscript T). 
d and d' are the matrices for the direct and the reverse piezoelectric effect, where the 
superscript means the transposed matrix. Thus, (2.1) and (2.2) describe the reverse 
and the direct piezoelectric effect, respectively. Without the coupling term d' - E, 
(2.1) is simply Hooke's law relating strain and stress. Likewise, (2.2) is simply the 
dielectric equation when the coupling term d - T is neglected. The detailed matrix 
notation of (2.1) and (2.2) for a material of the 6 mm crystal class such as PZT 
ceramic 1s shown in (2.3) and (2.4): 
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The subscripts 1, 2 and 3 are related to the x, y and z axis in the cartesian 
coordinate system. Due to symmetrical reasons, many matrix elements are zero, and 
only few are independent. Generally, D and E are vectors, i.e. tensors of rank 1, 
and the permittivity €^ is a tensor of rank 2. In principle, strain S and stress T are 
also tensors of rank 2, but due to the symmetry of the strain and stress tensors, the 
subscripts are conventionally re-labelled according to the so-called Voigt notation 
[13]: 11 —> 1; 22 > 2; 33 — 3; 23 —^ 4; 13 — 5; 12 — 6. Thus, they can 
be written as vectors or rank-1 tensors in (2.3) and (2.4). Whereas the first three 
elements are related to the force applied on the surfaces perpendicular to the x, y 
and z direction, respectively, the other three elements denote shear forces. Due to 
the Voigt notation, the order of the compliance tensor s^ can be reduced from 4 to 
2, resulting in the 6 x 6 matrix in (2.3). 


2.2 Piezoelectric Harvester Design Configurations 


Since this work focuses on piezoelectric material used in energy harvesting appli- 
cations, the following explanations refer to the case where the piezoelectric mate- 
rial is used as a sensor rather than an actor. As shown in the constitutive (2.4), 
there are only three independent elements within the piezoelectric coupling tensor d. 
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Fig. 2.3 Piezoelectric material operated in a 33 and b 31 mode 


A 





Fig. 2.4 Piezoelectric bimorph configured for a series and b parallel operation 


Whereas d15 is related to shear stress which is not practical for energy harvesting, it 
is worth describing d33 and d3; more detailed. 

Fig. 2.3 illustrates the two corresponding modes 33 and 31 in which piezoelectric 
material is generally used. The first number (3) indicates that the voltage is gener- 
ated along the z axis in both modes, i.e. the electrodes are attached to the surfaces 
perpendicular to the z axis. The second numer indicates the direction of the applied 
stress. In 33 mode, the stress acts in the same direction as the voltage appears. When 
the piezoelectric material is operated in 31 mode, the stress is applied along the x 
axis, whereas the voltage appears in the z axis. For 31 mode, the constitutive (2.3) 
and (2.4) can be simplified to 


(2.5) 


Sı — sh -Ti - dai: E3 
D3 = dq - T) + €3, - E3. 


The 31 mode is commonly used for energy harvesters, e.g. for bimorphs where two 
piezoelectric layers are attached together as shown in Fig. 2.4. When the cantilever 
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€ Piezo layer 
Tip mass m 


m NP 


Fig. 2.5 Cantilever structure converting the displacement z in the 3 direction into a elongation € 
in the 1 direction 


bends down, the surfaces of the layers where the electrodes are attached are stretched 
and compressed for the top and the bottom electrode, respectively. The opposite 
characteristic can be observed for a cantilever bending upwards. When the two layers 
are poled in the opposite direction, the layers are connected in series. Thus, compared 
to one layer, the voltage doubles, the capacitance is divided in half, and the current 
stays the same. This situation is also shown in Fig.2.4a. The parallel connection of 
the layers as illustrated in Fig. 2.4b is usually favored because the current and the 
capacitance are doubled whereas the voltage stays the same. Note that the type of 
layer connection only affects the voltage to current ratio, not the deliverable output 
power. 

Piezoelectric material is usually very stiff (or low-compliant) resulting in very 
high, impractical resonant frequencies when it is attached directly to the frame. 
Thus, the cantilever configuration shown in Fig. 2.5 is commonly used where a large 
deflection z in the 3 direction is transformed into a small deformation & in the 1 
direction. This leads to a higher compliance, up-scaling the strain and lowering the 
resonant frequencies to the 100—200 Hz range according to the frequency spectrum 
of practical ambient vibrations. Often, a tip mass is added to the end of the beam, 
further increasing the compliance and the strain and thus decreasing the resonant 
frequency. Although the d3; coefficient correlating to the cantilever configuration is 
lower than the d33 coefficient (d31 ~ 0.5d33 for PZT [17]), the 31 mode configuration 
is usually preferred in energy harvesting applications due to these advantages. 


2.3 Modeling of Kinetic Energy Harvesters 


In this section, the modeling of the general vibration-based kinetic energy harvester 
is discussed, mainly corresponding to [9]. Besides the dependency of the excitation 
frequency, it is described how the transducer affects the characteristics of the system, 
and what is the upper limit for the extractable output power. In addition, the electrical 
equivalent circuit is presented, allowing the computer-aided simulation. 
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Fig.2.6 Modelofa kinetic energy harvester structure composed of lumped elements with connected 
electrical harvesting circuit 
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2.3.1 The General Mechanical Model 
of Kinetic Harvesters 


Inetrial-based kinetic energy harvesters are modelled as second-order spring-mass- 
damper systems. The general model of kinetic energy harvesters was first developed 
by Williams and Yates [21]. Figure2.6 shows the general model of a kinetic energy 
harvester composed of lumped elements with a transducer connected to an electrical 
interface circuit [14, 18, 21]. The frame represents the harvester which consists of 
a seismic mass m suspended on a spring with the stiffness ks, generating a resonant 
spring-mass system. An external vibration which is assumed to exert a sinusoidal 
force on the harvester frame causes the frame to move harmonically with 


y(t) = Ysin (at), 


where y indicates the amplitude of the frame motion, and w denotes the angular 
vibration frequency. The relative movement of the seismic mass m with respect to 
the frame is given by 


z(t) — Zsin (ot + o), 


where Z indicates the amplitude of the mass motion, and ¢ points out that there is a 
phase difference between y(t) and z(t). 

Mechanical damping due to friction, air resistance etc. is represented by the 
damper d. The movement of the mass can be used to deform a transducer, con- 
verting the mechanical energy into usable electrical energy. There are many possible 
conversion mechanisms, e.g. the capacitive, the inductive or the piezoelectric (see 
Sect. 1.1). Due to electromechanical feedback, the transducer exerts a restoring force 
F, on the seismic mass if an interface circuit is connected, causing additional elec- 
trical damping de. The work which is done by the seismic mass against the restoring 
force is being converted into electricity. 
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2.3.2 Transfer Function 


In order to derive the governing equation of the general kinetic energy harvester, 
it is assumed that the mass of the excitation source is much larger than m, so the 
excitation source is not damped by the energy harvester and thus it is assumed to 
provide infinite power. The force balance is given as 


ma=mz+dz+ksz+ Fo, (2.6) 
where 
a(t) = (t) = —w* sin (wt) = Tsin (of) 
indicates the acceleration acting on the harvester frame. So, ma represents the exter- 


nal force which is being exerted on the harvester frame. If the restoring force is 
considered as a damping force Fe = dez, (2.6) can be rewritten as 


ma=mz+(d+d.)z+ksz. (2.7) 


In order to derive the transfer function, the Laplace transform of (2.7) is written 
as 


ms? y = ms^z4 (d + de) sz + ksz. (2.8) 


Using the dimensionless terms for the mechanical and the electrical damping, 


d d 
and ĉe = — (2.9) 
2mog 








fq = 


2mog 


where c, = ,/ ks denotes the natural frequency of the mechanical system, (2.8) can 
be rearranged to provide the transfer function 


KORE : 


y(s)  s2-F2os (£a + be) s + v2" any 


From (2.10), the displacement amplitude of the seismic mass can be derived as 


ls 


y 
i (2.11) 
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Fig. 2.7 Equivalent circuit of d 1/ks 
M A l Z m 
the kinetic harvester including Z 
electrical damping 
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2.3.3 Equivalent Circuit 


The mechanical equivalent model discussed before is important for understanding the 
physical behavior of piezoelectric harvesters. However, since the the main scope of 
this work is the development of an interface circuit, an accurate electrical equivalent 
model of the harvester has to be found that can be used in simulation software like 
(P)Spice [11, 20] or Spectre [11]. Therefore, generally, a mechanical force can be 
represented as an electrical voltage, and a mechanical velocity (i.e. the first derivative 
of the displacement) can be replaced by an electric current [19]: 


| 1 £z 
(2.12) 


|| 


VF 


Using this analogy, each term in (2.7) represents a voltage, and thus the electrical 
equivalent circuit as shown in Fig. 2.7 can be derived. The mass is represented as an 
inductor with the value m, the spring is modelled as a capacitor with the value E 
and the mechanical and the electrical dampings are described by the resistors d and 
de. The voltage source my = ma models the external vibration. See Sect. 2.4.3 for 
a more detailed derivation of the lumped element values. 


2.3.4 Output Power 


The power dissipated in the transducer, i.e. in the resistor de, can now be calculated 
as follows: 





P(w) = ~/-d, = do. (2.13) 


Using (2.11), equation (2.13) can be expanded as 
3 
m (2) c y C. 


E - Opi + [2a + to) "d 


P (c) = (2.14) 
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When the kinetic harvester is driven at resonance, i.e. @ = @, holds, (2.14) 
simplifies to 


P (o) n 


———— 7 2.15 
Ater m 


From (2.15), it can be calculated that the power dissipated in the transducer max- 
imizes for £; = £4, which means that the mechanical damping equals the electrical 
damping. For this condition, (2.15) can be finally written as Plim, representing the 
absolute maximum power extractable from a kinetic harvester: 

m mo; y^ 
lim — 16 tj , 
or, written in terms of the mechanical damping d = 2mq@,¢q and the acceleration 
amplitude @ = «^y = cy: 
(ma) 
Piim = ———. 2.16 
lim gd ( ) 


Sometimes, the quality factor Q is used instead of the damping ¢, where the two 


are related by Q = Te Thus, the quality factors due to mechanical and electrical 
damping can be expressed as 


1 1 
and Qe = —., (2.17) 


Ca 3 2. 


and the total quality factor combining mechanical and electrical damping can be 
written as 


|. ] " 1 
Q Qa Qe 
Combining (2.9) and (2.17), the quality factors can also be expressed as 


Qq = d and Qe = dde (2.18) 





Expression (2.16) points out that the maximum extractable output power only 
depends on the excitation force and the mechanical damping of the kinetic harvester. 
Since the maximum extractable output power is proportional to L, it seems desirable 
to try to keep the mechanical damping as low as possible in order to achieve a large 
output power. Decreasing the damping towards zero results in the maximum output 
power tending towards infinity, under the prerequisite that the input power provided 
by the excitation source is infinitely large. In terms of the quality factor, it seems 
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Fig. 2.8 Normalized maximum harvested power as a function of the excitation frequency, for 
different quality factors 


desirable to achieve large values, as the quality factor and the damping are reciprocal 
to each other, see (2.17). 

Fig. 2.8 illustrates the normalized power dissipated in the electrical damper, i.e. the 
harvested power, for different quality factors. As expected, the higher the quality fac- 
tor, the larger the maximum extractable power, but the sharper the peak. This means 
that the resonant frequency of the harvester has to perfectly match the excitation 
frequency, otherwise the extractable power drops significantly. 


2.4 Modeling of Piezoelectric Harvesters 


In our case, the piezoelectric harvesting is of special interest, so this transducer 
mechanism will be investigated in detail in the following. The conclusions made in 
Sect. 2.3 are therefore taken as a reference. 


2.4.1 Mechanical Modeling 


Let's assume the piezoelectric material shown in Fig. 2.9. The force F acting on 
the material causes an elongation e in the 1 direction. Using the relations 


g 
j 24 


Ee 2.19 
F (2.19) 


F E 
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Fig. 2.9 Beam of piezo material used as energy harvester in the 31 mode 


where q denotes the charge, the constitutive (2.5) can be rewritten in terms of the 
macroscopic variables F, €, V and Z instead of the local variables S, E, D and T [3]: 


F =kpe +I Vp 
L = l'é — Cp Vp [en 
where 
bh d2, X bl daib 
bam xul Ll PS 2:4 


In (2.21), kp denotes the stiffness ofthe piezoelectric material, C p is the piezoelec- 
tric output capacitance, and J" represents the generalized electromechanical coupling 
factor (GEMC) [16]. The first equation of (2.20) shows that the force F 1s composed 
of the “spring” force kpe which depends on the material stiffness and the coupling 
force Vp which depends on the voltage across the piezoelectric material. Due to 
the balance of forces, F can be considered as the restoring force F acting on the 
seismic mass, as defined in Fig. 2.6. 

Usually, piezoelectric materials are very stiff, so very high resonant frequencies 
would result if the piezo beam would be suspended directly within the harvester 
frame. Thus, the cantilever configuration shown in Fig. 2.5 is preferred where a large 
deflection z in the 3 direction leads to a small elongation e in the 1 direction. For this 
configuration, the constitutive (2.20) still hold [3], and hence they can be written as 


Fe = kpz + I'Vp 
and thus (2.6) can be written as 
ma = mz + dz 4-kz - I'Vp, (2.23) 


where k = kp + ks is the sum of the stiffnesses of the piezoelectric and the mechan- 
ical structure. This can be done under the assumption that a perfect bond exists 
between the beam and the piezoelectric material [8]. 
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Finally, the spring mass damper system shown in Fig. 2.6 can be modelled by the 
differential equations 
ma=mz+dz+kz+I Vp 
[Parente 00 224 


The energy balance of the vibration harvester system can be derived by multiplying 
(2.23) with the mass velocity z(t) and integrating over the time variable f [6]: 


1 1 
/ maz dt = zm? + J dż? dt + ake" + / I Vpz dt (2.25) 
Wem — 00 eee” 
total kinetic damping elastic electrical 
where 
I 2 
I’ Vpzdt = "ad + | Vpl dt. (2.26) 


Expression (2.25) states that the energy injected into the system is composed 
of the kinetic energy, the mechanical damping losses, the elastic energy and the 
energy converted into electrical energy. According to (2.26), the energy converted into 
electrical energy is separated into the energy stored on the piezoelectric capacitance 
and the energy absorbed by the electrical load. The latter energy is the part which is 
actually being harvested. 


2.4.2 Electromechanical Coupling and Damping 


The piezoelectric coupling factor is an important parameter because it is a mea- 
sure for the part of injected mechanical energy which is converted into electrical 
energy. High coupling means that a large amount of energy can be harvested from 
a given environmental mechanical energy, i.e. the last term in (2.25) is large. Vice 
versa, the same mechanical energy leads to only little useable eletrical energy for 
a lowly coupled harvester. Since the actuator piezoelectric effect is present in the 
reverse direction, an electrical load drawing power out of the harvester induces a 
feedback on the deflection of the piezoelectric beam which induces electrical damp- 
ing, corresponding to d, as defined in Sect. 2.3.1. This means that a highly coupled 
piezoelectric harvester is damped more intensely than a lowly coupled one. Hence, 
the coupling factor is of outstanding importance because it determines whether the 
SECE technique discussed in this book extracts more power out of a piezoelectric 
harvester compared to a standard bridge rectifier. 
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For the piezoelectric beam shown in Fig. 2.9, according to the IEEE Standard on 
Piezoelectricity [13], the squared coupling factor is given by 


2 
h = BL. (2.27) 
234351] 


Note that ks. exclusively depends on piezoelectric material properties. The mea- 
sured squared material coupling factor KS. has been reported to be up to 0.7 [3]. In 
order to link the material coupling factor i and the generalized electromechanical 
coupling factor (GEMC) I” which depends on the geometry of the used cantilever 
as defined in (2.21), (2.27) can be expressed as 


2 
d =—— 
31 7 . 
kpCp 
In order to describe the total harvester structure (e.g. the resonator geometry) instead 
of the piezoelectric layer properties, the effective coupling factor can be used [13]: 


2 _ 
— (2.28) 
w 


kert = 
where Woc and cx. denote the angular resonance frequency of the open-circuited and 
short-circuited structure, respectively. The effective coupling factor is a convenient 
parameter allowing to compare different harvesting structures using any piezoelectric 
material. According to [6, 15], realistic values of koe are in the range of 0.01—0.1. 
Due to the electromechanical feedback, two resonant frequencies exist for the 
piezoelectric harvester, depending on the electrical loading. The fundamental reso- 
nance frequency is determined for the case that the piezoelectric harvester terminals 


are short-circuited: 
Lk 
Msc — mm" (2.29) 
m 


The anti-resonance frequency is higher than the fundamental resonance frequency 
and shows at open circuit configuration: 


Woe = Ggc4/ 1 + ke. 


To avoid confusion, it is emphasized that k denotes the total stiffness composed 
of the mechanical structure and the piezoelectric material, whereas ka indicates the 
squared effective electromechanical coupling factor. For both the fundamental and 
anti-resonance frequency, the electromechanical damping exerted by the respective 
electrical load (i.e. Rg = O0 for the fundamental frequency and Ry; = oo for the 
anti-resonance frequency) is zero [4]. 
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Fig.2.10 Uncoupled equivalent circuit of the piezoelectric harvester a neglecting and b considering 
the dielectric losses 


The following equation correlates the GEMC to the effective electromechanical 
coupling factor [16, 18]: 


For the equivalent circuits presented in the following, the GEMC J” will be used 
in order to describe the coupling phenomenon. 


2.4.3 Equivalent Circuits 


As a first approach, the second equation of (2.24) can be converted into the uncoupled 
electrical equivalent circuit indexUncoupled electrical equivalent circuit shown in 
Fig. 2.10a [1, 18]. The sinusoiodal current source represents the vibrating beam with 
Ip = Iz, and the piezoelectric output capacitance is modelled by Cp. By means 
of (2.11), the amplitude of the current source Íp can be expressed in terms of the 
external acceleration à = Ş and several generator parameters as follows: 


~ rm 
IP oc = a+ —— —, (2.31) 
d 
assuming excitation at resonance (i.e. @ = œn) and no load (i.e. £; = 0). For 


the Synchronous Electric Charge Extraction (SECE) technique which is the basic 
principle of the proposed interface circuit (see Sect. 3.3), in a first approximation it 
can be assumed that ¢e = Z4, and hence the following equation holds at resonance: 


T, "E (2.32) 
P,SECE — a od . . 


The voltage across the piezoelectric terminals is indicated by Vp . As depicted 
in Fig. 2.10b, sometimes the resistor Rpar is added to consider the dielectric losses 
since the piezoelectric capacitance is not a perfect insulator. Since Rpar is usually 
very high (>10 M2), it is neglected throughout this book. This equivalent circuit 
completely neglects the mechanical resonator structure, but it is sufficient for the 
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Fig. 2.11 Coupled electromechanical equivalent circuit of the piezoelectric harvester [10, 16] (The 
controlled sources are sometimes replaced by a transformer.) 


design process of interface circuits, where it is irrelevant which external vibration 
force generates the actual voltage across the piezoelectric terminals. 

However, if the feedback of the interface circuit on the piezoelectric harvester or an 
accurate prediction of the output power for a given vibration acceleration is important, 
a more sophisticated electrical equivalent circuit has to be found. According to (2.12), 
the mechanical velocity z can be represented by an electrical current: 


Z = E. 


and hence the differential (2.24) can be rewritten as: 
ma = mIn +dIn +k f Imdt +I Vp (2.33) 
I =TI,—CpVp 


Note that each term in the first equation of (2.33) can now be represented as an 
electrical voltage instead of a mechanical force, as stated in (2.12). Thus, the coupled 
electrical equivalent circuit composed of lumped elements, as shown in Fig. 2.11, can 
be derived [10, 16]. The electrical domain is a representation of the second equation 
of Fig. 2.33, so the value of the capacitor can be directly defined as C p. An additional 
parallel resistor modelling the leakage losses is omitted since its resistance is usually 
very high (>10 M2). 

A little more effort has to be made in order to determine the lumped element para- 
meters of the electrical domain. Therefore, the method of equating the coefficients 
between the first equation of (2.33) and the Kirchhoff Voltage Law (KVL) of the 
mechanical domain is applied: 


; I 
m 
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Fig. 2.12 Simplified electromechanical equivalent circuit of the piezoelectric harvester, where the 
coupling factor is considered in the lumped element parameters 


As a result, the mass is modeled by an inductor with inductance L,, — m, the 
stiffness of the piezoelectric beam is represented by a capacitor with capacitance 
Cm = 1/k, the parasitic damping is modeled by a resistor with resistance Rm = d, 
and the alternating input force is modeled as a voltage source V, = my = ma. 
Controlled voltage and current sources using the GEMC (see Sect. 2.4.1) are used 
to couple the mechanical and the electrical domain, alternatively a transformer with 
the winding ratio 1 : J” often found in literature has equal functionality. 

Finally, the controlled sources/the transformer can be eliminated as shown in 
Fig. 2.12, simplifying mathematical calculation. Therefore, the first equation of (2.33) 
is multiplied by 45: 
oe eee | (2.34) 

I =I Im— CpVp 


Now, introducing the new symbol 
Imc — I Im, 


and performing the method of equating the coefficients between (2.34) and the KVL 
of the simplified equivalent circuit: 





. 1 
Vmc — LmcImc T Rmce Imc t C 


mc 


J eate v 


the lumped element parameters of the simplified equivalent circuit are given by 


I 
Rme = —5  Cme = T (2.35) 


Both the equivalent circuits shown in Figs.2.11 and 2.12 accurately model the 
electromechanical feedback, i.e. the resonance frequency shift and the piezoelectric 
voltage change 1f an electrical load is connected (see also Sect. 2.4.2). However, the 
simulation time using these equivalent circuits can be long since a startup oscilla- 
tion takes place until the steady state is found. Since the accurate behavior of the 
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piezoelectric harvester in terms of feedback is rather unimportant for the design of 
the interface circuitry, the basic model shown in Fig. 2.10 can be used during the 
design process. The simulation time using the simple model is much shorter than 
using the complex model since no steady state has to be found. Hereby, the excitation 
is modeled by a sinusoidal current source which represents the vibrating beam. In 
most cases, it is sufficient to set its current amplitude such that a desired piezoelec- 
tric voltage is generated. For the design of the proposed Pulsed Synchronous Charge 
Extractor (PSCE) interface circuit, the current amplitude can be calculated using 
(2.32). The piezoelectric capacitance should have a somehow realistic value since it 
determines how much energy can be extracted from the harvester during one transfer 
cycle. 

After the successful chip design and fabrication, it may be interesting to know 
how the chip performs connected to different piezoelectric harvesters, especially in 
terms of the power that can be extracted. For that reason, it is important to have an 
accurate model that behaves similar to a real piezoelectric harvester. 


2.4.4 Maximum Output Power 


The main figure of merit of interface circuits for any kind of generator is the power 
which they are able to extract. Therefore, it is helpful to know the absolute maximum 
power which can be extracted from the harvester. Although this is only a theoretical 
value which cannot be achieved in real applications, it is nevertheless a means for 
classifying the actual extracted power. 

Basically, a generator with internal impedance Z outputs a maximum power if 
it is terminated by a load impedance equal to its conjugate complex impedance Z. 
Using the simplified model from Fig. 2.12, the internal impedance of a piezoelectric 
harvester is given as 


Z = 
jæCp + 


1 
KRmc + C + JoLmc 


By substituting the imaginary j by — j, the conjugate complex impedance can be 


written as 1 
VAP 


—joCp + 


1 


1 . 
mc — JoCa — JO Lmc 


Thus, in order to get an electrical equivalent circuit representing Z, each capacitor 
within Z has to be substituted by an inductor and vice versa, leading to the equivalent 
circuit shown in Fig. 2.13. The conjugate complex pair Cp — Lp represents an infi- 
nitely high impedance, and the other pairs C — L and L — C create a zero impedance, 
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Fig. 2.13 Piezoelectric model with conjugate complex load 
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Fig. 2.14 Simplified piezoelectric model with conjugate complex load 


if the values of the load impedance components are set as follows: 


= l = 1 — ] 
Lp = ——— = ———— Cm = ———. 2.36 
P w2C p mc wC me mc ( ) 


Qus 

In this case, the circuit from Fig. 2.13 reduces to the circuit depicted in Fig. 2.14, 
which represents a simple resistive voltage divider. The root mean square power 
dissipated at the load resistor represents the maximum extractable power, and it can 
be calculated as 





cl 

y2 

Pim = — 5 = A Y2 / (2.37) 
Kme Rinc 


where Vp rms indicates the root mean square voltage at the load resistor, and Vac 
denotes the amplitude of the voltage source Vme. For Vme = ma/T and Rme = d/I? 
as given in (2.35), (2.37) can be finally expressed as 


Pim = ——. (2.38) 


Note that the derivation of Piim shown above and the derivation based on the 
mechanics of kinetic harvesters explored in Sect. 2.3.4 lead to identical (2.38) and 
(2.16). 
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Chapter 3 
Analysis of Different Interface Circuits 


As described in the introduction, there are many ways to extract power from a 
piezoelectric harvester. On one hand, the diode rectifier is a simple interface circuit 
applicable to any power generator with an alternating output. Since the diodes are 
passive elements, there is no power supply necessary, which can be advantageous 
when the generated power is low. On the other hand, active power processing circuits 
like the formerly presented Synchronous Electric Charge Extraction (SECE) tech- 
nique can increase the extracted power and decouple the load from the harvester. 

In the following, the basic theory behind the resistor load, the commonly used 
full-wave rectifier load, and the SECE load is described. At the end of this chapter, 
the characteristics of these different interface circuits are compared with regard to the 
electromechanical coupling and the excitation modes at resonance and off-resonance. 
Whereas due to its linearity, the analysis of the purely resistive load can be done 
straightforward using simple Laplace transformations, the analysis of the full-wave 
rectifier load and the SECE load requires complex math to get general closed-form 
expressions for the output power. Since the fundamental analysis of the different 
interface circuits is not the main goal of this work, for the two latter load circuits, it 
is distinguished between excitation at resonance and off-resonance. 


3.1 Resistor Load 


The resistor load as shown in Fig. 3.1 represents the simplest "interface circuit" for 
piezoelectric harvesters. Due to the absence of a rectifying stage, the voltage and 
the current at the load resistor is alternating. In addition, there is no storage element 
where the harvested charge can be accumulated for later usage. Thus, the resistor load 
is inappropriate for real energy harvesting applications. Instead, the resistor load is 
very often used to characterize the piezoelectric devices with respect to their output 
power cababilities. In this book, it is used to verify the correctness of the simulation 
model compared to a real piezoelectric device. 
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Fig. 3.1 Piezoelectric I 
harvester with resistor load 


3.1.1 Calculation of Harvested Power 


In the following, an expression for the output power using the resistor load is derived. 
As shown in Fig. 3.1, the output current J = a is flowing through the load resistor 
RL, producing the voltage Vp. Since the whole system is linear, the constitutive 
equations in (2.24) can be transformed into the frequency domain using Laplace 
transformation, resulting in 


ma = (ms? +ds +k) z e T V» 


V : (3.1) 
E = ['sz — SCpVp 
Rr 


which is equivalent to the matrix reprentation 


ma z \ (ms +ds+k r 
(5)-()( I's Loc): (3.2) 


where s = jæ denotes the complex frequency variable. Using Cramer's rule, the 
system of equations (3.2) can be solved for the piezoelectric voltage Vp: 


mal's 


ee, 3:9 
(sCp + $) (ms? + ds + k) + T?s e» 


Vp = 


Using (3.3), the RMS power dissipated in Rz can now be calculated as follows: 


7 (vel /v2) _ VpVp 








ee oe (3.4) 
_ (ma)? o^ Rir? 

2. (k— (m 4- Ced Ri) o?) + o? (d + Rz (F2 + Cp (k —mo2)))- 

(3.5) 


where V p is the complex conjugate of the piezoelectric voltage. Note that (3.5) 
models the resonance characteristics of the harvester properly, hence it is not required 
to distinguish between the excitation at resonance and off-resonance. 
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The optimal load resistance which gives maximum output power is 


d?a? + (k — mo?) 


= E E S r (3.6) 
C? dot + (T?w + wCp (k — wm) ) 


Rres,opt = 
and the optimal squared coupling factor which gives maximum output power is 


res,opt — 


1 
j2 — (1+ C$ R20?) (da? + (k — mo?) . (3.7) 
wR L 
If both the load resistance and the squared coupling factor are chosen optimally, 


the maximum power 


(ma)? 
Sd 





(3.8) 


Pres,max = Piim = 


can be achieved. 


3.1.2 Discussion of the Harvested Power 


Figure 3.2a shows the 3D surface plot of the harvested power as a function of the load 
resistor and the squared coupling coefficient, normalized to the absolute maximum 
power Pj. First of all, it can be observed that for a given load resistance, optimum 
power can be extracted for one single coupling factor I ae as calculated in (3.7). 
This phenomenon is described in detail in Sect. 3.4.1. Foraloadresistance in the order 
of 1 kQ, the coupling factor must be very high in order to harvest optimum power. 
In this regime, the optimum power is slightly below Plim. If the load resistance has 


a higher value in the order of 100 kQ2, a much smaller coupling factor is sufficient 
to extract optimum power. At a certain pair (Rr, r°) = (Rres,opt a) the 


optimum power equals the absolute maximum power Pj. For a load resistance 
higher than Rres,opt, it is not possible to extract the maximum power, even for I” 2 — 
i 

As illustrated in Fig. 3.2b, (3.5) considers the resonance characteristics properly. 
The resonant frequency is somewhere around 175 Hz, which will be verified by sim- 
ulation and measurement in Sect. 6.2.4.1. The fact that the optimum harvested power 
is lower than Pj indicates that the squared coupling factor does not correspond to 
the optimum one. 

For a given coupling factor value, there exists exactly one single value Rg = 
Rres,opt as provided in (3.6), at which the harvested power exhibits a maximum. As 
opposed to the exact calculation shown in Sect. 3.1.1, the approximation made in [2, 
Sec. II-A-1] leads to the appearance of two different optimal load resistances for a 
high coupling factor. This characteristic comes from the assumption that the external 
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Fig. 3.2 a Function of Ry, and I°’. b Function of Rz and f normalized harvested power using the 
resistor load (harvester values taken from 6.1) 


vibration force ma and the velocity z are always in phase, i.e. the device is always 
at resonance, independently of the coupling and the load resistance. In fact, this is 
not true, since the resonant frequency is a function of the coupling factor and the 
attached electrical load, as discussed in Sect. 2.4.2. 


3.2 Full-Wave Rectifier with Capacitor 


Full wave rectifiers represent a simple way to convert the AC type output voltage 
of vibration based energy harvesters into a DC type voltage which can be used to 
supply common electronics. Figure 3.3 depicts the piezoelectric harvester with the 
full-wave rectifier interface and the load composed of a buffer capacitor and a parallel 
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Fig. 3.3 Piezoelectric harvester with a full-wave rectifier composed of diodes and load composed 
of a buffer capacitor Cpuf and a resistor Rr, 
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Fig. 3.4 Typical waveforms of the output voltage of the piezoelectric harvester, Vp, the movement 
of the mass, z, and the current flowing into the buffer capacitor, [pur 


resistor. For the following analysis, the buffer capacitor Chur is assumed to be very 
large so that the voltage Vbuf can be seen constant over several excitation periods. The 
diodes are supposed to be ideal, i.e. they show no forward voltage drop. The typical 
waveforms of Vp, the output voltage of the harvester, together with the movement of 
the seismic mass, z, and the current flowing into the buffer capacitor, /pyr, are given 
in Fig. 3.4. The diodes conduct only if Vp reaches Vbuf; since Vpuf is fixed, Vp is 
clipped during power transfer to the buffer capacitor. When the mass displacement 
starts to decrease after its maximum Z, Vp falls proportionally to the mass velocity z 
until — Vpu¢ is reached, according to the second equation of (2.22) for 7 = 0. During 
the time where — Vbuf < Vp < Viuf, the diodes are blocking and no energy transfer 
occurs, meaning that the piezoelectric harvester is driven in open circuit. 
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3.2.1 Calculation of Harvested Power 


In the following, the power which can be harvested using the rectifier-capacitor load 
is discussed, as it was originally derived by Guyomar [2]. The charge accumulated 
on the buffer capacitor during two mass displacement extrema can be calculated by 
the integral over the current puf. In steady state, this charge is equal to the charge 
removed by the load resistor. Thus, the charge balance can be written as 


Hie Vour T 
| ET eee (3.9) 
0 





where = = i indicates one half-period of the excitation. If the current /py¢ flowing 
into the buffer capacitor is substituted by the second expression of (2.22), (3.9) can 


be rewritten as 
Vbuf a 





T/2 
| (Cz — Cr Vp) dt = LP 


leading to an expression for the buffer capacitor voltage: 


Vi ar, (3.10) 
= ———— MZ, x 
as R,Cpo + 3 : 
Using (3.10), the extracted power can be calculated as 
2 
Prect = “out 
Rr 
Rir- T 
= — OR ex. (3.11) 


In contrast to the purely resistive load, the recitifer represents a non-linear load 
which requires complex math in order to get a closed-form expression for the dis- 
placement amplitude 7. Therefore, in order to simplify things, a differentiation is 
made between excitation at resonance and at off-resonance, as originally shown in 
[2]. A more extensive analysis deriving closed-form equations as a function of the 
excitation frequencyis given in [1]. 

If the harvester is driven at off-resonance, the resulting output voltage Vp is low, 
and the same is true for the feedback force "Vp. So, the beam displacement is not 
strongly affected by the load, and hence it can be assumed to be constant. As a result, 
the output power at off-resonance (superscript 0^) can be written as 


pes Pacis (5.12) 


rect *— 


Under the assumption that the harvester is driven at resonance, i.e. the driving 
force ma and the velocity z are in phase, the mass displacement amplitude Z can 
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be calculated using the energy balance expression given in (2.25). Since the kinetic 
energy at the instants tı = T/2 and t? = 3T /2 (see Fig. 3.4) and the variation of the 
elastic energy between the same instants are both zero, the provided energy is the 
sum of the mechanical losses and the harvested energy: 


T /2 ig e V2. T 
| maidt = | dz? dt + 2L. (3.13) 
0 0 Ry 2 


From (3.13), the mass displacement amplitude can be calculated under the 
assumption that the displacement remains sinusoidal, which is a good approximation 
for structures with low viscious losses [2]: 

" ma 
ee (3.14) 


2Rrol2 
d |. 2RLo01 ^ — 
w + (RL Cpo+mn/2)? 


where à denotes the amplitude of the external excitation (t). Using (3.11) and 
(3.14), the harvested power at resonance (superscript R) can be expressed as 


2 
Ri r? ma 
R L 
Prect = (R C mr FP ZR r? (3.15) 
Hp 2 (Ry Cpoz/2) 


For the off-resonant excitation, the optimal load resistor at which maximum power 
is extracted can be calculated as 


OR "m 
R rect, opt — Cpa’ (3.16) 


and the maximum power is given by 


OR J ~ 
P on 2 m : (3.17) 





In the following, the resonant excitation is considered. When the expression for 
the output power (3.15) is investigated for extrema, it becomes evident that there two 
solutions possible for Rz: 


R OR 7 
Frect,opt — Rrect, opt m 2650 (3.18) 


and 


RR 2r? — n Cpdo +2r yT? — 5 Cpdo 


rect,opt, 1,2 — Cdo? (3.19) 
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Hence, there are two corresponding values of the maximum extractable power: 








R I (ma) 
Prectmax = anCiO], mW (3.20) 
(a T a) 
and . -- 
Prect,max = Pim = gd (3.21) 


Since (3.20) is dependent on the coupling factor, it is worth calculating the opti- 
mum value of the coupling factor for which the first possible solution of the extracted 
power becomes maximum: 


d(m + 2CpRro)? 

2 

 rect,opt = SRI. —. (3.22) 
Combining (3.18) and (3.22) provides an expression for the border between the 

two regimes, which is no more a function of Rz: 


r2 


rect,border 


= m doC p. (3.23) 


3.2.2 Discussion of the Harvested Power 


Figure 3.5 illustrates the characteristics of the harvested power as a function of 
the squared coupling coefficient and the load resistance for excitation at resonance. 
For the following analysis, it must be distinguished which one of the variable pair 
(R Lad A is kept constant. For a given load resistance, there exists exactly one squared 
coupling factor calculated in (3.22) at which the optimum power can be extracted. 
A physical explanation of this effect is given in Sect. 3.4.1. 

For a given coupling factor, there are two different regimes. If the electromechani- 
cal coupling is small, i.e. I" 2 < P Ade there exists one optimum for the extractable 
power at one resistance, as given in (3.20) and (3.18). For high electromechanical 
coupling, i.e. T? > I rein the resistance given in (3.18) corresponds to a local 
minimum of the harvested power, and one can distinguish between two possible 
optima for the extractable power at two resistances, as provided in (3.21) and (3.19). 
According to [2], this particular behavior with two optimal load resistances is char- 
acteristic for strongly coupled systems, such as piezoelectric transformers [3]. 

From Fig. 3.5, it becomes obvious that the absolute maximum power Pj can be 
harvested if Rz equals the optimum load resistance, but only for the regime of high 


coupling. For low coupling, it is not possible to achieve Piim. 
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Fig. 3.5 Normalized harvested power using the rectifier-capacitor load, assuming excitation at 
resonance (harvester values taken from Table 6.1) 


3.3 Synchronous Electric Charge Extraction 


Figure 3.6 illustrates the energy harvesting system using the Synchronous Electric 
Charge Extraction (SECE) interface. It consists of a rectifier, three switches S1, S2, 
S3 and an inductor L. The NVC flips the negative half-wave of the sinusoidal input 
signal to the positive domain. The switches are used to store the energy drawn out of 
the piezoelectric harvester temporarily in the inductor L, before it is released into the 
buffer capacitor. By implementing this intermediate step, the SECE interface is able 
to completely extract the energy stored on the piezoelectric capacitance Cp during 
each cycle. 

Figure 3.7 shows the corresponding waveforms of the output voltage of the piezo- 
electric harvester, the mass displacement, the inductor current and the digital signals 
controlling the switches. The SECE interface operates only at specific instants, dur- 
ing short time intervals compared to the excitation period, where power is extracted 
from the piezoelectric harvester. These instants are triggered by the maxima and 
minima of the mass displacement which are in phase with the maxima and minima 


SECE Interface 





Fig. 3.6 Piezoelectric harvester with SECE interface and load 
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Fig. 3.7 Typical waveforms of the output voltage of the piezoelectric harvester, Vp, the movement 
of the mass, z, the current flowing through the inductor, Iz, and the digital signals controlling the 
switches 


of the voltage between the piezoelectric harvester terminals. Apart from these power 
extraction cycles, the SECE interface keeps the harvester in open circuit, which 
means that no current is flowing into the SECE interface. 

When the harvester voltage Vp reaches a maximum Vp, the energy accumulated 
in the piezoelectric capacitance is also maximum. At that instant, the energy is 
instantaneously transferred into the inductor by turning on switch S, resulting in 
a fast decay of Vp. When Vp crosses 0 V, S, is turned off, and S5 and 55 are 
turned on, releasing the energy into the buffer capacitor. After the energy transfer 
process has completed, the SECE interface is shut down, again turning into infinite 
input impedance. A more detailed discussion of the transfer process using different 
switching techniques based on SECE is given in Chap. 4. 


3.3.1 Calculation of Harvested Power 


In the following, the output power which can be achieved using the SECE interface 
is derived according to [4]. Therefore, it is assumed that there are no losses inside 
the SECE interface. Since no current is flowing for a large part of time, the second 
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equation of (2.22) can be integrated using 7 = 0, which results in an expression for 
the amplitude of the harvester output voltage: 


P I. 
Vp —2——z. (3.24) 
Cp 


The energy stored on the piezoelectric capacitance is extracted once completely 
during a half-period of the excitation. So, using (3.24), the average extracted power 
can be calculated as: 


1 v2 2 
;CpV I' ^0 32 
Psgcg = “p= 2. (3.25) 
2 ids 


Similarly to the rectifier load, since the displacement remains approximately con- 
stant for excitation at off-resonance, the corresponding expression for the extracted 
power is given by 


POR. = Psgcg. (3.26) 


The energy balance is used to derive the displacement amplitude Z, when the piezo- 
electric harvester is driven at resonance. Since the power extraction happens twice 
per excitation period, the energy balance after a half-period is taken into account: 


P 15 r2 
| mazdt = / d a2. (3.27) 
0 0 Cp 


Again, the given energy balance implies that the kinetic and the elastic energy vanish. 
Assuming that the displacement remains sinusoidal, the displacement amplitude can 
be calculated as 


~ 


ma 


AI? ’ 
dæ + nC 


z= 


(3.28) 


where à denotes the amplitude of the external excitation ÿ(t). Thus, the extracted 
power for excitation at resonance can be written as 


2 
PÈ =2 re 9 (3.29) 
SECE — 1Cp p " ADT à ' 
TEP 


In order to find the optimum coupling factor, (3.29) has to investigated for extrema. 
It turns out that the maximum power 


(ma)? 
PsECE,max = Plim = E (3.30) 


can be achieved for 
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Fig. 3.8 Normalized harvested power using the SECE load, assuming excitation at resonance 
(values taken from Table 6.1) 


I SECE,opt m 4 : (3.51) 


3.3.2 Discussion of the Harvested Power 


Figure 3.8 shows the harvested power with connected SECE interface, assuming 
exctitation at resonance. As expected from (3.29), the harvested power is independent 
of the load, under exclusion of any parasitic losses. This is due to the nonlinear 
switching which decouples the piezoelectric harvester from the load. This sets the 
SECE technique apart from the other extraction techniques such as SSHI [5] or 
the previously described diode rectifier, because hence it is possible to keep the 
harvesting system in an optimum operating point independently of the load. 

Another important observation is the fact that there exists a value of the squared 
coupling factor at which an absolute maximum of the harvested power is achieved. A 
physical justification for this effect which has already been observed for the resistor 
and the rectifier-capacitor load is given in Sect. 3.4.1. For the SECE technique, the 
optimum coupling factor is given in (3.31). 


3.4 Comparison Under Coupling Considerations 


In the following, the interface circuits discussed up to now are compared, taking 
the electromechanical coupling into account. The basics of the electromechanical 
coupling phenomenon present in piezoelectric energy harvesters are explained in 
Sect. 2.4.2. 
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Fig. 3.9 Comparison of the harvested power with connected rectifier and SECE interface as a 
function of the squared electromechanical coupling factor, assuming operation at resonance (Ry, = 
10 kQ2, harvester values taken from Table 6.1) 


3.4.1 High/Low Coupling Factor 


Figure 3.9 shows a comparison of the harvested power of the resistance load, the 
rectifier-capacitor load and the SECE load as a function of the coupling factor. A 
load resistance value of Rg = 10k has been chosen arbitrarily in order to ensure 
comparable conditions. The first important observation is the fact that the optimum 
extracted power using the rectifier-capacitor load and the SECE load equals the 
absolute maximum power Piim, as predicted in the theoretical analysis, whereas the 
optimum power with the resistance load is below Pj. This is due to the fact that using 
the resistor load, there exists only one dedicated pair (R Bs E" » = (a ia 
for which the optimum output power can be achieved. Obviously, the chosen load 
resistance Rz = 10k{2 does not correspond to the optimal value. 

Secondly, it is evident that there exists an optimum squared coupling factor for 
each interface circuit, as given in (3.7), (3.22) and (3.31). This can be explained by 
means of the damping of the beam displacement due to the electrical load as discussed 
in Sect. 2.4.2. Usually, for a higher coupling factor, more mechanical energy is 
converted into electrical energy. This is true at least for the case that the coupling is 
small since in this regime, the electrical load has only a small damping influence on 
the beam displacement, because the backward coupling force amplitude represented 
by I'Vp in (2.24) is low compared to the external vibration force amplitude ma. 
Above a certain coupling factor, the electrical load affects the beam displacement in 
a stronger way so that an increase of the coupling factor leads to a reduction of the 
extractable output power. Hence, all the curves show an optimum output power at a 
certain coupling factor. 

Now, the question arises why the optimum squared coupling coefficient of the 
SECE interface is much smaller compared to the purely resitive and the rectifier- 
capacitor load which both exhibit maximum power extraction at similar coupling 
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coefficients. There exists a crossing point between Psgce and Prect at 


1 IT Cpo 
DR m 144 zy E (x t 2CpRro), (3.32) 


indicating the coupling factor where the extracted power using the SECE load and 
the rectifier-capacitor load are equal. Whereas the SECE load performs better than 
the rectifier-capacitor load for F? < I4... the rectifier load outperforms the SECE 
load for r? > b This non-intuitively understandable phenomenon has been 
discussed widely in literature [4, 6—8]: As described above, when the coupling factor 
is low, the mechanical structure of the piezoelectric harvester is damped only very 
weakly by any electrical load. Hence, in this regime, a higher value of I? leads to 
a higher harvested power, because more power is converted from the mechanical 
into the electrical domain. In the regime of high I? values, the electromechanical 
feedback can no longer be neglected, i.e. the backward coupling force amplitude 
I Vp becomes dominant compared to the external vibration force amplitude ma. 
Thus, if I? increases, the piezoelectric beam is being more and more damped, and 
above some turnover point, this results in a decreasing harvested power. For the 
SECE load, the turnover point is obviously reached for lower I"? values than for 
the rectifier load, which is due to the fact that the SECE load enforces higher output 
voltages and thus a stronger backward coupling force. 


3.4.2 Excitation at Resonance/Off-Resonance 


Figure 3.10 shows a comparison of the harvested power using the different interface 
circuits, for operation at resonance and off-resonance,! normalized to the maximum 
extractable power Piim (see Sect. 2.4.4). The plots represent a sectional view of the 
3D plots shown in the previous sections at a fixed coupling factor. At resonance, the 
extracted power using the SECE interface is almost as high as the absolute maximum 
power Plim. In the following, only the output power at the optimum load resistance is 
considered. Whereas the resistor load harvests maximum 80 % compared to Piim, the 
rectifier-capacitor extracts maximum 70 % compared to Piim. So, the SECE interface 
is capable of increasing the harvested power to about 124 % of the rectifier-capacitor 
load. 

When the harvester is driven at off-resonance, on one hand, the extracted power 
reduces greatly compared to the excitation at resonance. But on the other hand, the 
power harvested using the SECE interface is improved to 400 % of the maximum 
power harvested using the rectifier, which is exactly the ratio of the corresponding 
equations (3.26) and (3.17): 


! For the operation at off-resonance, the mass displacement amplitude has been chosen such that the 
optimum output power using the resistor load corresponds to the value Pre; max at far off-resonance 
as given in Table 6.4. 
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Fig. 3.10 Normalized harvested power with resistor load, rectifier-capacitor load and SECE load, 


for operation at a resonance and b off-resonance (r? = 0.56 x 10-9 (As/mP^, z = 45um at 
off-resonance, harvester values taken from Table 6.1) 


PSECE _ 

OR 

rect,max 

Similarly to the discussion in Sect. 3.4.1, the dependency of the gain achievable 

by SECE compared to the rectifier load has its origin in the backward coupling force 
amplitude I” Vp. An excitation at resonance results in a high piezoelectric voltage 
amplitude Vp, hence the backward coupling force is high, leading to a small beam 
deflection. Since SECE damps the beam deflection stronger than the rectifier, the 
harvested power using SECE can only be increased slightly, at least for the given 
coupling factor. For devices with higher I? values, SECE would even reduce the 
harvested power, as depicted in Fig. 3.9. For excitation at off-resonance, Vp is much 
lower compared to the excitation at resonance, hence the backward coupling force is 
also smaller. In that regime, the power improvement achievable by SECE 1s higher 
(400 % in theory). In Chap. 6, this phenomenon will be verified by measurements. 


3.4.3 Conclusion 


As a result, it is impossible to make a general statement which interface circuit per- 
forms best. Actually, according to (3.32), this decision depends on many parameters 
like the piezoelectric capacitance, the stiffness, the damping, the load resistor and 
the excitation frequency. In addition, it is important if the piezoelectric harvester is 
driven at resonance or at off-resonance. It can be concluded that using the SECE load 
only makes sense for piezoelectric harvesters with low electromechanical coupling. 
For harvesters with high coupling, the simple diode rectifier represents the better 
choice. 
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Chapter 4 
Theory of the Proposed PSCE Circuit 


In this chapter, the theoretical background of the proposed Pulsed Synchronous 
Charge Extractor (PSCE) circuit is discussed. The principle way of power extraction 
used in the PSCE is based on the Synchronous Electric Charge Extraction (SECE) 
technique, which has the potential to increase the power harvested from piezoelec- 
tric devices significantly compared to the well-known diode rectifier. However, one 
important point has been (deliberately) left out in the corresponding theoretical dis- 
cussion: the influence of parasitic resistances on the harvested power. Since the main 
property of the SECE technique is to transfer the energy in a very short time slot 
compared to the excitation period, the current flowing during this short time can 
easily get as high as several ten milliamperes for transferred power in the order 
of few milliwatts. Subsequently, parasitic series resistances in the current path can 
produce significant power losses which reduce the harvested power. The techniques 
presented in the following have the goal to diminish these undesired power losses. 
Besides numerical analysis, a linear approximation is discussed, allowing to make a 
statement about the superiority of these techniques for any circuit parameter set. 


4.1 Basic Operation Principle 


The main ideas and methods shown in this chapter have been developed during 
a student research project [1]. Figure 4.1 shows the functional schematics of the 
PSCE interface, composed of a NVC, three switches, and an inductor. The NVC 
"rectifies" the piezoelectric voltage, although it is no real rectifier in the sense that 
it is capable of blocking reverse current. By means of proper controlling of the 
switches, the inductor temporarily stores the energy which has been extracted from 
the piezoelectric harvester during the transfer process, before is 1s released to the 
storage capacitor and dissipated by the load. See Chap. 5 for a detailed description 
of the circuit blocks of the PSCE interface. 
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PSCE Interface 





y, (D 


Transfer Process 


z-— — Circuit 











— Veo — |Vp| 0 to to +T /2 t 


Fig. 4.2 Basic SECE operation principle. a State diagram. b Transient characteristics 


A similar interface circuit was initially presented in [8], whereas the switch S3 was 
not present. The functionality of this switch was adopted by the NVC, in the sense 
that current was able to flow through the body diodes of the MOSFET transistors 
inside the NVC, of course producing undesired power losses and the danger of latch- 
ups. For that reason, S3 is introduced in the proposed design in order to bypass the 
NVC when necessary. 

The basic operation principle can be explained by means of Fig. 4.2. The switching 
techniques are completely determined by the piezoelectric harvester voltage Vp and 
the inductor current /7,, which are called state variables. The transfer process starts at 
t = to when the piezoelectric voltage has reached a maximum value Vp = Vpo = Vp. 
The energy stored in the piezoelectric capacitance and the inductor at this moment 
are 


1 

Epo := Ep (to) = 5 VPoCp (4.1) 
1 

Ero := EL (to) = 5Llto — 0Ws. (4.2) 


During the transfer process, the energy Epo is distributed between the energy 
stored on the piezoelectric capacitance, Ep, the energy stored on the the inductor, 
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EL, and the energy transferred into the buffer capacitor, Epuf: 
Epo = Ep + EL + Eput (4.3) 


with 
t 


Euf := / Voutlfput dí for fo <t < tg, (4.4) 


to 


where tg denotes the end of the transfer process. When the energy Epo has been 
transferred completely into the buffer capacitor at the end of the transfer process, the 
following conditions hold due to (4.3): 


Eput (te) = Epo 
Vp (tg) = OV (4.5) 
Ir (tg) = 0A. 


Hence, the trajectory for the transfer process must start at 
Zo := (Vp, IL) (to) = (Vpo. OA) 


and stop at 
Zg :— (Vp, IL) (te) = (OV, OA). 


One possible shape which fulfils this condition is shown in Fig. 4.2. The trajectory 
shape is arbitrary and depends on the control of the switches, but in general, at the 
beginning of the transfer process, the inductor current must increase somehow so that 
the inductor can store energy, and after a certain instant, it has to decrease somehow 
in order to release the stored energy to the buffer capacitor. 


4.2 Energy Loss Approximation Approach 


Assuming an arbitrary interface circuit with Pin and Pout being the average input and 
output power, its efficiency is calculated according to 





n= == = = n, (4.6) 


where Pw denotes the average power loss inside the interface circuit. Generally, the 
average power can be written as 
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tE 
— 1 
P P (t) dt 
tE — 10 J Z 
to 
E (t 
— Er) vith E(f) = 0Ws. (4.7) 
fg — fo 


Using (4.7), Eq. (4.6) can be rewritten as 


.. Ew (tg) 
Ein (fg) 


dec. (4.8) 


Since the functional schematic shown in Fig.4.1 is modelled through ideal com- 
ponents, the energy loss is zero, and thus the efficiency equals 1. However, the energy 
loss in real interface circuits is always non-zero due to the parasitic series resistances 
inherently present, resulting in a efficiency smaller than 1. The energy loss in the 
NVC, the inductor, the switches and the buffer capacitor can be modelled by their 
parasitic series resistances RNvc, Rac, Rs1, Rs2, Rs3 and Rpur, as shown in Fig.4.3. 

Since it can be difficult to obtain an analytical expression for Ew, an approximation 
method is presented in the following. Let Rs be the sum of all parasitic resistances 
in an arbitrary path through the interface circuit and let J; be the current flowing 
through Rs. The energy loss caused by Rs during the transfer process is calculated 
through the integral over its power loss: 


tg (Rs) 
Ew(Rs) = J P (t, Rs) dt 

Os 
PSCE Interface — — 
| NVC he L Rig S; Rs | 
mA A E a ! 
: R 93 ai Rout 
NVC Y " i 

Sg ve Jm RS Vus 
: Cour 
: Rss Rs, ! 
© a 


Fig. 4.3 Schematics of the PSCE interface including parasitic series resistances 
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te (Rs) 
—Rs. J I? (t, Rs) dt (4.9) 


to 


It is important to mention that both the duration of the transfer process and the 
current through the resistor are functions of the resistance Rs. For that reason, a new 
variable p(Rs) substituting the integral is introduced: 


[E (Rs) 
p(Rg) — J I? (t, Rs) dt. (4.10) 
Os 


Hence, the energy loss in (4.9) and its first derivative can be written as: 


Ew (Rs) = Rs - p (Rs) (4.11) 


La (Rs) Sp) H-Rs- P (Rs) (4.12) 
IRs S) =P S S UR. S). . 


Since the parasitic resistances are assumed to be small, i.e. Rs ~ (0-10 2), a first 
order Taylor approximation around Rs = 0 Q is performed using (4.11) and (4.12): 


QOEw 
Ew jin (Rs) := Ew (0 2) + Rs - ORs (0 $2) 


= Rs-p(0Q) with Ew(0Q) —0Ws. (4.13) 


By means of the approximation Ew (Rs) ~ Ew. iin (Rs), it is possible to calculate 
the energy loss of a small parasitic resistance Rs by simply calculating the current 
through the resistance for the assumption that Rs = 0 Q, i.e. itis sufficient to calculate 
p (0 $2) according to (4.10). That means that the resistance is assumed to have no 
influence on the characteristics of /7, which greatly simplifies the calculation of the 
current for any switching technique in contrast to the lossy case where Rs > OQ. 
This method allows to compare the efficieny of different switching techniques for 
various variables, under the condition that the approximation error is reasonably 
low. Moreover, it gives insight into the loss mechanisms, as will be discussed in 
section 6.4.3. For the sake of simplicity, the abbreviation p :— p (0 $2) is used in the 
following. 


4.3 Switch Configurations 


By means of three appropriate switch configurations SC1, SC2, and SC3, the interface 
circuit shown in Fig. 4.1 performs the energy transfer from the piezoelectric harvester 
to the inductor and to the buffer capacitor. In this section, the equivalent circuits of 
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the switch configurations are analyzed, and for each switch configuration indicated 
by n € (1, 2, 3}, the respective state equation 


o Vp 
Ot 

fn (Vp, IL) := for ne {1, 2, 3} 
olr, 


at 


and its solution 


F, (t, Zon) :— (i) (t, Zon) | for ne {1,2,3} 
L 


n 


is calculated. The arbitrary starting point of the solution is indicated by 
ZOn = (Vpon, Iron) for ne ae 2; 3h: 


The trivial switch configuration where all three switches are turned off is not con- 
sidered since this does not imply an energy transfer. For the following analysis, a 
commonly used abbreviation is 


1 
JLCp 


In order to simplify the analysis, the interface circuit depicted in Fig. 4.1 can be 
simplified twofold: First, the NVC can be omitted under the assumption that Vp > OV 
during the transfer process. Second, the piezoelectric harvester can be substituted by 
its capacitor Cp because the current flowing into Cp due to external vibration is 
assumed to be much smaller than the current flowing out of Cp during the short 
transfer process. 





Wo := 


4.3.1 SCI 


Figure 4.4a shows the equivalent circuit of SC1 when switch S, is turned on. This 
represents a Cp — L tank circuit. From the constitutive equations of the storage 
elements 


o Vp 
.— — —7] 
P 3f L 
ol 
‘Pe UL — Vp, 
Ot 


the following state equation can be derived: 
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(b) 





Fig. 4.4 Illustration of switch configuration SC1 (S, on, S2 and S3 off). a Equivalent circuit. 
b Trajectories in the state space 


OVE TL 
Ot C 
=|" pev s (4.142) 
a. | (v 
at E 
with the initial condition 
"P Oey a). (4.14b) 
Ir Iroi 


The solution of the initial value problem (4.14) is taken from (A.13) as derived 
in Sect. A.1: 


L " 
Vpo1 cos(oot) — Iro14/ C5 sin(wot) 


; (4.15) 
Itoi cos(wot) + Vroi4/ >> sin(cot) 


V 
Fi €, zo) = (9) mw = 


For a fixed starting point Zo), the two transient functions t > Vp, (t, Zo1) and 
t +> Ir (t, Z91) represent harmonic oscillations with the amplitude 


L 
0 Cp 


Cp 


Vpi (201) := ,/ Vài + 2 and Tri (zo) :— Ioi Vio > (4.16) 


and the average values Vp; := OV and Iz; := 0A. Due to a phase difference between 
both oscillations, in the state space, they represent concentric elliposidal trajectories 
with the center point (Vpi. I 11) — (0V,0A). The semi-axes of the ellipses are 
defined by (4.16). Figure 4.4b shows the trajectories for two different starting points. 
For each starting point (Vpo1, /z01) = (Vpo, OA), exactly one trajectory exists which 
can be used to start the transfer process. 
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OV Vout 2Vouf Vp 


Fig. 4.5 Illustration of switch configuration SC2 (S2 on, S, and S5 off). a Equivalent circuit. b 
Trajectories in the state space 


4.3.2 SC2 


Figure 4.5a shows the equivalent circuit of SC2, where S» is turned on. The following 
corresponding state equation can be derived: 


d 4 
ot C 
= Y =: f» (Vp, Ir) (4.172) 
ofr, Vp — Vout 
ot L 
with the initial condition 
YP) (Os) = zoz :— ( P2}. (4.17b) 
Ir Iro2 


The solution of the initial value problem (4.17) is taken from (A.14) as derived 
in Sect. A.1: 


Vp? Vout + (VP02 — Vout) cos(wot) — Ir02,/ c sin (wot) 
F? (t, Z02) := |, J (65202 = = l 
s (Vpo2 — Vout) 4/ 7 sin(wot) + Iro» cos(cot) 
(4.18) 


For a fixed starting point Zo2, the two transient functions t œ> Vp (t, zo?) and 
t +> Ir» (t, Zo») represent harmonic oscillations with the amplitude 


Z^ L 

Vp2 (202) :— |f (Vpo2 — Vout) + loc. and 

T, 2 CP 2 

I12 (Zo2) :— 4/ (Vpo2 — Vout) PA + Trop (4.19) 
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and the average values Vp := Vou and Iz2 := OA. Due to a phase difference 
between both oscillations, in the state space, they represent concentric elliposidal 
trajectories with the center point (V po, I 12) = (Vpbur, OA). The semi-axes of the 
ellipses are defined by (4.19). Figure 4.5b shows the trajectories for three different 
starting points. In contrast to SC1, SC2 can be used both to start and to finish the 
transfer process. With the starting point zo? = (Vpo, OA) for all Vpo > 2Vbuf, SC2 
can be used to start the transfer process. For all Vpo < 2Vbur, the corresponding 
trajectories miss the origin, so the prerequisite as defined in Sect. 4.1 is violated. In 
order to finish the transfer process, the corresponding trajectory must go through the 
origin, so only the trajectory with the semi-axis Vp? (Zo2) = Vpuf can be taken for 
that purpose. In Fig. 4.5b, this condition is met by the middle trajectory. 


4.3.3 SC3 


Figure 4.6a shows the equivalent circuit of the SC3, where S5 and S3 are turned on. 
The particularity about SC3 is the fact that it is only valid for Vp = OV. Thus, the 
following state equation can be derived: 





“r OV /s 
t 
ET fs (Vp, Ir) (4.20a) 
T : 
with the initial condition 
"P (Os) = Zo3 := d (4.20b) 
Ir 1103 
The solution of the initial value problem (4.20) is 
Vp3 OV 
F; (t, Zo3) := t, Z03) := 4.21 
3 (£, Z03) e 03) (noo Ya (4.21) 
Fig. 4.6 Illustration of switch (b) 
configuration SC3 (S2 and i; 
$5 on, Sı off). a Equivalent 
circuit. b Trajectories in the 
state space (a) 
dL 
0A 
Cyr |n 
0V Vp 
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As shown in Fig. 4.6b, the trajectories representing (4.21) are straight lines heading 
downwards on the /z axis. For all /ro03 > OA, the trajectory goes exactly through the 
origin, and thus, SC3 can be used to finish the transfer process. 


4.4 Switching Techniques 


A switching technique means the sequence of two switch configurations (or two 
phases) such that the trajectory in the state space starts in Zo := (Vpo, OA) and ends 
in Zg :— (OV, OA). The first switch configuration is denoted by SCa € {1, 2} and the 
second switch configuration is denoted by SCb c 12, 3}, so that the transition point 
between both switch configurations is indicated by tap and the end of the transfer 
process is indicated by tgap. The transfer process starts at fo. 


4.4.1 Energy Loss Approximation 


Each switching technique is characterized by its approximated energy loss. For this 
purpose, a resistor R, models the accumulated resistive losses in the current path of 
each swich configuration n € (1, 2, 3}: 


Ry := Rac + RNvc + Rsi 
Rz := Rac + RNvc + Rs2 + Rout (4.22) 
R3 := Rae + Rs3 + Reo + Rouf. 


For each switch configuration, the current through R, equals the current through 
the inductor. Since the characteristics of the current during SCb depends on the 
current during the preceding SCa, a piecewise current /7;5 is defined describing 
the characteristics for the whole switching technique. The total energy loss of the 
switching technique SECEab, Ewap, is the sum of the energy loss portions of SCa 
and SCb. From the definition (4.13), the approximated energy loss is calculated as 


tab tEab 


E Wab,lin ‘= Ra | Tha (t) dt + Rp J i (t) dt 
to fab 
= RaPa,ab + RbPb,ab, (4.23) 


whereas the integrals over the squared inductor current have been abbreviated as 
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lab tEab 
Pa,ab = | is (f) dt and Pb,ab = J i (t) dt. 
fo fab 


Substituting the energy loss Ew in (4.8) by the approximated energy loss Ewap iin 
from (4.23), the approximated efficiency can be written as 


a,ab R Pb,ab 


a b 
Epo Epo 








Nab lin = 1 — R (4.24) 


Hence, the goal of the following analysis is to find closed-form expressions of 
Pa, ab and pp ab for each switching technique. In the following analysis, commonly 
used abbreviations are 

i V 1 
wo i= and x:= = and Epo := = p Vg. (4.25) 


/LCp Vout 2 








4.4.2 SECEIS 


Xu et al. [8] were the first to present a CMOS-integrated interface circuit based on 
SECE13. Simulation results of an improved CMOS circuit have been presented by 
the author of this book at the PowerMEMS 2008 conference [3], and measurement 
results have been shown at the PowerMEMS 2010 conference [5]. In SECE13, the 
transfer process starts with SC1, so that the state variable pair (Vp, IL) propagates on 
the trajectory from zo :— (Vpo, 0A) in the left direction around the origin, according 
to Fig. 4.4b. When the capacitor voltage has reached Vp = OV, SC3 is enabled, 
causing the capacitor voltage to decrease, as depicted in Fig. 4.6b. When the origin 
is reached, which is equivalent to the de-energized state, all switches are turned off 
and thus the transfer process ends. 

The state equation of SECEI3 is piecewise defined by the state equations of SCIT 
and SC3, as derived in (4.14) and (4.20): 


o Vp 
Ot fı (Vp, Ip) : V, 0V 
_ Jf Ve, Ir) P (4.262) 
olr f5 (Vp, IL) : Ir > 0A 
Ot 


with the initial condition 


V V. 
( ^ ) (to) = zo :— ( m and Vpo > OV. (4.26b) 
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Fig. 4.7 Capacitor voltage and inductor current characteristics of SECE13. a State diagram. 
b Transient diagram 


Figure 4.7 shows the numerical solution of the initial value problem (4.26). In 
Fig.4.7a, the piezoelectric voltage is plotted against the inductor current, and in 
Fig.4.7b, both state variables are plotted against the time. SC1 is active from tọ until 
t13, and from £13 until tg13, SC3 is enabled. The following subsection discusses the 
calculation of these time intervals. 


4.4.2.1 Analytical Solution 


Using the solutions of SC1 and SC3 as derived in (4.15) and (4.21), the analytical 
solution of the initial value problem (4.26) is as follows: 


Mid pus B (t — to, Zo) : to Et E tia (4.27) 


Ir13 F; (t — 113,213) : t13 < t < te13, 


Z3: Mur (t13) 
Iria 


denotes the state at the transition point t13. Using (4.15) and (4.21), the functions 
F; (t — to, zo) and F; (t — t13, Z13) can be evaluated as follows: 


where 


" Vpo cos (co (t — to)) 
F, (f — to, zo) = t — fo, Zo) = pC 
1 ( 0. Zo) [e ) ( 0s 2 Vpo T sin (wo (t — to)) 


I35(f — 113,213) =E; (f — fg13, ZE) 


Us 0V 
= ( ) (t — tg13, ZE) = Vout (4.28) 
173 CL (t — fg13) 
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whereas the equivalence relation from (A.18), 


Nea (C= 1935713) — B3 01513: Zn). 


has been applied to F3. This equivalence relation avoids the calculation of the state 
variables at the starting point of SC3, z13, which is more difficult than evaluating the 
state variables at the end of SC3, zg. The times t13 and tz13 can now be calculated 
using the continuity condition at the transition instant t = 13: 
! 
Fy (413 — to, Zo) = F3 (tis — tE13, ZE). 
The continuity condition of the capacitor voltage: 


! 
Vpi (f13 — to, Zo) = Vp3 (tia — tE13, ZE) 
Vpo cos (wo (t13 — to)) = OV 


leads to the transition time = 
f13 = to + / LCp - 2 (4.29) 


In the same way, using the result (4.29), the continuity condition of the inductor 
current 


! 
Iri (t13 — to, Zo) = Ira (t13 — tE13, ZE) 


Cp. Vout 
Vpoy/ T; sin (wo (f13 — to)) = e (fg13 — £13) 


leads to the end of the transfer process 


[E13 = t13 + V LCp - x sin (o (t13 — to)) 


Tm 
= f13 + y LCp - xsin 5 


^ 
zl 


= nue E 7 dest m (4.30) 


4.4.2.2 Approximation of the Energy Loss 


The next step is the calculation of the approximated energy loss 


Ew13,in = Rip1,13 + Kapa,13 
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according to (4.24). From the inductor current in SC1, 


Cp . 
Ipi3 (t) = fri (t — to, Zo) = Vero, T sin (wọ (t —19)) for tọ €t X fia, 
(4.31) 


the integral over the squared inductor current is calculated as: 


t13 
C 
P13 = 12,4 (t) dt = Veo E aj (1 — cos (2a (t — to))} dt 


10 


= — — sin (zo) E 


Cp JU 
—E = 4.32 
Poy 5 (4.32a) 


where t13 has been substituted by the result found in (4.29). 
From the inductor current in SC3, 


buf 
Iri = Ir} (t — tE13, ZE) = —L (t — tg13) for ti3 <t < tf13 


the integral over the squared inductor current is calculated as: 


tE13 


Dade | ts (t) dt = Vout i (t — tg13)^ dt 


t13 [13 


2 
“but 





731 (tg13 — t13)^ 
Ba 5 VP 5 [ue 
Cp 2 
EE UR ecd 4.32b 
Poy y. 3* ( ) 


where t13 and tg135 have been substituted by the results found in (4.29) and (4.30). 
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4.4.3 SECE23 


The theory behind SECE23 has been presented at the Eurosensors 2009 conference 
[4], and first simulation results of a CMOS-integrated interface circuit using this 
technique have followed at the PowerMEMS 2009 conference [2]. The fully autono- 
mus CMOS chip implementing SECE23 which is discussed in Sect. 5.8 has been 
presented at the ESSCIRC 2011 conference [6], and published in the Journal of 
Solid-State Circuits [7]. 

When SC2 is used to start the transfer process, the state variable pair (Vp, IL) 
propagates on the trajectory from zog :— (Vpo, OA) in the left direction around the 
point (Vpu¢, OA), as shown in Fig. 4.5b. The capacitor voltage can reach the origin only 
for the case that for the semi-axis of this trajectory, the condition Vo (Zo) > Vbuf 1S 
met, which is equivalent to Vpo > 2Vyyr from (4.19). After the trajectory has reached 
the /7, axis under that condition, SC3 is enabled. According to Fig. 4.6b, the capacitor 
voltage thus decreases, until the origin is reached and the transfer process ends. 

The state equation of SECE23 is piecewise defined by the state equations of the 
switch configurations 2 and 3, as derived in (4.17) and (4.20): 


QVp 
ot u b (Vp, Ip) : Vp » OV (4.332) 
al, |  |fs(Vp,Ip) : Ip» 0A ' 
3t 
with the initial condition 
V V 
( E ) (fg) = zo :— ( E and Vpo > 2Vyt. (4.33b) 


Figure 4.8 shows the numerical solution of the initial value problem (4.33). In 
Fig.4.8a, the piezoelectric voltage is plotted against the inductor current, and in 








/ 





OV = Vout 2VoutVpo |Vpl to 3  lE23 t 


Fig. 4.8 Trajectory of SECE23. a State diagram. b Transient diagram 
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Fig. 4.8b, both state variables are plotted against the time. SC2 is active from fo until 
t23, and from f23 until tg23, SC3 is enabled. The following subsection discusses the 
calculation of these time intervals. It becomes directly evident that the maximum 
inductor current is smaller compared to SECE13, which causes less losses in the 
parasitic series resistances, as discussed in Sect. 6.2.4.4. 


4.4.3.1 Analytical Solution 


Using the solutions of switch configurations 2 and 3 as derived in (4.18) and (4.21), 
the analytical solution of the initial value problem (4.33) is as follows: 


F> (t — to, : lo <t<t 

P Oe 2 (t — to, Zo) Q0 €t (4.34) 
F3 (1 — 123,203) : f93 E 1 < te, 

where 

Ira 


Z53 :— (m (t23) 


denotes the state at the transition point t23. Using (4.18) and (4.21), the functions 
F- (t — to, zo) and F3 (t — t23, Z23) can be evaluated as: 


Vout + (Veo — Vout) COS (co (t — to)) 
F2 (t — to; Zo) = (e) (f — to, Zo) = Cp . 
L2 (Veo — Vout) rf Sin (wo (t — to)) 
F3 (f — 123, Z23) = F3 (t — tg23, ZE) 
0V 
Vpa 
— [dj Zp) = V i 
(15) (f — tgo3, ZE) ( buf ee " 


L 


whereas the equivalence relation from (A.18), 

F; (t — t23, 223) = F3 (t — tẸ23, Ze), 
has been applied to F3. This equivalence relation avoids the calculation of the state 
variables at the starting point of SC3, z23, which is more difficult than evaluating the 


state variables at the end of SC3, zg. The times t23 and tg23 can now be calculated 
using the continuity condition at the transition instant t = t23: 


F2 (t23 — to, Zo) = F3 (t3 — tE13, ZE) - 


The continuity condition of the capacitor voltage: 
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! 
Vpo (t23 — to, Zo) = Vp3 (t23 — tk23. ZE) 
Vout + (Veo — Vout) COS (wo (t — t9)) = OV 


leads to the transition time 


t23 = to + v LCp - arccos 





—] 
(4.35) 
x— | 


In the same way, using the result (4.35) and the relation from (A.20), 


sin (arccos u) = v 1 — u?, 


the continuity condition of the inductor current 


! 
Ij2 (t23 — to, Zo) = Ira (t23 — tE23, ZE) 
Cp. Vout 
(Veo — Vout) T; sin (wo (t23 — t9)) = —L (t23 — tg23) 
leads to the end of the transfer process 


fE23 = t23 + V LCP - (x — 1) sin (wo (t23 — to)) 


1 
nic Huc al | aaa TO) 


—|] | 2 
= to + LC? - arccos 1 + /LCp-x,/1—-. (4.36) 
x— x 


4.4.3.2 Approximation of the Energy Loss 





The next step is the calculation of the approximated energy loss 
Ew23,lin = R2p2,23 + R3p3,23 (4.37) 
according to (4.24). From the inductor current in SC2, 
1123 (t) = Ir2 (t — to, Zo) = (Veo — Vout) [Sin (wo (t — t9)) for to < f < 193, 


the integral over the squared inductor current is calculated as: 
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[53 123 
2 2 Cp 1 
D2,23 = | fp53 (t) dt = (Vpo — Vowt) T {1 — cos (2o (t — to))} dt 
to 10 
Éro (4 _ Vout i t3 — t m (209 (t23 — to)) 
= — 4 1 — fy — —— sin (20 — 
5 Vo 23 — to Jong 0 (t23 — to 








Epp | Cp 1\? —] 1. =| 
= — ,/— J 1 — — arccos — — sin | 2 arccos 
L E X x— 1 2 x— 1 
Cp jx 1 1 2 
= Epo, | — ]—-] arccos +-,/l——-}. (4.38a) 
L x l1—x x x 


where t23 has been substituted by the result found in (4.35), and the trigonometric 
relation from (A.21), 





I 
5 Sin (2 arccos u) = uv 1 — u? for —1 <%= 1, 


has been used. 
From the inductor current in SC3, 


buf 
Iro3 = [73 (t — tgo3, ZE) = D (f—fgo3) for $3 €t < tm23, 


the integral over the squared inductor current is calculated as: 








1E23 v2 tE23 
pios = | Gos (t) dt = 1» J (f — tg23)^ dt 
123 123 
— l Var 3 
-737]2 (tg23 — 123) 





=F J2 2a 2y! (4.38b) 
= EPO T. 3* " , . 


where t23 and tg23 have been substituted by the results found in (4.35) and (4.36). 


4.4.4 SECEIZ 


When the transfer process starts at zo :— (Vpo, OA) using SC1, SC2 can be enabled 
at an appropriate instant t12, leading to SECEI2. The transition instant £j» has to 
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be chosen such that in SC2, the origin can be reached in order to finish the transfer 
process properly. Only the trajectory with the semi-axis Vp» (Z2) = Vout; depicted as 
the middle trajectory in Fig. 4.5b, can be chosen to finish the transfer process. Hence, 
in order to reach this trajectory using the SC1 trajectory, the condition Vpo < 2Vbuf 
has to be met for the starting point zo :— (Vpo, 0A). 

The state equation of SECE12 is piecewise defined by the state equations of the 
switch configurations | and 2, as derived in (4.14) and (4.17): 


o Vp 
Ot fi (Vp, Ip) : V Vp (t 
_ Jf (^p. In) P > Vp(ti2) (4.392) 
olr f> (Vp, IL) ; Ir > 0A 
ot 
with the initial condition 
"Wiwen Od and 0 < Vpo < 2Vpuf. (4.39b) 
Ir, OA 


The capacitor voltage at the transition point, Vp (t12), can be derived using the 
relation from (4.5), 


Pout (fg12) = Epo, (4.40) 


where tg12 indicates the end of the transfer process. (4.40) means that the energy 
transferred into the buffer capacitor is equal to the energy initially stored on the piezo- 
electric capacitance. Based on (4.4), the energy transferred into the buffer capacitor 
can also be written as 


tE12 TE12 
Eput (te12) = | Voutlput dt = / {—Voburlp} dt 
t12 12 
= —V»ut : Cp (Vp (te12) — Vp (t12)] , (4.41) 


where —/p = wf is the current flowing out of the piezoelectric harvester and into 
the buffer capacitor. Thus, combining (4.40) and (4.41), the capacitor voltage at the 
transition point can be calculated as 


2 


V _ V Po 4.42 
P2) = 577 (4.42) 
u 





Figure 4.9 shows the numerical solution of the initial value problem (4.39). In 
Fig.4.9a, the piezoelectric voltage is plotted against the inductor current, and in 
Fig. 4.9b, both state variables are plotted against the time. SC1 is active from fo until 
t12, and from fj»? until tg12, SC2 is enabled. The following subsection discusses the 
calculation of these time intervals. 
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Fig. 4.9 Trajectory of SECE12. a State diagram. b Transient diagram 


4.4.4.1 Analytical Solution 


Using the solutions of switch configurations 1 and 2 as derived in (4.15) and (4.18), 
the analytical solution of the initial value problem (4.39) is as follows: 


F; (f — to, Pd pee f 
(ia) a ee DS ME (4.43) 
Tri F; (t — 12,212) : tra € f < try, 


where 
V 
Zjo i (52) (112) 
L12 


denotes the state at the transition point t12. Using (4.15) and (4.18), the functions 
F; (t — to, zo) and F>? (t — t12, Z12) can be evaluated as: 


Vpo cos (wo (t — t9)) 


" 
F; (f — to, zo) = (i) (t — to, Zo) = 


Cp. 
Vpo T; sin (wo (t — to)) 


Fo (t — t2, Z12) = Fo (t — te12, ZE) 


V 
= Dn (t — tE12, ZE) = 
L2 


Vout — Vout COS (Wo (t — tg12)) 
Cp . 
Vpo Vout EN (wo (t — tg12)) 
whereas the equivalence relation from (A.18), 
Fy (t — t12, Z12) = Fo (t — tgi2, Ze), 


has been applied to F2. This equivalence relation avoids the calculation of the state 
variables at the starting point of SC1, z15, which is more difficult than evaluating the 
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state variables at the end of SC2, zg. The times fj» and ftg4» can now be calculated 
using the continuity condition at the transition instant t = t12: 


| 
Fy (t12 — to, Zo) = Fo (2 — tg12, ZE) . (4.44) 


Further, the transition condition (4.42) and the continuity condition (4.44) can be 
combined, only considering the capacitor voltage: 


2 





V 1 Vio 4.45 
pi) Vout (4.45a) 
u 
! 
= Vp, (f12 — to, Zo) (4.45b) 
! 
= Vp (ti2 — tg12, ZE) . (4.45c) 


From (4.45b), the transition instant £1» can be derived: 





y2 
P0. — Vpo cos (œo (t12 — to)) 
2 Vbuf 
ti2 = to + y LCp - arccos : (4.46) 


From (4.45c), the end of the transfer process, tg12, is calculated: 





2 
Vbo — u = 
= Vout — Vbur COS (wo (t12 — fg12)) 
2 Vbut 
x x? 
tg12 = to + y LCp - arccos 5 + y LCP - arccos (: — =) ; (4.47) 


where t12 has been replaced by the result found in (4.46). 


4.4.4.2 Approximation of the Energy Loss 
The next step is the calculation of the approximated energy loss 


Ewi24in = Ripi,12 + Rop2,12 


according to (4.24). From the inductor current in SC1, 


Cp . 
[112 (t) = fpi (t — to, zo) = Vpoy/ T; sin (wo (t —t9)) for fo X t X ti, 


the integral over the squared inductor current is calculated as: 
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t12 t12 
Cp 1 
P1122 = ie (f) dt = Veo 7- = fu — cos (2m (t — t9))] dt 
to 10 
Epo L|. 
= — 1112 — to — —— sin (20 (t2 — to)) 
L 200 
Epo Cp X 1 . X 
= — ,/ — {arccos — — — sin (2 arccos z) 
L L 2 2 2 
Epo |C 2 
E arccos m ]— = l (4.48a) 
g L 2 2 4 


where tj2 has been substituted by the result found in (4.46), and the trigonometric 
relation from (A.21), 


1 
5 Sin 2 arccos u) = uv 1 — u? for Lege 


has been used. 
From the inductor current in SC2, 


Cp. 
Ip12 (t) = Ir (t — te12, ZE) = Vouty/ PR (wo (f —tg12)) for tig <t <tr, 


the integral over the squared inductor current is calculated as: 


tE12 tE12 
P2,12 = J I7, (t) dt = Var ] J {1 — cos (2o (t — tg12))} dt 
f12 f12 
Epo 1 1, 
me i re — fo m sin (2w9 (t1? — 2) 


Cp 1 x 1. x 
= Epo,/ — — 1arccos (1 — — ] — = sin f 2 arccos į 1 — — 
L x? 2 2 2 
Emp E- 1% HN: EE E (4.48b) 
= — 4—arccos{l1— —]—-—|{-- = ——1. ; 
PVT | x2 2 x 2 4 


where 11» and tgı2 have been substituted by the results found in (4.46) and (4.47), 
and the trigonometric relation from (A.21), 


1 
5 Sin (2 arccos u) = uv 1 — u? fo —1l<u<1, 


has been used. 
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4.4.5 Combined Switching Techniques SECE1223 and SECE1323 


In the previous sections, three different switching techniques were analyzed. Because 
of the initial condition (4.26b), SECEI3 can be applied for any voltage ratio x :— 
Vpo/Vbur Within the interval 


x13 :— (x e R|x > 0]. 


In contrast to SECEI3, due to the initial conditions (4.39b) and (4.33b), the usage 
of the novel switching technique SECEI2 is limited to 


x2 :={xER|O<x< 2}, 
and the usage of SECE23 is limited to 
x23 := {x E R| x > 2}. 


Because of the property 
X12 U x23 = x13, 


SECE12 and SECE23 can be combined to the new switching technique SECE1223. 
One can also combine SECEI3 and SECE23 to the new switching technique 
SECE1323, whereas SECEI3 is activated for x € x1? and SECE23 is enabled for 
X € X23. These two new combined switching techniques are evaluated in the follow- 
ing section. 


4.5 Evaluation 


In this section, the switching techniges SECE1223 and SECE1323 are evaluated with 
regard to their approximation error and their efficiency improvement compared to 
the commonly used SECE13 technique. All the plots have been made using the para- 
meters shown in Table4.1, which have been determined roughly from simulations 
of the PSCE chip. 


4.5.1 Generic Efficiency Approximation 


As derived in Sect. 4.2, the approximated efficiency of the switching technique 
SECEab € 113, 23, 12] can be written as 


Pa,ab Pb,ab 
a — Rp 
Epo Epo 








Nablin = | — R ; (4.49) 
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Table 4.1 Circuit parameters Parameter Value 

used for efficiency evaluation Rac 43Q 
Rnvc 15.1 Q 
Rs} 4.6 €2 
Rs2 7.8 €2 
Rs3 4.5 €2 
Rout 0 2 
Cp 19.5 nF 
L 10 mH 
Vout 25V 


The closed-form expressions of p; ab and pp ap given in (4.32), (4.38) and (4.48) 


are all functions of the coefficients Epo/ fp and x. Hence, 4.49 is a function of 


e the parasitic resistances R4 and Rp, 


‘ox — LP 
e the voltage ratio x = us 


e and the coefficient ,/ e. 


The energy loss within the parasitic series resistances is the main reason for an 
efficiency gap < 1: the larger the parasitic resistance R, and Rp, the higher the energy 
loss and thus the lower the efficiency. Their influence on the approximated efficiency 
depends on the voltage ratio x. Since the current is flowing through the inductor in 
each switch configuration, the inductor resistance Rac has the largest influence on 
the approximated efficiency. So, Rac should be made as small as possible in order to 
achieve a highly efficient transfer process. 


Due to the coefficient FE , the approximated efficieny increases for decreasing 
Cp and increasing L. For a given harvester, the piezoelectric capacitance Cp is fixed, 
so the only possibility to optimize the approximated energy is to increase the inductor 
inductance L. However, since the inductance of a physical inductor can be only 
increased by using more wire windings, the inductor resistance increases at the same 
time, unless a larger wire diameter is chosen. Thus, the only way to increase the 
inductance without affecting the parasitic series resistance is to increase the inductor 
size as well, probably resulting in unpractical circuit implementations. 

The approximation performed in (4.49) is only suitable for small parasitic 
resistances, otherwise the approximation error gets too large. However, under the 
assumption that the approximation errors of the three switching techniques are equal, 
the efficiency improvement of SECE1223 and SECE1323 in relation to SECEI3 can 
be reasonably calculated even for a large approximation error. In the following, this 
issue will be investigated in detail. 
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4.5.2 SECE 1223 


In this subsection, the SECE1223 technique is evaluated regarding its approxima- 
tion error and efficiency improvement with regard to the commonly used SECEI3 
technique. The goal is to make a statement about the quality of the approximation. 


4.5.2.1 Approximation Error 


In order to evaluate the quality of the efficiency approximation, the approximation 
error of the SECEI223 technique is investigated in the following. Therefore, the 
exact efficiency is determined by means of numerical solution methods which are 
explained in section A.4. 

Figure 4.10a shows the approximated and numerically solved efficiencies of 
SECE1223 and SECE13 as a function of x. For x < 2, the two curves of SECE1223 
are almost identical, indicating that the approximation is almost perfect in that regime. 
Forx > 2, the two curves of SECEI223 spread more and more for increasing x which 
means that the quality of the approximation drops. The two curves of SECEI3 start 
spreading already beginning from x = 0. At a first glance, looking at x > 2, the 
numerically solved efficiency curves of SECEI223 and SECEI3 seem to have an 
almost parallel characteristic; the same is true for the approximated curves. For 
almost all values of x, the approximated efficiency is smaller than its numerically 
solved counterpart, which is indicated by the approximation error being negative, as 
shown in Fig. 4.10b. However, even for x = 14 (which is the maximum ratio at which 
measurements have been conducted), the approximation error of both SECEI3 and 
SECE1223 is still below 8 % which is quite acceptable. 

For the curves in Fig. 4.10, the voltage ratio x — vae had been used as the function 
variable, whereas all the other parameters were assumed to be constant as given in 
Table 4.1. Since the inductance L of the external inductor has a great influence on 
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Fig. 4.10 Comparison of SECE1223 and SECE13 techniques as a function of x. a Numerically 
solved and approximated efficiency. b Approximation error 
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Fig. 4.11 Comparison of SECE1223 and SECE13 techniques as a function of L (x = 4). a Numer- 
ically solved and approximated efficiency. b Approximation error 


the characteristics of the SECE techniques, it is also interesting to know how the 
efficiency and the approximation error depend on L, whereas the voltage ratio x is 
fixed. 

Figure 4.11 illustrates the corresponding curves. From Fig.4.1la, it can be 
observed that for a very large inductance, the efficiency tends to 1, and the efficiency 
improvement of SECE1223 compared to SECEI3 tends to zero. This correlation 
becomes directly obvious from the fact that all the p coefficients in (4.23) are pro- 
portional to 1/./L, as calculated in Sect. 4.4, which means that the average power 
loss is also proportional to 1 /4/L. A demonstrative explanation for this coherence can 
also be given: The duration of the transfer process increases as L increases, because 
the time constant of the resonance circuit, J/ LCp, increases. This means that the same 
amount of energy is transferred in a longer time, leading to a smaller current J; and 
thus finally to a smaller average power loss in the parasitic series resistance. This is of 
course only true if the power consumption of the control circuitry of the PSCE chip 
is not taken into consideration. Realistically, there exists an optimum for a certain 
value of L, because the average power consumption of the control circuitry which is 
turned on during the transfer process depends on the duration of the transfer process. 
Also, a realistic view should incorporate the dependency of the parasitic wiring resis- 
tance from the inductance L, because Rc is a function of L. Since this correlation 
depends on many parameters like the length and the diameter of the inductor wire 
which varies for different types of inductors, it is impossible to provide a general 
expression. So, here Rac is assumed to be constant with regard to L. 

The curves representing the approximated efficiency are always below the curves 
representing the numerically solved efficiency. As shown in Fig. 4.11b, the approx- 
imation error dramatically increases when L decreases. Whereas the error for both 
switching techniques is below 15 96 for L = 1 mH, the error of SECEI3 is 70 96 for 
L = 200 pH. For an inductance of about L = 160 uH, according to the approxi- 
mation, SECEI3 reaches zero efficiency, whereas the numerical solution gives an 
efficiency of about 45 %. The reason for this behavior is the fact that the duration of 
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the transfer process decreases for decreasing L, resulting in a higher current /7, and 
hence a higher average power loss. In that regime, the assumption that the duration 
of the transfer process is approximately independent of the power loss is no longer 
valid, thus the approximation is getting worse. It can be concluded that a sensible 
approximation can be made for inductors with L > 1 mH. 


4.5.2.2 Efficiency Improvement 


The approximated efficiency improvement of SECE1223 compared to SECEI3 can 
be calculated as 


411223,lin(X) = N1223,lin C). — 13,1in CÓ 
PEIG) n P2,12(X) 











: = 
Lg PLBO) |p, P3813) _ PE DER EXD 
= : i 
Epo Epo Ry P223 ) 4 Ry P23! ) "— 
Epo Epo 


_ 1 [Ri (P1,13 — P112) Œ — Rop2,12G) + R3p3,130%) : x € x12 
Epo |Ripi,130) — Rop2,23(x) + R3 (p3,13 — p3,23) (x): x € x93 . 
(4.50) 


After substituting the resistances R1, Ro and R3 as defined in (4.22): 


R1 := Rac + RNvc + Rs1 
Ro := Rac + RNvc + Rs2 + Rout 
R3 := Rac + Rs3 + Rs2 + Rouf, 


Eq. (4.50) can be written as 


Rac(pi,13 + P3,13 — P1,12 — P2,12) Œ) 
+ Ruvc(p1,13 — P1,12 — p2,12) 2 
+ Rsı (p1,13 — pi, 12) 0) 
+ Rsapa,13(x) 
A | |. ol + (Rs2 + Rout) (p3,13 — p2,12) 0) [hex 
111223,lin(X) = —— 
Epo | Rooil (P1,13 + P3,13 — P2,23 — P3,23) (x) 
+ RNvc(p1,3 — p223) Œ) 
+ Rsipi,13(x) 
+ Rs3(p3,13 — p3,23) (x) 
+ (Rs2 + Rout) (93,13 — P223 — P3,23) (x) : x € x23. 
(4.51) 


Finally, the approximated efficiency improvement can be expressed as 
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1 

A1223,lin (X) = Epo (RacPac + RNvcPNvc + Rsipsi + Rs3ps3 + (RsoRowt) Ps2,but)» (4.52) 
P 


whereas the following coefficient functions are used: 


(p1,13 + 3,13 — P1,12 — 2,12) (x) : x € xi 


Pác(x) := 
(P1.13 TPD — 9293 — p3,23) (x) : else 
PNvc(X) — (p1,13 = P1.12 — p2,12) (x) : o4 Eus 

(71,13 — p223) (x) : else 

psi (x) —_ (p1,13 — p1,2) (x) : dq bu 
p1i,13 (x) : else 

mee p3,13(X) |: X EX 
(p3,13 — p3,23) (x) : else 

_ X Ae Rs 
piura = (p3,13 — 2,12) (x) Do ER 


(p3,13 —P2.23 — pa23) (x) : else 


In order to evaluate the coefficient functions regarding their sign, they are simply 
plotted against x, as illustrated in Fig. 4.12. Since all of the coefficient functions have 


got the factor Epo/ CP in common, they are normalized to this factor. According to 


the figure, all of the coefficient functions are greater than zero for any value of x. 
In addition, the figure suggests that for x — oo, they asymptotically tend to values 
greater than or equal to zero. 

Figure 4.13a shows the numerically solved and the approximated efficiency 
improvement between SECE1223 and SECE!I3 as a function of x. In the range 
0 < x < 2, the two curves are almost identical, indicating a very small approxima- 
tion error, as observed in Sect. 4.5.2.1. The largest efficiency improvement of 5.5 % 


Pcoil 
PNVC 
PS1 
Ps3 


PS2 buf 





Fig. 4.12 Graphical representation of the normalized coefficient functions 
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Fig. 4.13 Approximated and numerically solved efficiency improvement of SECE1223 compared 
to SECE13. a Function of x. b Function of L (x = 4) 


can be achieved at x = 2. For x > 2, the numerically solved efficiency improvement 
is slightly smaller than the approximated one, but the basic shape is the same. 

A similar obervation can be made for the efficiency improvement as a function of 
L, as shown in Fig. 4.13b. The two curves show a similar shape, however the difference 
gets worse for a small inductance. E.g., For L = 100 pH, the approximation gives 
an efficiency improvement which is almost twice as large as the numerically solved 
efficiency improvement. For L = 1 mH, the error in the efficiency improvement is 
around 33 %, and increasing the inductance further brings the error into a negligible 
range. These findings correspond to the approximation error as depicted in Fig. 4.1 1b. 
From Fig. 4.13, it seems that both the numerically solved and the approximated 
efficiency curves tend to zero for x — oo and L > oo. 


4.5.2.3 Conclusions 


From the findings made in Sects. 4.5.2.1 and 4.5.2.2, it can at least hypothesized 
that A71223 Jin (X) > 0 and Amj223 (x) > O for all x > 0, which means that the 
approximated and the numerically solved efficiency of the new SECEI223 is larger 
than the approximated efficiency of the formerly known SECE13, independently of 
the voltage ratio x :— Vpo/Vypyr and the parameters 


Cp, L, Rac, Rsi, Rs», Rs3, Rout, Veo and Vout. 


Compared to the time-consuming numerical approach, a great benefit of the effi- 
ciency approximation is the possibility to evaluate the efficiency quickly for different 
variables. Figure 4.14 illustrates the efficiency improvement between SECE1 223 and 
SECE13 as surface plots as a function of the voltage ratio x, the inductor resistance 
Rac and the inductance L. In both plots, there is a “ridge” at x = 2, indicating 
the points of highest efficiency improvement. From Fig.4.14b, it can be seen that 
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(a) (b) 





Fig. 4.14 Surface plot of the approximated efficiency improvement of SECE1223 compared to 
SECE13. a Function of x and Rac. b Function of x and L 


the ridge is most distinct at low inductances, and there the curve drops rapidly for 
x < 2. At high inductances, the overall level of achievable efficiency improvement is 
low, and there are no distinct areas. Focusing on Fig. 4.14a, it becomes evident that 
SECE1223 achieves the highest benefit compared to SECE13 at high parasitic series 
resistances for any voltage ratio. Looking at (4.52), the efficiency improvement is 
directly proportional to the parasitic series resistances. In contrast to Fig. 4.14b, the 
ridge at x = 2 has no point of outstanding distinction. 


4.5.3 SECE1323 


For x > 2, the SECE1323 and the SECE1223 techniques are identical, which can be 
verified by comparing Figs.4.15 and 4.10a. Hence, the corresponding conclusions 
for SECE1223 drawn in Sect. 4.5.2.3 are also valid for SECE1323. Due to realization 
aspects discussed in Sect. 4.6, the SECE1323 technique is actually applied for the 
PSCE chip (see Sect. 5.8). Therefore, in the following, the SECE1323 technique will 
be discussed while focusing on practicability. 


4.5.3.1 Approximation Error 


As it is the preferred load variable for the measurements shown in Sect. 6.4, the 
SECE1323 characteristics will be predominantly discussed in terms of the buffer 
voltage Vbuf instead of x in the following. The piezoelectric voltage amplitude is 
fixed to Vpo = 12 V. Figure 4.16 shows the numerically solved and the approximated 
energy loss of SECE1323 as a function of the buffer voltage for two different inductor 
series resistances. The first important obeservation is the discontinuity at Vbuf = 
Uu which indicates the border between SECE13 and SECE23, according to (4.57): 


Whereas SECE23 is active for Vpyp < "o, SECE13 is enabled for Viuf > "o. 
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Fig. 4.15 Numerically solved and approximated efficiency of SECE1323 and SECEI3 techniques 


As discussed before, the SECE23 technique causes less losses under the presence 
of parasitic series resistance than the SECEI3 technique. This is why the energy 
loss changes significantly at the border between SECE23 and SECE13. When the 
parasitic series resistance is higher, as can be seen in Fig.4.16b, the overall energy 
loss is higher compared to Fig. 4.16a, but there is the potential to save more energy 
due to the usage of SECE23 compared to SECE13, as already seen in Fig. 4.14b. In 
Fig.4.16, this observation finds its graphical correspondence in the different jump of 
the power loss at Vyyr = T. 

Generally, it can be seen that the energy loss increases for decreasing buffer 
voltage, and for Vbuf — O, the energy loss tends towards infinity. This is obvious 
when we look at the corresponding equation for the approximated energy loss of 
SECE23, according to (4.37): 


(a) | 


Ew (uU Ws) 





Vout (V) Vout (V) 


Fig.4.16 Energy loss of SECE1323 as a function of the buffer voltage (Vpp = 12 V). a Rac = 4.3 Q. 
b Rac = 23.7 Q 
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Fig. 4.17 Efficiency of SECE1323 as a function of the buffer voltage (Vpo = 12 V). a Rae = 4.3 Q. 
b Rae = 23,7 G2 


£w23lin = Rop2,23 + K3p3,23, 


Cp ix* 1 1 2 
p223 = Epo,/ — 31|41— - ] arccos TN e 
L x ]—-x x X 
E [Cp 2 1 2 5 
= —— >œ —X —— i 
P3,23 PO L 3 z 


according to (4.38a) and (4.38b), and R and R3 are the total parasitic series resis- 
tances for the first and the second phase of the transfer process. In order to observe 
what happens when Vpyr tends to zero, we have to determine the limits of p? »3 and 


p3,23 for x — oo: 
Cp 
lim p2,23 = Epo, —— - 
X—-> CO 


lim p3,23 = oo. 
X—00Q 


where 





N| A 


As a result, the limit of Ew23 Jin is infinity for Vbuf — O. At the same time, the 
approximation error is increasing. 

Figure 4.17 shows the efficiency improvement as a function of the buffer voltage. 
Again, the discontinuity in the energy loss at Vous = m indicates that a certain 
amount of energy can be saved using SECE1323 compared to SECE13. For a larger 
parasitic series resistance depicted in Fig.4.17b, the approximation error increases 
compared to Fig. 4.17a. The harvested energy tends to zero when the buffer voltage 
tends to zero. This means that in practical applications, a certain minimal buffer 
voltage must be present, otherwise the ohmic losses may exceed the output power 
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of the harvester, leading to the fact that no energy can be harvested or even that the 
storage capacitor is drawn empty. 

4.5.3.2 Efficiency Improvement 

For the SECE1323 technique, the efficiency improvement can be written as 


411323,lin Œ) ‘= 71323,lin GO) — 113,lin (X) 


0 : 4 € X]2 
p1,13x) P3,13(x) 

= Rj —_— RR3—01—— — 1. 222309 p3,23 x) 

Epo Epo R ——_ + R3 = xX EX3 
Epo Epo 
= I 0 s: weal 
Epo |Ripi,136) — R2p2,23(x) + Ra (p3,13 — P323) C. : x € X3. 
(4.54) 


After performing all the calculation steps shown above, the approximated effi- 
ciency improvement can be finally expressed as 


1 
A1323,lin (X) = Epo (RacPac + Rnvcpnve + Rsipsi + Rs3ps3 + (Rs2Rvuf) Ps2,but), (4.55) 


whereas the following coefficient functions are used: 


( ) >: X XD 
Pac) :— 
i (P1.13 q P3,13 — P223 — 3,23) (x) : else 


0 . X € X12 
PNvc(x) := 
(p1,13 — p2,23)(x) : else 
0 Xx E X12 
psi(x) := 
Pi,13(x) : else 
0 à X E X]? 
ps3 (x) := 
(p3,13 — p3,23) (x) : else 
0 : X S AD 
PS2,buf (X) :— 
(p3,13 — p2,23 — pa23) (x) : else 


(4.56) 


The corresponding efficiency plots are shown in Fig.4.18. As can be observed, 
the combined SECE1323 offers an improvement only for x € x23 (see Sect. 4.4.5), 
which is equivalent to one of the following expressions: 
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Fig. 4.18 Approximated and numerically solved efficiency improvement of SECE1323 compared 
to SECEI3 
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Fig. 4.19 Efficiency improvement of SECE1323 compared to SECE13 as a function of the buffer 
voltage (Vpo = 12 V). a Rac = 4.3 2. b Rac = 23.7 Q 
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In that regime, the results shown in the previous section about SECE1223 can 
be directly taken to explain SECE1323. For x € x12, SECE1323 has no benefit 
compared to SECE13. 

Figure 4.19 shows the efficiency improvement of SECE1323 compared to SECE13. 
Whereas no efficiency improvement can be achieved for Vbuf > “po because SECE13 
is active in this range anyway, the largest efficiency improvement can be observed 
at Vine = m. At this point, the approximation error is at minimum. For decreasing 
buffer voltages which is similar to increasing values of x — ve, the approximation 
error increases, corresponding to Fig. 4.18. Also, it becomes evident that the approx- 
imation error increases for increasing series resistance, as discussed in Sect. 4.5.1. 
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The figures suggest that peak improvements of 5 % for Rac = 4.3 Q or even more 
than 10 % for Race = 23.7 Q are possible. So, the best operating point for the buffer 
voltage is around Vpyp = “pO whereas it has to be ensured that this point is not 
exceeded. For Vbuf — O, the numerically solved efficiency improvement tends to 
zero, whereas the approxmated efficiency improvement is approaching a fixed value. 


4.5.3.3 Conclusions 


It can be concluded that under the prerequisite that x is not too large, the approxi- 
mation shows very good agreement with the numerically solved curves, even for the 
large parasitic series resistance. As x increases towards infinity, the approximation 
error approaches infinity as well. In practical applications, a buffer voltage of no 
smaller than 1.5 V and piezoelectric harvester voltages up to 20 V can be expected, 
resulting in x < 14, for which the approximation error is acceptable, as shown in 
Fig.4.18. 


4.6 Realization Aspects 


Due to the promising results presented in Sect. 4.5, it is worth considering a practical 
realization of the new switching techniques SECE1223 and SECE1323. In order to 
implement SECE1323, it is sufficient to evaluate the condition 


Veo > 2Vbut (4.58) 


in order to exactly know whether to use SECE13 or SECE23, according to (4.33b). 
Equation (4.58) can be realized using a voltage comparator. The switching voltage 


Vp (t13) = Vp (t23) =O V (4.59) 


which determines where the change between the two switch configurations should 
happen is the same for SECEI3 and SECE23 and can be easily detected by means 
of a zero crossing detector. 

SECE1223 also requires that (4.58) is evaluated in order to select the appropriate 
SECEI2 or SECE23. In case SECEI2 is activated, the switching condition 


2 


" 
Vpi2 :— Vp (t12) = — (4.60) 
2 Vbut 





defines where to change between switch configurations 1 and 2, as initially shown in 
(4.42). Compared to (4.59), (4.60) is much more complicated to implement. First, due 
to the parasitic series resistances present in practical realizations, energy will prob- 
ably remain in the piezoelectric capacitance after the transfer process has finished. 


92 4 Theory of the Proposed PSCE Circuit 


Second, the initial voltage Vpo :— Vp (to) has to be stored on an auxiliary capacitor 
for later evaluation, which implies that the circuit evaluating (4.60) has to consume 
extremely low power in order to preserve the voltage on that auxiliary capacitor. 
Third, a multiplier which is commonly known for its high power consumption is 
necessary for the calculation of Veo: 

Another possibility to realize the SECEI2 is an active control of the stored transi- 
tion voltage Vp12, which permanently checks whether the voltage on the piezoelectric 
capacitance, Vp, and the inductor current /7, get null at the same time, indicating a 
proper operation. If this is not the case, the transition voltage Vp;2 can be varied by 
a small amount. The transition voltage Vp12? can be stored on a capacitor or digitally 
in a register, which would require an additional analog-to-digital converter. 

To conclude, the implementation of the SECE1223 is very problematic. One goal 
of the developed PSCE chip is the energy automous operation without supplying extra 
power which exceeds the power delivered by the piezoelectric harvester. Another goal 
is the ability to accept piezoelectric harvesters with a low output power of 10 uW — 
30 uW, which limits the PSCE chip power consumption to several microwatts. Since 
these goals contradict the usage of power-consuming signal processing or control 
electronics, the author has focused on the realization of the SECE1323. 
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Chapter 5 
Implementation of the PSCE Circuit 
on Transistor Level 


Itis always good practice to verify the findings of a theoretical discussion by means of 
practical measurement results. Hence, the Pulsed Resonant Charge Extractor (PSCE) 
circuit derived in the previous chapter has been implemented on transistor level in 
a 0.35 um technology. In the past years, this technology has well established in the 
field of energy harvesting, and the fabricated chips have shown a good performance. 
For the circuit design process, the electrical equivalent circuits of the piezoelectric 
harvester as derived before haven been used. One of the main goals was to reduce the 
average power consumption within the circuits into aregime where it can be neglected 
compared to the power converted by the mechanical energy harvester. One general 
strategy to achieve this goal is to aggressively scale down power consumption in the 
blocks which have to be active all the time, accepting slow operation. In constrast, 
higher power consumption can be allowed for blocks which have to react faster, but 
can be deactivated most of the time. 

In the following, the operation principle of the circuit blocks 1s explained individ- 
ually in detail. Wherever possible, the blocks are grouped by their functionality, and 
the explanations refer to their functionality at top level. After all the blocks have been 
discussed, the top level implementation shows how all these blocks are combined to 
a top-level design. At the end of the chapter, photographs of the fabricated chip are 
shown. 


5.1 Power Switches 


Since the SECE1323 technique relies on redirecting the current flow through different 
paths, power switches play a major role in the PSCE circuit. In the following, all 
important aspects regarding the used power switches are highlighted. 
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5.1.1 Realization as MOSFETs 


Generally, the switches used within the PSCE circuit are realized as enhance- 
ment metal-oxide semiconductor field effect transistors (MOSFETs) driven in linear 
region. Depending on the doping of the MOSFET channel, two types of MOSFETs 
exist: The n-channel MOSFET placed within a p-well or p-substrate (substrates are 
commonly p-doped), and the p-channel MOSFET placed within an n-well. In order 
to turn MOSFETs on, 1.e. to establish a conducting channel between the drain and 
the source terminal, the following condition must be fulfilled: 


Ves] > IVrl. (5.1) 


where Vgs denotes the voltage between the gate and the source terminal (also called 
gate-source voltage), and Vy indicates the threshold voltage of the MOSFET device. 
The gate-source voltage and the threshold voltage are both positive for n-channel 
MOSFETs and negative for p-channel MOSFETs. In contrast to the less commonly 
used depletion MOSFETs, enhancement MOSFETs can be used as normally-off 
switches, which means that they are open (turned off) with no applied voltage at its 
gate. In fact, a gate-source voltage 


|Vas| = 0 (5.2) 


has to be applied to turn the MOSFET off safely, otherwise significant undesired 
sub-threshold conduction may result. Figure 5.1 shows the different configurations 
of n-channel and p-channel MOSFETS used as switches. 

Since one terminal of the switches S; and S3 is connected to ground (GND = 0 V) 
potential as shown in Fig. 4.1, they are realized as low-side n-channel MOSFETs Ms; 


GND —] 


) 


(a) (b) Vpp Vpp 
IVasl = Vpp 
Voo — c GND-] 2&2 
Vos= a 
Vop Vdp 
| IVesl=0 V 
An ol OV 
= EE 
Vgsz0 vw 


Fig. 5.1 a N-channel MOSFETs used as low-side switches and b p-channel MOSFETs used as 
high-side switches 
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and Mss. Hence, they can be turned off with a gate voltage equal to GND and turned on 
with a gate voltage equal to the supply voltage Vpp > Vr, where Vpp = Vpbur holds 
when the PSCE chip is active. So, GND and Vpp are the logic low (0) and high (1) 
levels to turn these switches off and on. In contrast, the switch S5 is a high-side switch 
because one terminal is connected to Vbuf. Of course it would be possible to take an 
n-channel MOSFET for S» as well, but this would imply two important drawbacks: 
First, under the assumption that a standard device is taken without isolated p-well, 
a body effect would occur resulting in an increased threshold voltage. This penalty 
would increase the minimum required supply voltage and reduce the achievable 
efficiency due to higher ohmic losses. Second, an additional voltage booster would 
be necessary in order to generate a voltage of minimum (VDD -- Vr,s2) in order to 
turn on the switch. In order to have a well conducting channel, this voltage would 
required to be much larger, especially when the body effect is an issue. The voltage 
boosting would need significant amounts of power which is very valuable especially 
in energy harvesting systems. 

Thus, the usage of the p-channel MOSFET Msz is preferred to realize S2. Accord- 
ing to (5.1), a gate voltage equal to 0 V turns S2 on when Vot > Vr,s2 is given, 
assuming the source terminal being connected to Vpyr. In order to turn Ms» off, the 
gate voltage has to be set to Vpp corresponding to (5.2). So, the logic low and high 
levels turn Ms» on and off in a reverse way compared to Ms; and Ms3. 

In the linear region, the parasitic channel resistance can be approximated as [6] 


1 


== (5.3) 
Un Cox T- (Ves — Vr) 


Kon 


where Uun, Cox and w denote the mobility of negative charge carriers, the oxide 


capacitance and the aspect ratio of the MOSFET, respectively. In order to obtain 
a highly efficient power transfer from the harvester to the buffer, a small channel 
resistance is desired. According to (5.3), this can be achieved by a large aspect 
ratio or a high gate-source voltage, provided that the other parameters are process 
dependent and cannot be changed. So, the length of the power switch MOSFETs, 
L, is set to the minimum value provided by the CMOS process, and the width W 
is chosen orders of magnitude larger. Simulations have shown that above a certain 
point, there is no significant increase in the efficiency or the efficiency even reduces. 
The reason for this behavior is the existence of the parasitic gate-source capacitance 
which has to be charged and discharged during each cycle, resulting in inevitable 
energy loss. Due to the fact that the parasitic gate-source capacitance is proportional 
to the aspect ratio, this energy loss increases for increasing aspect ratio. 
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5.1.2 Bulk Regulation 


In general, the bulk connection of PMOS devices used as a switch is critical in power 
circuits where the drain and the source terminal voltages are not well defined with 
respect to each other. If the bulk is hard-wired to one of its two terminals, forward- 
biased bulk diodes may occur, leading to excessive power loss and/or latchup, result- 
ing in circuit malfunction or even destruction [6, ch. 17.8]. 

For the specific PSCE circuit, the bulk connection of the PMOS Ms» is very 
critical in terms of energy harvesting efficiency. Hard-wiring of the bulk terminal 
to Vbuf would result in a body diode across Mg» as shown in Fig. 5.2a, creating a 
conducting path between the piezoelectric harvester and the buffer capacitor when 
the piezoelectric voltage is higher than the buffer voltage. This would result in a 
failure of the SECE technique, because the piezoelectric harvester cannot be kept in 
open circuit condition as required. Figure 5.2b shows the situation if the bulk terminal 
is connected to Vz—. In that case, a body diode in the reverse direction would appear 
across Ms», discharging the buffer capacitor when there is no piezoelectric activity. 

So, the bulk of Ms» has to be dynamically connected to the terminal with the 
higher potential in order to avoid the occurrence of bulk diodes. This can be achieved 
by means of the bulk regulation circuit [2] shown in Fig. 5.3. The relation of Vj. 
and Vyyr to each other determines the position of the source terminal of Mp; and 
Mp2. For Vr... > Vius, the left terminal of Mp; connected to V;_— 1s the source of 
Mp}, and its gate-source voltage is given by Vos; = Vbuf — Vr. For Vgsi > Vr, 
Mp; turns fully on, leading to Vy = Vj... Since the gate-source voltage of Mp2 is 
given by Vgs2 = Vr. — Vx = Q in this case, Mp2 is turned off. Vice versa, for 
VL- < Vy, Mp; turns off and Mp? is conducting, connecting the bulk of Mg» to 
Vout, SO that Vy = Vbuf holds. To conclude, the bulk regulation circuit implements 
the maximum function max (Vbuf, Vr, .). Proper regulation is only ensured when 
the potential difference between Vz— and Vpyr exceeds the threshold voltage of the 
regulating MOSFETs M p; and Mp2. Otherwise, only a weak sub-threshold current 
is available to regulate the bulk, extending the regulation time and thus increasing 
the probability for an undesired power loss within the body diode. 





MT MeT 


Fig. 5.2 Position of the body diodes of the p-channel MOSFET switch Ms», with the bulk terminal 
connected to the a buffer voltage and the b rectified piezoelectric voltage terminal 
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Fig. 5.3 Schematic of Bulk Regulation 
the bulk regulation circuit 
connecting the bulk of Ms2— 
max (Vbur, VL—) in order to 
eliminate bulk diodes 





5.1.3 Level Shifter 


Let's first motivate why a level shifter is mandatory in the PSCE circuit. When the 
rectified pizoelectric voltage exceeds the buffer voltage during the SECE charge 
process, i.e. VL+ = VL- > Viuf, the terminal of Mg? connected to Vr... becomes the 
source. Since the logic high level driving the gate of Ms» cannot rise above Vbuf, Ms? 
starts conducting in sub-threshold region since |Vas| > 0. When the piezoelectric 
voltage increases further, i.e. Vr... — Vbuf > Vr, Msz turns fully on and connects the 
piezoelectric harvester directly to the buffer capacitor. Since the buffer capacitor is 
now charged passively, the piezoelectric harvester can not be kept in open circuit as 
required by the SECE technique, leading to failure of this technique. 

A workaround for this problem is to set the logic high level of the Mg» control 
signal to Vz— for Vi... > Vbuf and to Vbuf otherwise. This functionality can be 
achieved by means of a common level shifter, where the high supply voltage is 
connected to max (Vbuf, Vr, -) which can be easily generated using the bulk regulation 
circuit discussed in Sect. 5.1.2. Since neither the bulk regulation nor the level shifter 
circuit consume static power, this method provides a power efficient work around for 
keeping Ms» turned off safely for piezoelectric voltages below or above the buffer 
voltage. 

The level shifter shown in Fig. 5.4a shifts the high level of the digital input signal 
Di, from Vpp to Vpp,max and outputs the result as Dou. The circuit operates as 
follows: For Di, being at high level Vpp, n-channel MOSFET My; is turned on, 
pulling the Vy node to GND. Since n-channel MOSFET Mmwo is turned off due to the 
inverted gate voltage, p-channel MOSFET Mp4 pulls Vy to the high supply voltage 
Vpp.max- Since the Vy node is connected to the gate of p-channel MOSFET Mp3, 
this transistor is turned off so that it does not interfere with the Vy node pulled to 
GND. The output inverter finally converts the GND potential to Vpp max. If Din is at 
low (GND) level, Mw; is off and My2 is on, setting the Vy potential to GND which 
turns Mp3 on and Mp4 off. Thus, Vy = Vpp,max and Vy = GND which leads to Dout 
level being at GND. As a result, the input level Vpp is shifted to Vpp, max, whereas 
the input GND level is fed through to the output. 

Figure 5.4b shows the configuration which is commonly used in the presented 
harvester chip. By connecting the Vpp, max terminal of the level shifter to the out- 
put of the bulk regulation circuit, the high level of the Din signal is converted to 
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(a) Vpp, max 






18/5.5 200/3.5 








Fig. 5.4 a Schematics of the level shifter circuit, shifting the high level from Vpp to Vpp max, and 
b common configuration with the bulk regulation output connected to the Vpp max terminal of the 
level shifter 


max (Vbur, V_—). As a result, the output signal Dout can safely keep several tran- 
sistors, e.g. Ms», turned off for the rectified piezoelectric voltage Vz— being at any 
level. The dimensioning of the level shifter transistors is critical in terms of the slew 
rate of the control signal when SECE23 is enabled. For achieving low power loss, 
Ms? has to be turned on instantaneously after the transfer process has been initiated. 
This means that the level shifter output must provide a fast decay from V;,.. to GND 
shortly after the peak of Vz— has been detected. At the end of the transfer process, 
the level shifter must set its output rapidly to Vbuf in order to turn off Ms2, avoiding 
unnecessary reverse currents discharging the buffer capacitor. 

Since Vz— can be much higher than the buffer voltage Vor, the level shifter tran- 
sistors MN; and My? have to be dimensioned asymmetrically, as shown in Fig. 5.4a. 
By designing Ms? much wider than Msp, the Vy node is pulled to GND stronger 
than the Vx node to Vpp, max, resulting in a fast decay of the output. Of course, the 
output signal is not fed directly to the gate of Ms» since this would exceed the drive 
capability of the level shifter by far. Instead, a cascade of CMOS inverters is used to 
drive the large power switch devices (see drivers subsection in Sect. 5.1). 


5.1.4 MOSFET Drivers 


Large power switch MOSFETs have a high parasitic gate-source capacitance Cas 
which has to be charged and discharged each time they are turned on and off. Thus, 
in combination with the parasitic series resistance of the gate connection line, Reger, 
an Rser — Cos low pass filter is established, causing a poor slew rate of the edges of 
the control signals. Since flat edges create high cross current in the CMOS inverters 
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Fig. 5.5 Low-voltage switch driver for Msı, Ms3 and the NVC MOSFETs 


that are often used to generate control signals, drivers are used to increase the slew 
rate. 

Usually, chains of CMOS inverters with exponentially increasing aspect ratio are 
used. On one hand, there are drivers required for driving the low-voltage power 
switches Msi, Ms3 and the NVC MOSFETs. The corresponding control signals are 
Ds, Ds3, DNvca and Dwvcp. Since the high level of these transistors equals Vor, 
the low voltage driver shown in Fig. 5.5 is used. The inverters are standard inverters, 
composed of standard low-voltage transistors. 

On the other hand, the power switch Ms» requires a more complex high voltage 
driver as depicted in Fig. 5.6, because the high level of its control signal Ds» has to be 
shifted to Vpp, max, as described in Sect. 5.1.3. So, each inverter has to be composed 
of high voltage MOSFETs which are supplied by Vpp,max. The n-channel MOSFET 
(NMOS) Mw ties Vx to GND during startup (see Sect. 5.5.1) in order to safely turn 
off Ms». Since the voltages of the nodes between the first five inverters are undefined 
during startup, the supply of these inverters has to be decoupled by the p-channel 
MOSFET Mp, otherwise extensive cross-currents may occur. After the PSCE chip 
has started up successfully, Dig is pulled low, turning off Mw and thus releasing 
Ds» for normal operation. 


5.2 SECE Selector 


The task of the SECE selector is to determine whether to activate SECE13 or SECE23, 
therefore it is the main component for the implementation of the combined switching 
technique SECE1323 (see Sect. 4.6). 
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HV Switch Driver 





Fig. 5.6 High-voltage switch driver for the power switch Ms? 
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Fig. 5.7 Schematics of the SECE selector 


As written in Sect. 4.6, SECE13 is activated for 


y <2 


€— 0< Vpo < 2Vwf 


V 
c— 0 < — < Var. 


and SECE23 is enabled for 
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Fig. 5.8 Trajectory of SC2 
for R = 0 and for R > 0 





YSZ 
4> Veg > 2Viuf (5.4) 
Vpo 
EN > Vout; (5.5) 


where x :— Vpo/ Vour. Equation (5.4) says that the piezoelectric voltage amplitude 
Vpo must be at least twice the buffer voltage Vbuf in order to invoke the efficient 
SECE23. An equivalent statement is that Vbuf must not be larger than half of Vpo in 
order to activate SECE23, according to (5.5). For any other condition, SECEI3 is 
enabled. 

A problem in the practical realization can occur when Vpo is only slightly (some 
millivolts) larger than twice V5yr. Mathematically, SECE23 would be activated in 
this case. But due to the losses R in the parasitic circuit resistances, it might hap- 
pen that in the corresponding switch configuration 2, the energy cannot be removed 
completely from the piezoelectric capacitance. The corresponding situation is illus- 
trated in Fig. 5.8. This means that no zero crossing of Vz, can be detected, and 
as a result, the PSCE chip stops operating. For this reason, in order to avoid those 
situations, the lower limit of x for enabling SECE23 is increased by 10 % in the real 
implementation: 


X522 
<=>  Vpo > 2.2 Vbuf (5.6) 
Vpo 
—— > Vig: a) 
—À 25 = Vt (5.7) 


The SECE selector shown in Fig. 5.7 checks for this condition and outputs a low level 
digital signal Dsgcg if SECE23 can be activated, otherwise it is high. The voltage 
divider R;-R» generates a voltage V; ,./2.2 which is compared to Vpyp permanently 
by a low-power voltage comparator. Since its structure is well-known, consisting of 
a differential stage, a common-source single-stage amplifier and a digital inverter, 
the comparator is not shown in detail here. The bias voltage Vista for this comparator 
is generated by the supply independent bias generator BGI (Fig. 5.26). 
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In case Dsgce has gone low during one excitation half period, it is reset to high 
after the transfer process by means of the RCD output signal Dgcp. The signal Dsgcg 
is evaluated within the switch control logic immediately after the peak of Vz has 
been detected. For testing purpose, the SECEI3 technique can be forced externally 
by engaging Den sEce13; thus Dsgcg is permanently high. 


5.3 Oscillation Cancellation 


Figure 5.9 depicts the equivalent circuit of a real inductor. In parallel to the inductance 
L, there is the frequency dependent core resistance Rpe which is composed amongst 
others of the eddy current losses and hysteresis losses. The wire resistance is modelled 
as the series resistance Ring. The element Cpar summarizes the parasitic capacitances 
between the terminals and between the inductor windings. 

Especially Ring and Cpar have non-negligible influences on the behavior of the 
entire harvesting circuit. Whereas Ring reduces the total efficiency due to conduction 
loss, Cpar induces oscillation at the dangling Vz+ node after the transfer process has 
finished. This is because the current /7, 1s slightly non-zero at the moment when all 
switches are turned off, causing the remaining energy 0.5 L7 5 to oscillate within the 
tank circuit formed by the parasitic components of the real inductor. The voltage 
amplitude of the resulting oscillation at the V; ,. node can reach several hundred mil- 
livolts and thus trigger several undesired transfer processes after the actual transfer 
process has taken place. In case the ZCD is not able to react on these small peaks fast 
enough, switch Ms; remains turned on permanently, causing the PSCE chip to stop 
operation. For that reason, an oscillation cancellation transistor Moc is connected in 
parallel to the inductor, as shown in Fig. 5.9. By switching Moc on immediately after 
the transfer process, the residual energy stored on the inductor is removed immedi- 
ately after the transfer process has finished. As a result, the oscillation vanishes and 
thus no instability problems can be observed for any buffer voltages. 


Fig. 5.9 Equivalent circuit Dos 
of a real inductor including 
shorting transistor Moc 







Inductor 


5.4 Negative Voltage Converter 103 


(a) I NVC, act V; à (b) INC. pass VL. 





Fig. 5.10 a Active and a b passive path of the negative voltage converter (NVC) connected in 
parallel (all transistors have a dimension of 10000/1), and c symbol of the NVC 


5.4 Negative Voltage Converter 


Since the piezoelectric harvester produces an alternating output voltage, a rectifi- 
cation stage is necessary in order to charge the buffer capacitor, although recent 
publications have proposed rectifier-free harvesting structures as well [3]. For the 
PSCE circuit, the rectification stage is implemented as a parallel connection of an 
active and a passive path as shown in Fig. 5.10a and b. Unlike real rectifiers using 
diodes based on the p-n (silicon) junction or the metal-semiconductor (Schottky) 
junction, this stage is not able to block current flowing in the reverse direction, thus 
it will be referred to as negative voltage converter (NVC) from here on. The basic 
design has been taken from [7] and adapted to the special needs of high voltage 
piezoelectric voltage conversion. 


5.4.1 Circuit Design 


The passive path is composed of two cross-coupled inverters connected such that 
the involved transistors are turned on and off by the piezoelectric terminals Vp} and 
Vp... Hence, the passive path is able to process the negative piezoelectric voltage 
during the startup phase where the buffer capacitor is insufficiently charged. The 
passive path starts conducting when |Vp+ — Vp_| > Vr, assuming that Vr is the 
threshold voltage of the n-channel and the p-channel MOSFETs. If Vp+ is greater 
than Vp_ by about one threshold voltage, Myo and Mp, turn on whereas Mws and 
Mp» are off, creating a conducting path from the Vp, to the Vz, terminal and 
from the GND to the Vp_ terminal. Assuming that a load is connected externally 
across Vz, and GND, the current circuit is closed and current can flow from one 
piezoelectric terminal to another. If Vp_ is greater than Vp by about one threshold 
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Fig. 5.11 Simulated NVC output voltage and output current of the passive and the active NVC, 
indicating which part of the NVC is conducting for which NVC output voltage a Vpyg = 1.5 V 
b Vou =3 V 


voltage, Mys and Mp» turn on whereas Mye and Mp; are off, establishing a current 
circuit from Vp_ to Vp. across an external load. 

Once the buffer capacitor voltage Vor is large enough to drive the control circuitry, 
the active path starts operating. If Vp., is greater than Vp_, the logic signal DNvca 
is set to high level, turning on the switches MN» and My3. At the same time, the 
switches MN; and Mw; are turned off by pulling the control signal Dyycp to GND. 
Thus, a conducting path from Vp. to Vp_ is generated, assuming an external load 
across the NVC. From the transient output current of the passive and the active NVC 
during the transfer process, as shown in Fig. 5.11, it can be seen which path of the 
NVC is conducting for which piezoelectric voltage. It becomes evident that the active 
path only conducts for low piezoelectric voltages, i.e. for Vp+/— < Vout — Vr, MN3/4; 
whereas Vr, mn3/4 denotes the threshold voltage of the two upper MOSFETs My3 
and My4. Due to the bulk effect My3 and My4 are exposed to, Vr Mw3/4 1S supposed 
to be quite large, reducing the input voltage swing. For a very low buffer voltage as 
depicted in Fig. 5.11a, the passive NVC conducts most of the time during the transfer 
process, and the active NVC takes over when the piezoelectric voltage is almost zero. 
In contrast, for a high buffer voltage as illustrated in Fig. 5.11b, the passive NVC 
starts conducting, but soon the active NVC takes over and conducts more current 
than the passive NVC. 

Hence, both parts of the NVC complement one another: On one hand, the passive 
NVC is responsible for conducting at the beginning of the transfer process when 
the piezoelectric voltage is large and is of course essential for passive startup. On 
the other hand, the active NVC cares about small piezoelectric voltages such that the 
NVC output voltage can reach zero, which is important for the zero crossing detector 
(ZCD) to be able to finish the first phase of the transfer process properly. 
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Fig. 5.12 Control circuitry for the active path of the negative voltage converter (NVC) 


5.4.2 NVC Control 


Figure 5.12 depicts the circuit used to generate the control signals for the active path 
of the NVC. Each piezoelectric terminal is connected to a source follower shifting 
the common mode up by the threshold voltage of the high voltage transistors Mp3 
and Mps, which is about 1.1 V. By means of the common mode up-conversion, the 
comparator is quite sensitive and offers a fast reaction time when the difference input 
signal Vp = Vp, —Vp.. changes sign. In order to avoid extensive oscillation and thus 
instable operation when Vp and Vp- are very close to each other, a comparator with 
positive feedback creating hysteresis is used [1, ch. 8.4], as shown in Fig. 5.13. The 
hysteresis window is defined by the dimensions of transistors Mpg through Mpo. In 


Vpp 





Fig. 5.13 Hysteresis comparator used in the NVC control circuitry [1] 
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the current application, a hysteresis window of 4-40 mV has shown to be sufficient to 
guarantee stable and robust operation. From simulations, the slew rates of the output 
have been determined as 2.9 MV /s for the rising edge and 3.6 MV /s for the falling 
edge. 

The source follower of the NVC control circuit and the hysteresis comparator are 
both biased by a supply independent voltage Vs: generated by the biasing generator 
BGI, as discussed in Sect. 5.7. Thus, persistent conditions for varying supply voltages 
is ensured. Since the piezoelectric voltage signal is low-frequent (50 — 500 Hz), as 
mentioned in Chap. 1, the biasing current can be in the nanoampere range (around 
55 nA). In order to drive the very large switch transistors of the NVC, MOSFET 
drivers are used to generate the control signals DNvyca4 and Dnvcp. 


5.5 Startup 


In this section, the components used to provide cold-startup capability are explained. 


5.5.1 Startup Trigger 


The startup trigger shown in Fig. 5.14 decouples the buffer capacitor from the supply 
of the control circuitry during the startup phase, i.e. when the buffer voltage is too 
low to supply the control circuitry. The minimum supply voltage or trigger voltage 
is mainly defined by the threshold voltage of the power switch transistors, because 
the supply voltage determines the high level of the switch control signals. In order to 
ensure proper startup functionality, the trigger voltage is chosen to be around 300 mV 
higher than the threshold voltage of the switch transistors, leading to a trigger voltage 
of around 1.4 V. The circuit design is based on [4] which has been extended and 
adapted such that the desired trigger voltage is achieved. 

The operation principle of the startup trigger circuit can be described as follows 
by means of the schematics shown in Fig. 5.14 and transient signals depicted in 
Fig. 5.15. Let Vr, and Vr,» be the threshold voltages of the involved n-channel and 
p-channel transistors. Starting from Vpy¢ = O V, Vy follows Vpur since Mpz is off and 
Vpp is zero so that Mp3 is slightly conducting in the sub-threshold region. So, Myı 
and My2 are conducting in the sub-threshold region as well, pulling Vx to GND. As 
Vout reaches Vr n, Mn; and Mnz are set into triode region since Vy keeps pulled to 
GND and Vr,» > Vr,n holds. At the same time, Dig goes high, actively blocking 
Mp4. As Ver reaches Vr, p, Mp2 enters saturation region. At some point of Vpur, 
Mp» gets "stronger" than Mw; so that Mp» enters triode region whereas Mw; enters 
saturation region, clamping Vx-(Vout — Vr, 4) due to the diode-connected M p,. The 
voltage where this happens is also called trigger voltage. 

As a result, the chain of inverters Ij—Is5 pulls Dyig-GND, activating Mp4 and thus 
coupling Vbuf—Vpp. At the same time, Vy goes low, disabling My2. One one hand, 
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Fig. 5.14 Startup trigger decoupling Vpuf and Vpp when Vpuf S 1.4 V 
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the current through the Mp;-Mpo-Mwi;-Mnmwpo branch is thus blocked. On the other 
hand, a hysteresis is introduced which prevents the startup trigger from oscillating 
immediately after the startup event has occured. This hysteresis can be observed 
by the fact that Vpp is decoupled at the shutdown trigger voltage which is smaller 
than the startup trigger voltage. As can be seen in Fig. 5.15, the low level of Vy 
is not perfectly O V due to the diode-connected My4, but that is sufficient to block 
the current. As a result, no current is consumed after startup has occured and the 
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PSCE chip is operating. The proper dimensioning of Mp2 and My; is very critical 
for achieving the desired trigger voltage. The diode connected transistors Mp; and 
My4 are necessary to attain trigger voltages higher than the threshold voltage of the 
involved transistors, as desired in the actual design. Simulations have shown that Dis 
has a very weak capability to switch to GND at the trigger point. Thus, five inverters 
with a very asymmetric drive capability (high/low level) have been introduced to 
afford this switching. In simulations, proper startup could be verified for fast and 
slow slew rates of the buffer voltage Vpuf. 


5.5.2 Bypass 


As explained in Sects. 5.1.2 and 5.1.3, the power switch transistor Ms? has to be 
strictly turned off all the time, except during the power transfer from the piezoelectric 
harvester to the buffer capacitor. On one hand, this characteristics is crucial for proper 
operation of SECE and avoiding buffer capacitor discharge, but on the other hand, 
passive charging during the startup phase when the PSCE chip is not yet operating is 
also prevented. For this reason, an additional bypass is necessary which establishes a 
path from the harvester to the buffer capacitor. The structure of this bypass is shown 
in Fig. 5.16a. Two p-channel transistors Mp; and Mp» are connected in parallel 
to Ms». Individual bulk regulations prevent the occurrence of parasitic bulk diodes. 
Whereas the diode-connected transistor Mp» allows the buffer capacitor to be charged 
passively, the series connected transistor Mp, is controlled externally by the signal 
Dig generated from the startup trigger in order to enable or disable the current flow 
through the bypass. During startup, Dyrig is high, and the inverter generates a low 
signal activating Mp;. As Dyig switches low when the control circuitry of the PSCE 


(a) Bypass 





Fig. 5.16 a Schematics and b symbol of the bypass which charges the buffer capacitor passively 
when Vpuf - 1.4 V (for simplicity, the bulk regulation of Ms» is not shown.) 
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chip starts operating, a high signal max (Vr, —., Vbuf) generated by the inverter and the 
level shifter turns Mp; off, disabling the current flow through the bypass. Without the 
level shifter, Mp; is turned on when V;,_ exceeds Vpuf and thus a positive gate-source 
voltage |VGs| = VL- — V» establishes. 


5.6 Switch Timing 


The timing of the involved switches is a crucial issue in order to provide proper 
operation. In the following, the circuit blocks used to derive the timing from the 
piezoelectric voltage signal are explained. 


5.6.1 Pulse Generator 


Figure 5.17a shows the architecture of the pulse generator. This element is used by 
the peak detector and the zero crossing detector (ZCD), amongst others, in order to 
generate well-defined pulses that can be further processed by a control logic. The 
pulse is generated by means of an XNOR gate, whereas one input of the XNOR 
gate is connected directly to the digital signal, and the other input is connected to 
the same signal delayed by an R-C element. Thus, the length of the pulse is defined 
by tg, = RC = 500k&$2 - 400fF = 200 ns for the given values. By the use of an 
additional NAND gate, the pulse is triggered by a rising edge only. As can be seen 
in Fig. 5.17b, a pulse with fixed length is generated for any length of the input pulse. 
Figure 5.17c shows the symbol used for the pulse generator throughout this book. 


(a) 





Fig. 5.17 a Architecture, b transient characteristic and c symbol of the pulse generator with rising 
edge detection 


110 5 Implementation of the PSCE Circuit on Transistor Level 


10/0.35 
= 


Pulse gen. 
: Dpp 


ure 


V stat 1 





4/10 


Fig. 5.18 Schematics of the proposed peak detector circuit 


5.6.2 Peak Detector 


The peak detector circuit has the task of determining the instant when the rectified 
piezoelectric voltage V; 4 has reached a maximum value. The proper detection of this 
instant is of great importance since it initiates the power transfer from the harvester 
into the buffer capacitor. Formerly used peak detector implementations mostly rely 
on voltage comparators, obviously restricting the maximum tolerated input voltage 
to the supply voltage, which is the case e.g. in [7]. Figure 5.18 shows the peak 
detector circuit proposed in this book using a current comparator instead of a voltage 
comparator, whereas a capacitor is used to convert the input voltage into an input 
current, i.e. the input voltage is derived. By doing this, the functionality is decoupled 
from the absolute value of the input voltage, as described in the following by means 
of Fig. 5.19 which shows the simulation of important signals. 

A capacitor Cs with a capacitance of 50 pF is connected between the rectified 
piezoelectric voltage V; , and the input of the current mirror Mx?-Mw3. When a 
positive voltage Vr, is applied, Cs is being charged by the current 


Vi+ — Vx VERS 
io de creen ee 5.8 
S S 3 5 " (5.9) 


whereas Vx is equal to the gate-source voltage of My2. Since /5 is in the nanoampere 
range, Mwo is driven in sub-threshold regime, and thus Vy can be neglected compared 
to VL+. The current mirror MN?-Mw; draws Zs out of the node Vy. At the same time, 
a very small fixed current Ibias pp ^ 23nA generated by the supply independent 
biasing circuit is injected into Vy via the current mirror Mp?-Mpj. Thus, both current 
mirrors form a current comparator with the inputs Z/s and Jpias, pp and the output at the 
Vy node. Assuming that /s > /pias,pp. the node Vy is at GND potential, because Myg 
is "stronger" than Mp;. As Vr. approaches a maximum, /s is decreasing towards 
zero according to (5.8). When Zs is falling below pias pp shortly before the maximum 
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Fig. 5.19 Simulated transient 
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is reached, Vy is pulled to Vpp. A sharp rising edge is generated by the inverting 
amplifier and the CMOS inverter, which triggers the pulse generator to produce a 
200 ns pulse, as explained in Sect. 5.6.1. Immediately after the pulse has been exerted, 
the DEN.zcp signal is set high, pulling Vx to GND. Hence, during the first phase of 
the transfer process, the capacitor Cs is discharged since Vz is decreasing to zero, 
resetting the peak detector for a new cycle. The transistor Mp3 precharges the output 
node Vy once directly after startup, in order to avoid an undesired peak detection 
causing circuit malfunction. 

In the following, the dimensioning of Jpias pp and Cs is discussed. The peak 
detector has to be permanently active in order to be capable to react on a voltage 
peak anytime. Thus, the biasing current Ibias, pp must be chosen very small so that 
the power extraction can also work efficiently for a low piezoelectric output power 
and the buffer capacitor is not drained empty too fast. But on the other hand, the 
biasing current must be capable of pulling the Vy node to high in a sufficiently small 
time, otherwise the peak would be detected too late. A 23 nA biasing current has 
shown to be a good tradeoff. The proper dimensioning of Cs is also important. The 
charge current /s must be sufficiently large. Otherwise, according to (5.8), Zs falls 
below /pias.pp too early resulting in a premature peak detection. This is especially 
a problem for a small piezoelectric voltage amplitude, and also a small excitation 
frequency, since both will result in a small derivative dV; | /dt. At some point, it 
might even happen that /s does not exceed Jpias pp at all, resulting in circuit failure. 
So, Cs has to be chosen sufficiently large for proper peak detection. A value of 50 pF 
has been chosen to ensure that functionality, without requiring too much expensive 
chip area. Best accuracy in terms of peak detection at the maximum piezoelectric 
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voltage is given at a vibration frequency of 175 Hz, corresponding to the resonance 
frequency of the MideV22B harvester (see Sect. 6.2.3). 

For testing purposes, an additional variable external offset current Joffset,pp can 
be set in order to trim the *peak detection" event. 


5.6.3 Zero Crossing Detector 


After the peak of the rectified piezoelectric voltage Vz has been reached and hence 
the energy transfer process has started, Vz is decreasing towards zero volts. This 
Zero crossing instant has to be detected properly in order to accomplish the transition 
to the switch configuration 3, i.e. the second phase of the transfer process. 

The implementation of the zero crossing detector (ZCD) is shown in Fig. 5.20. 
Here, a source-input comparator topology is used to detect the zero crossing of 
Vr... The Vx node represents the actual input of the ZCD which is connected to 
V,_+ during normal operation. The transistor MN; clamps Vx-(Vop — Vr,MN1) if 
Vr. > Vpp — Vr.mn1 So that mainly normal 3.3 V transistors can be used in the 
comparator stage. When the ZCD is disabled, transistor Mw? pulls Vx to GND to 
avoid crosstalk on Vy through the drain-source capacitance of the disabled transis- 
tor My . In addition, the Vpp rail is decoupled in disable state by Mp4 to avoid 
unnecessary power consumption. 

The operation principle is as follows and can be explained by means of Fig. 5.21. 
The gate of My3 and My3 is biased with a fixed voltage Vs ^ 1.2 V, whereas 
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Fig. 5.20 Schematics of the ZCD 
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Fig. 5.21 Simulated transient signals of the ZCD 


the exact value is not important—the only requirement is that Vstat2 is higher than 
the threshold voltages of Myn and Mna. Vstat2 1s generated on-chip within the bias 
generator BGI, as shown in Fig. 5.26. Assuming Vogfset,zcp = 0 V and Vy = Vr, > 
0 V at the beginning of the transfer process, the gate-source voltage of My4 is fixed to 
Vstat1, and the gate-source voltage of Myg equals (Vstatt — VL+), i.e. Vas4 > Vasa. 
As a result, the Vy node is pulled to GND since more current is drained from Vy than 
injected via the Mp; -Mp» current mirror. After two times inversion, Dzcp = 0 V in 
that case. 

As VL+ is decreasing, Vasa is increasing, and at some instant, Vgs3 exceeds Vgs4, 
i.e. Vas4 < Vasa. Hence, the Vy node is pulled to Vpp, which triggers the pulse 
generator to produce a short pulse on Dzcp. Immediately after this event, Dgw is set 
to GND, disabling the ZCD. 

Simulations have shown a more robust operation when this point is set to approxi- 
mately 0.4 V instead of 0 V, so the “zero crossing” is detected slightly too early. This 
is because a higher V; ,. peak voltage (Vp) means a faster decrease of V; +, resulting 
in a delayed reaction time. When the zero crossing detection voltage is set too low, 
it might happen that the zero crossing event thus cannot be detected properly for 
large Vp voltages. This is For this reason, My4 is implemented as a 5 V high voltage 
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Fig. 5.22 Function of the reverse current detector (RCD) as an active diode 


transistor exhibiting a threshold voltage which is higher than the threshold voltage of 
My3. Furthermore, the Mp; -Mp» current mirror is dimensioned unequally due to this 
reason. The external terminal Voffset,zcp can be used to change the trigger voltage 
manually for testing purpose. For normal operation, it must be connected to GND. 


5.6.4 Reverse Current Detector 


As mentioned in Sect. 5.4, the NVC allows current flowing in both directions, thus an 
additional reverse current detector (RCD) is required to block the reverse current after 
the transfer process has finished. Without the RCD, the buffer capacitor would be 
discharged after each transfer process. Figure 5.22 shows how the RCD is connected 
on top level. In combination with the switch Ms?, the RCD forms a type of active 
diode as originally presented in [5]. However, the proposed RCD only turns the 
switch off, whereas the switch is turned on by other circuit blocks. In [5], a similar 
comparator turns the switch on and off. 

As shown in Fig. 5.24, the configuration of the RCD is similar to a single-stage 
comparator with a differential stage as input stage, and an inverting amplifier and 
a CMOS inverter as output stage. The functionality of the RCD can be explained 
by means of Fig. 5.23. The RCD is enabled after a zero crossing event of Vz has 
been detected. At that instant, the current /pyr is flowing through Ms» in the positive 
direction charging the buffer capacitor (see Fig. 5.22). Hence, due to the parasitic 
channel resistance of Ms», the Vz— potential at the inverting input is slightly higher 
than Vpur present at the non-inverting input, resulting in a low Vy potential. After two 
inversions are performed by the inverting amplifier and the CMOS amplifier, Drep 
is tied to GND. When /Jpuf approaches zero, the voltage difference (Vr... — Vpbur) also 
approaches zero, and at some instant, Vy flips to high potential, causing Dgcp to 
generate a rising edge. Immediately after that event, Den rcp is pulled low, disabling 
the RCD and causing Dgcp to trigger a falling edge. 
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Fig. 5.23 Simulated transient signals of the RCD 


In order to strictly avoid the occurence of a reverse current, a slight offset is 
introduced to the differential stage by unequal dimensioning of My; and My2. The 
bias voltage for the difference stage and the output stage is generated by the supply 
independent bias generator BG2 (Fig. 5.27), assisting to keep the reverse current 
detection event constant for varying conditions. Compared to the hysteresis com- 
parator used in the NVC control (see Sect. 5.4.2), the biasing current has to be mch 
higher, because the reaction time of the RCD has to be faster. In the current imple- 
mentation, the biasing current is set to 6.3 wA. Hence, the slew rates of the output 
are 36.7 MV /s for the rising edge and 38 MV /s for the falling edge. 

In order to tolerate high piezoelectric voltages, a decoupling stage is inserted 
between the inverting input transistor MN? and Vz—, as shown in Fig. 5.25. The 
decoupling stage shields the low voltage transistor Mw» from high voltages.! 
The node Vy— decoupled 1$ coupled to Vz- only when the RCD is enabled during 
the transfer process, when Vz- is in the range of Vpur. The control signal of Maecoupling 
uses the same level shifting architecture as the bypass in order to keep it turned off 
beyond the transfer process. 


l One could also imagine to replace the standard low voltage transistor Mw» with a high voltage 
transistor. But due to the higher gate-drain capacitance of high voltage devices compared to the 
standard low voltage devices, significant crosstalk is induced into Vy, making it impossible to 
achieve reverse current detection at equal instants for varying supply voltages Vpp and piezoelectric 
voltage amplitudes V». 
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The additional transistors MN4, Mneo, Mn7 and Mp; disable the RCD when it 
is not used, saving valuable power. Furthermore, during the disabled state, they tie 
floating nodes to a defined potential, otherwise false reverse current detections might 
happen after activating the RCD. For testing purposes, an additional variable offset 
can be set externally via the input /otrset RCD- 


5.7 Supply Independent Biasing 


The voltage on the buffer capacitor, Vpbuf, directly depends on the piezoelectric har- 
vester output power. When the environmental vibration is low, the output power is 
also low, resulting in a small buffer voltage. However, for large vibration amplitudes, 
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the piezoelectric harvester delivers lots of power, leading to a large buffer voltage. 
Since the timing of the control circuitry should not depend on the buffer voltage, 
there are two possible scenarios. The first approach uses a DC-DC-converter to pro- 
duce a constant supply voltage which should be around 3 V in order to achieve a 
sufficiently low on-resistance of the MOSFET power switches. Since the tolerated 
supply voltage range is supposed to be 1.5 — 5 V, a step-up and a step-down converter 
or a combined buck-boost converter is necessary. This approach consumes lots of 
power which cannot be harvested and thus is not considered in this work. Another 
scenario is to directly connect the supply voltage to the buffer capacitor and to use a 
supply independent biasing concept to keep the timing constant. This approach will 
be explained in the following. 


5.7.1 Circuit Design 


Two different bias generators BG] and BG2 are implemented, as illustrated in 
Figs. 5.26 and 5.27. Whereas BGI is used to supply the components which have 
to be enabled permanently, BG2 is activated only during the transfer process to sup- 
ply the components which are needed only during that time period. The two bias 
generators are based on the same beta multiplier circuit [6, ch. 11] described in the 
following by means of BGI. The NMOS current mirror Myw;-Mwo copies the current 
lı flowing through the left branch into 7» flowing through the right branch, and the 
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Fig. 5.26 Schematics of the bias generator BG1 
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Fig. 5.27 Schematics of the bias generator BG2 
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Fig. 5.28 Simulated output current of the bias generators as a function of the buffer voltage. a Bias 
generator BGl. b Bias generator BG2 


PMOS current mirror Mp; -Mp» copies J2 back to 71. Since the two PMOS transistors 
have identical dimensions, the currents have to be equal, i.e. 7; = Jy. The resistor R 
defines the current within the beta multiplier—without R, the current is indefinite. 
When the supply voltage Vpp varies, the Vy voltage varies in the same way so that 
the gate-source voltage of Mp; and Mp2 remains constant. 

Figure 5.28 shows that the output current of BG1 and BG2 is not perfectly stable 
as Vbuf 1s swept from 1.5 to 5 V which is due to second order effects like the channel 
length modulation. Compared to the current at Vpuf = 3 V, the deviations at the outer 
edges are —6.1 and 4-9.8 % for BGI, and —4.7 and 4-7.2 % for BG2, which is quite 
acceptable. 
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Since the desired currents are in the nanoampere range, the beta multiplier 1s 
driven in weak inversion. In that regime, the drain current can be written as [1] 


W Vas 
Ip = Ipo — ex — l; 5.9 
7 exp (; 7: ) (5.9) 
where the parameters are 
Ipo a process-dependent parameter, 
w the transistor dimensions, 


n = 1 + £2- the weak inversion slope factor, 
Cox p 
Cp, Cox the depletion and the oxide capacitance, 


Ur = “ar the thermal voltage, 

kp the Boltzmann constant, 
T the temperature and 

q the elementary charge. 


Using (5.9) and W2 = K - Wi, and neglecting the body effect, the Kirchhoff 
voltage loop (KVL) gives 


Vasi1 = Vas? + h R 





= n(2 ) nu (22 ) nur + Dok 
— — } nUr = In | — nUr +I 
Ipo Wi Ipo K - Wi 
U 
> R=- Ink. (5.10) 


In the implemented current sources, K = 6, Ur = 26mV at T = 300K, and 
n = 1.28 which has been extracted from the simulation. Using these values, (5.10) 
evaluates to Jy = 59.6nA for BGI and Jy = 596nA for BG2, which corresponds 
very well to the simulated values 54.7 and 600 nA at Vpu¢ = 3 V. Note that R = 1 MQ 
for BGI and R = 0.1 MQ for BG2. 

The output currents are converted to fixed voltage using diode connected NMOS 
transistors. The fixed voltage generated in BG1, Vstat1, is fed to the components which 
have to be enabled permanently, i.e. the peak detector, the negative voltage converter 
and the SECE selector. Although the ZCD is enabled only during the transfer process, 
the corresponding bias voltage Vstat2 is also generated in BG1. This is due to the small 
output current leading to a quite long time until Vstat2 has settled to a steady state. 
BGI can be disabled externally by setting DgN.cs1 = GND, which is only necessary 
if the external currents Text and exo are injected due to testing purposes. 

The fixed voltages generated in BG2, Vaynı and Vayn2, are used to bias the ZCD and 
the RCD, i.e. the components which are only activated during the transfer process. 
Although Vayn2 is enabled during the whole transfer process, the bias voltage within 
the RCD is released only during the second phase of the transfer process by pulling 
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the DgN.ncp signal to GND, as described in Sect. 5.6.4. This is because simulations 
have shown that the RCD generates undesired pulses during the first phase of the 
transfer process if the RCD is enabled throughout the complete transfer process. BG2 
can be disabled by setting Den, cs2 = GND, which is performed automatically after 
each transfer process until the next transfer process starts. The enable signal can be 
overridden externally for testing purpose in order to inject the external currents Jext1 
and Jext2. 

Generally, the bias voltages for the different components are chosen as small as 
possible, so that the corresponsing blocks still operate properly within the process 
corners. 


5.7.2 Startup 


Both current sources need a startup circuit to establish the desired currents. The 
startup characteristics of BGl works as follows (see Fig. 5.26). Initially, the supply 
voltage Vpp is zero. When the buffer capacitor is being charged, Vpp is increasing 
but Vs is still zero. Thus, at some instant, Mp3 is turned on and precharges the Vx 
node. When Vpp is further increasing, Mp; is activated at a later moment, charging 
Cs and hence turning off Mp3. From that moment on, the beta multiplier current 
can settle to the desired value. It is crucial that Mp3 turns on before Mp; does, 
which is achieved by dimensining Mp3 with a much larger aspect ratio than Mp;. 
The additional transistor Mpg has nothing to do with the startup, it is used only to 
deactivate the beta multiplier externally for testing purpose so that an external current 
can be injected. 

As shown in Fig. 5.27, the startup of BG2 is slightly different from that described 
above. Whereas BGI must be started only once, BG2 is activated at the beginning 
and shut down at at the end of each transfer process. Thus, the beta multiplier current 
is required to settle very fast when it is enabled in order to set the corresponding 
detectors in a ready state. Figure 5.29 depicts the transient behavior of BG2. A pulse 
generator turns on the startup transistor Mp3 for a well-defined duration—not too 
long in order to minimize the cross current through Mp3 and My. In addition, the 
transistor Mp4 releases the Vy node only during the activation time and otherwise 
turns off the output transistors Mpsand Mpg in order to avoid any unnecessary power 
loss. 


5.7.3 Power Consumption Estimation 


The power consumption of both current sources will be estimated in order to find out 
the dominating loss mechanism. Within BGI, significant current is flowing through 
both branches of the beta multiplier plus the two output transistors. The power con- 
sumption of BGI can be calculated as follows: 
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Fig. 5.29 Simulated transient signals of the bias generator BG2 at Vpuf = 3 V 


Pest: = (H1 + D + Icsia + Icsiv) © Vout (5.11) 
Posi = 51 - Vost. (5.12) 


Assuming a constant current /; = 55 nA for all Vpyr, the power consumption 
varies in the range 0.4uW < Pcs; < 1.4 4W for 1.5V < Viuf < 5 V. Since 
BGI is enabled permanently, this power consumption is equal to the average power 
consumption: 

Pogi = Pcsi- (5.13) 


The power consumption of BG2 is given by 


Pes2 = (h + h + Ics2a + Ics2b) - Vout (5.14) 
Poss = 611 = Vout: (5.15) 


Assuming a constant current /; = 600 nA for all Vbuf, the power consumption 
varies in the range 5.4 wW < Pcs? < 18 uW for 1.5 V € Vous < 5 V. This power 
is consumed only during the transfer process duration trp, so the average power 
consumption has to be calculated as 


Pcs2 = 2ftrp + Pcs. (5.16) 
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Assuming f = 175 Hz and a constant trp = 30 us, the average power consump- 
tion varies in the range 0.06 uW < Posy < 02 uW for 1.5 V € Vy < SV. 
Comparing these numbers to Pcs;, it becomes evident that BG1 is clearly dominat- 
ing the total power loss, although BG2 is consuming more power than BGI during 
the transfer process. It has to be mentioned that this power loss estimation is very 
rough and ideal, e.g. the assumption of a constant frp is unprecise since the duration 
of the transfer process is a function of Vbuf. But nevertheless, the conclusion can 
be drawn that especially BG1 and the components biased by it have to be designed 
carefully in order to achieve a low power consumption. 


5.8 PSCE Circuit Top Level 


This section describes the implementation of the PSCE circuit top level. In combi- 
nation with the startup trigger (see Sect. 5.5.1), the buffer capacitor can be charged 
passively through the bypass (see Sect. 5.5.2) if the buffer voltage is too low to supply 
the switch control circuitry. 

Several measures have been taken to ensure stable operation for varying excita- 
tion and load conditions, e.g. buffer voltage, piezoelectric voltage and frequency. The 
oscillation cancellation circuit described in Sect. 5.3 strongly reduces the undesired 
oscillations occuring at the end of the transfer process. A novel ZCD architecture 
shown in Sect. 5.6.3 and the RCD depicted Fig. 5.24 are robust in terms of varying 
excitation parameters. The supply-independent biasing presented in Sect. 5.7 pro- 
vides a biasing for the active components which is almost constant with regard to the 
buffer voltage. High-voltage transistors tolerating 18 V (20 V)? are used in order to 
allow piezoelectric harvesters with a wide output voltage range to be connected. 


5.8.1 Top Level Implementation 


Figure 5.30 shows the top level implementation of the PSCE chip. The piezoelectric 
harvester, the inductor, and the buffer are external components directly connected to 
the chip, everything else is integrated on-chip. The inductance of the inductor has to 
be in the millihenry range since a significant amount of energy has to be transferred 
during each transfer cycle, so it cannot be integrated on-chip. The same is true for 
the buffer capacitor which has to have a capacitance in the microfarad range. 


? The voltage in parentheses is the absolute maximum rating which could destroy the device when 
exceeded. According to the process specification document, proper functionality is only ensured 
when the maximum operating conditions are not exceeded. Operation within the range between the 
maximum and the absolute maximum ratings, i.e. for voltages between 18 and 20 V, may affect 
reliability, especially for exposure to these voltages for extended periods. Measurements have shown 
that the fabricated chips can withstand this maximum voltage without any problems. 
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Fig. 5.30 Energy harvesting system showing the top schematics of the PSCE chip 


The chip is composed of the negative voltage converter (NVC), the power switches 
Ms, Ms» and Ms3, and the switch control circuitry where the control signals for the 
switches are generated. The bypass ensures a conducting path between the piezoelec- 
tric harvester and the buffer capacitor, in case the voltage Vor is too small to supply 
the switch control circuitry. Once Vbuf has reached a certain level, the startup trigger 
gives the signal to deactivate the bypass and start normal operation (Fig. 5.31). 

In order to fulfill the requirements declared in the introduction, high voltage tran- 
sistors with thick oxide, tolerating up to 18 V (20 V) at their terminals have been 
used to realize all the transistors which may come into contact with the piezoelec- 
tric harvester voltage, i.e. the power switches, the NVC transistors and the bypass 
transistors. In contrast to standard 3.3 V transistors which have a minimal length of 
0.35 um, the 18 V high voltage transistors have a minimal length of 3.5 and 2.5 wm 
for the n-channel and p-channel transistors, respectively. In order to obtain compa- 
rable on-resistances, the high voltage devices have therefore a larger area than the 
standard devices and suffer from a higher gate-source-capacitance. Since the gate- 
source-capacitance has to be charged and discharged during each switching cycle, 
this leads inevitably to a higher dynamic power loss. Also, the occupied area is larger 
and thus more expensive. 


5.8.2 Switch Control Circuitry 


Figure 5.31 depicts the block diagram of the switch control circuitry, combining 
the blocks described before such that the control signals for the power switches are 
generated properly. Therefore, mainly standard 3.3 V transistors have been used. For 
those devices, the maximum tolerated gate-source voltage and drain-source voltage 


124 5 Implementation of the PSCE Circuit on Transistor Level 







lottset.PD 
Del | > 
D trig pu D 
P zx A EN, start 
V stat 


Bias Generator T 


BGI 
Zero 
Crossing 
D Detector 
EN,TP V2 








Vottset.ZCD | > D " Switch 
Bias Generator | Control 


BG2 Logic 


Reverse 
Current 
Detector 






SECE 
Selector 


V,, D 


DEN, sece13 O 


Fig. 5.31 Switch control circuitry of the PSCE chip (unless otherwise indicated, the supply rails 
of all blocks are Vpp and GND) 


are both 3.6 V (5 V). Hence, since the supply rail of the switch control circuitry, 
Vpp, is directly connected to the buffer capacitor, the buffer capacitor voltage is 
limited to 5 V. An additional external load resistor is necessary to restrain the buffer 
capacitor voltage to the upper limit, because there is no over voltage protection circuit 
implemented in the chip. 

The corresponding logic is shown separately in Fig. 5.32. Some RS latches and 
multiplexers combine the output signals of the detectors in order to generate the 
control signals for the power switches, as explained in the following. As described 
in Sect. 5.6.2, the startup trigger signal Dyrig initializes the peak detector for proper 
operation directly after the startup has occured. During the open circuit phase, the 
peak detector is waiting for a peak of the rectified piezoelectric voltage VL+, whereas 
the ZCD and the RCD are disabled. Ds; and Ds3 are low, turning off the NMOS 
switches Ms; and Ms3. As discussed in Sect. 5.1.3, the level shifter clamps Ds» to 
max (VL+, Vbuf) in order to turn off the PMOS Ms» properly. When a peak of VL+ 
is detected, the SECE selector output level decides whether SECE13 or SECE1323 
is activated. 

When the piezoelectric harvester excitation is low so that VL+ < 2Vpu¢ holds, the 
SECE selector outputs a logic high signal (Dsgcg = Vor). So, the multiplexers L3 
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Fig. 5.32 Switch control logic block of the PSCE chip (unless otherwise indicated, the supply rails 
of all blocks are Vpp and GND) 


and L4 feed through the lower input signals, leading to a SECE13 configuration. This 
means that when a peak of Vr; , is detected, Ds, is set high, triggering the energy 
transfer from the piezoelectric harvester into the inductor by turning on Msı. At the 
same time, the ZCD is activated which is now waiting for VL+ to reach 0 V. When 
this happens, the ZCD outputs a short pulse, resetting L;'s output to GND, turning 
off Ms; and the ZCD. At the same time, L» sets its output to high level, causing 
Ds3 to perform a rising edge and Ds» to perform a falling edge, turning on Ms? and 
Mss. Thus, the energy stored on the inductor is released to the buffer capacitor. At 
the same time, the RCD is activated in order to monitor the voltage drop across M59. 
When the current flowing into the buffer capacitor starts to get negative, the voltage 
drop also gets negative, triggering a Dgcp pulse and thus turning off Ms» and Ms. 
At that instant, the energy transfer from the piezoelectric harvester into the buffer 
capacitor has finished, and the peak detector is again waiting for a new peak of VL+. 

If VL+ exceeds 2 Vpuf during the open circuit phase, Dsgce is pulled low, activating 
SECE1323. In this case, the multiplexers Lọ and L5 feed through the upper input 
signals, resulting in a Ds; which is permanently low. When a peak of Vz is detected, 
RS latch Ls pulls it output to high level, resulting in a low level of Ds». At the 
same time, Den Tp is set high and Den. TP is set low, deactivating the oscillation 
cancellation transistor Moc (see Fig. 5.30) and activating BG2. In addition, the ZCD 
is enabled by setting Den,zcp high. When a zero crossing of Vr, is detected, the 
ZCD is disabled, and the RCD is enabled by setting Den rcp high by means of the 
RS latch L2. So, Ds3 is set high as well, initiating the second transfer process phase 
(Ds» remains low). When a reverse current is detected, Ds3 is pulled low through L2, 
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and Ds» is set high through Ls and the inverter Iı, terminating the transfer process 
and deactivating the RCD and BG2. At the same time, Den. TP is set high, turning 
on Moc and thus engaging the oscillation cancellation, and Dsgczg is set high in 
order to reset the switch control circuitry. In order to compare the SECE1323 and the 
SECE13 technique under the same conditions, the SECE13 technique can be forced 
externally by setting Den sEcE13 high. In that case, Dspcg is permanently high. 


5.8.3 Die and Packaging 


Figure 5.33a shows the micrograph of the PSCE chip which has been fabricated in 
a 0.35 um process. The total active area is 1.25 mm^. The NVC MOSFETs occupy 
0.55 mm, the control circuitry needs 0.36 mmĉ, and the power switch MOSFETs 
cover 0.34 mm^. In total, the die comprises 44 pads, of which 7 are essential: Two 
pads for the piezoelectric generator (Vp_, Vp. .), two pads for the inductor (Vr, 
Vr), one pad for the buffer capacitor (Vor), one pad for GND, and one pad for 
activating SECEI323 or SECEI3. The other pads are for monitoring several digital 
and analog voltages, and for applying offset on the detectors. 

Figure 5.33b depicts a photograph of the chip placed into a CLCC48 package 
having 48 pins. All the pads are bonded to the corresponding package pins. Mounted 
into a special socket, measurements can now be performed. 
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Fig. 5.33 PSCEchip photographs. a Die micrograph. b Photograph of the bonded die in a CLCC48 
package 
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Chapter 6 
Performance Analysis of the PSCE Chip 


The Pulsed Resonant Charge Extractor (PSCE) chip has been characterized thor- 
oughly in order to verify the theoretical discussion of the operation principle. This 
chapter shows the corresponding measurement results, which have been gained using 
real piezoelectric energy harvesters under laboratory conditions. A laboratory shaker 
table has been used as a vibration source, in order to guarantee reproducible mea- 
surement conditions. After the corresponding setup has been explained, the two 
piezoelectric energy harvesters used as a power source are characterized with respect 
to different load circuits, and the procedure of extracting the values of the electri- 
cal equivalent circuit parameters is explained step by step. Then, a demonstration 
platform showing the effectiveness of the PSCE chip by means of a blinking light 
emitting diode (LED) is described. 

The main focus of this chapter is the performance analysis of the PSCE chip, 
describing in detail how the output power and the overall chip losses can be measured. 
Simulations are used to extensively analyze the different loss mechanisms within the 
chip. After that, the operation limits of the PSCE chip are discussed, ending up 
with the conclusion that the presented chip can be used universally in a wide range 
of environment conditions. The chapter ends with an analysis of the cold-startup 
capability. 


6.1 Vibration Setup 


Figure 6.1 shows aschematic view of the vibration setup which is used to characterize 
the PSCE chips. The piezoelectric harvester is mounted on a shaker (S 51110 from 
TIRA GmbH), which is excited by a sine wave from a signal generator amplified 
through a power amplifier (TIRA BAA 120). The output signal from an acceleration 
sensor (4513-B-001 from Brüel and Kjær) attached to the vibrating shaker plate is 
fed to the controlling unit (VibrationView® VR8500 from Vibration Research Cor- 
poration) comprising the signal generator and a controller. Using the Vibration View 
software running on a personal computer connected to the controlling unit, amongst 
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Fig. 6.1 Schematic view of the vibration setup 
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Fig. 6.2 Photograph of the vibration setup 


many others, the vibration frequency and the acceleration amplitude can be set. The 
controller regulates the frequency and the amplitude of the integrated signal generator 
such that the acceleration amplitude of the shaker plate remains constant, indepen- 
dent of any forces acting on the shaker plate. So, it is possible to connect different 
electrical load circuits to the piezoelectric harvester, without affecting the vibration 
amplitude due to their electromechanical feedback, as described in Sect. 2.4.2. A 
second accelerometer is attached to the shaker plate in order to visually check the 
proper value of the vibration amplitude. Special care has to be taken regarding the 
the fixture of the piezoelectric harvester and the accelerometer. Everything has to 
be properly fixed by screws, since other types of mounting, e.g. double-side tape, 
dissipates energy and thus causes the actual acceleration to deviate from the specified 
value. 

Figure 6.2 shows a photograph of the vibration setup while the shaker is running, 
and the PSCE chip is operating. The oscilloscope depicts the V; ,. signal (yellow) 
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Fig. 6.3 Eichhorn harvester. a Schematic view [2] (reproduced with kind permission of Christoph 
Eichhorn). b Photograph 


and the output signal of the second accelerometer (green). The digital multimeter 
displays the current buffer voltage. 


6.2 Characterization of Piezoelectric Harvesters 


Since the scope of the fabricated PSCE chip is the application in systems where vibra- 
tions are present, the measurements are mainly performed using real piezoelectric 
harvesting devices. For this work, one commercially available harvester (V22B from 
Mide Technologies, Inc. [4]) and one custom-made harvester (Eichhorn harvester) 
developed by Christoph Eichhorn [2] have been used. In this section, these two har- 
vesters will be evaluated regarding their output power performance with different 
load circuits. 


6.2.1 Eichhorn Harvester 


Figure 6.3a shows the schematic view of the Eichhorn harvester. Three layers of PZT 
are bonded on top of each other. Whereas the top and the bottom layers are used for 
the harvesting process, the middle layer only acts as a carrier. The harvesting layers 
are connected for parallel operation, so the output current is doubled compared to a 
single piezoelectric layer whereas the output voltage remains constant. As described 
in [2], the original intention of this harvester is the ability to tune its resonance 
frequency to the vibration frequency the harvester is exposed to. Therefore, the 
two outer arms can be used as piezoelectric actuators, inducing an axial preload 
on the main beam which is the actual harvester. Here, the harvester 1s used without 
triggering the frequency tuning mechanism, so the actuator arms simply are additional 
masses without any functionality. The photograph of a fabricated Eichhorn harvester 
is depicted in Fig. 6.3b. 
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Fig. 6.4 Measured open circuit amplitude of the Eichhorn harvester as a function of the frequency 
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Fig. 6.5 Measured harvested power using the Eichhorn harvester for a resistive load and a diode 
rectifier load as a function of the load resistance 


The collaboration with Christoph Eichhorn developed during the common time 
in the graduate program Micro Energy Harvesting. Using the Eichhorn harvester as 
a power source for an early version of the PSCE chip led to good results regarding 
the output power improvement compared to a common full bridge Schottky diode 
rectifier. So, the two projects were combined, resulting in a common demonstrator 
described in Sect. 6.3. 

Figure 6.4 shows the open circuit amplitude of the Eichhorn harvester as a function 
of the frequency. The shaker was excited with an acceleration amplitude of 5 m/s?. 
The curve shows a very symmetrical shape around the resonance frequency which is 
118 Hz, therefore the harvester 1s assumed to show linear behavior. At the resonance 
frequency, a maximum amplitude of 3.1 V is achieved. 

According to Fig. 6.5, the Eichhorn harvester dissipates a maximum output power 
of 50uW at a resistive load. Using a full bridge rectifier composed of BAT 46 
Schottky diodes, the optimum harvested power is 30 uW. 
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Fig. 6.6 Mide V22B piezoelectric harvester. a Schematic view [4] (dimensions given in inch, 
reproduced with kind permission of Midé Technology Inc.). b Photograph of the harvester including 
clamping 


6.2.2 Mide Harvester 


Due to the custom design and the brittle nature of the PZT material, the maximum 
tolerated deflection of the Eichhorn harvester beam is limited, and thus the output 
voltage is limited to a few volts. Since the PSCE chip is optimized for input voltages 
up to 20 V, it became evident that another harvester had to be searched for, which 
is more robust and allows a higher tip displacement. So, the commercially available 
Mide V22B harvester [4] schematically shown in Fig. 6.6a was chosen. The harvester 
is made of two piezoelectric wafers inside a plastics compound. As the wafers are 
electrically isolated against each other, they can be either connected in parallel or 
in series externally. Compared to a single wafer and assuming the same excitation 
conditions, the series connection implies half the electrical output capacitance, the 
same output current and twice the output voltage. For the parallel connection, one 
gets twice the capacitance, twice the current, and the same voltage. Since in most 
applications it is desirable to have a high output current, the parallel connection is 
chosen for the measurements shown in this book. 

Due to the fact that piezoelectric material is not very elastic, the tip deflection 
is not very high, and thus the output power is limited. So, a small additional mass 
of 1 x 10^? kg is attached to the tip of the cantilever. Here, the tip mass is made 
of two permanent magnets made of neodymium iron boron (NdFeB) with a mass 
of 0.5 x 10-7? kg each, in order to be capable to replace the mass if necessary. 
Figure 6.6b shows a photograph of the Mide V22B harvester including clamping. Its 
parameters are determined in Sect. 6.2.3. 


6.2.3 Determination of Piezoelectric Parameters 


In this section, it is described how the the values of the lumped elements of the 
piezoelectric model depicted in Fig.2.12 are determined. This procedure is only 
applied to the Mide V22B harvester, because the Eichhorn harvester is exclusively 
used for the demonstrator as discussed in Sect. 6.3 and thus will not to be evaluated 
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in detail. For the consideration described in the following, assume the vibration setup 
as described in Sect. 6.1, and the shaker acceleration amplitude set to à = 10 m/s”. 

In the literature, some publications can be found which present methods to deter- 
mine the lumped element parameters by means of finite element analysis (FEA) 
[12], or by means of the properties of the involved piezoelectric material [9]. The 
first method implies knowledge in FEA simulation tools and requires that the shape 
of the more or less complex piezoelectric harvester is implemented in such a tool, 
which can be very time-consuming. The other procedure depends on the knowledge 
of the exact piezoelectric parameters, which is problematic since datasheet values 
are very often only valid for a certain operation condition or sometimes even wrong. 

Hence, another method which has been demonstrated partly in [7] and [1] is used 
here to extract the lumped element parameters from a given piezoelectric harvester. 
The total effective mass is determined as 


m = Mbeam,eff + Mtip = 140 "beam + Mtip, (6.1) 
where Mpeam 1S the beam mass, mMpeam,eff 1S the effective beam mass [5], and Mtp is 
the tip mass. The deflection of a one-side clamped beam is not uniformly distributed 
across the beam: whereas the end of the beam has a large deflection, the displacement 
of the regions near the clamping is very small. Thus, the effective beam mass is much 
smaller than the actual beam mass. Since it 1s difficult to measure the beam mass 
of a piezoelectric harvester including clamping structure, the beam mass should 
be calculated from its geometry and the density of the involved materials. For the 
Mide V22B transducer, a tip mass of mip = 1.0 x 107° kg is taken. The beam is 
determined as myeag = 1.0 x 107? kg. Thus, according to (6.1), the total effective 
mass is m = 1.236 x 107° kg. 

Then, the open circuit piezoelectric voltage amplitude is measured as a function 
of the frequency, and the same is performed using a 1 KQ resistance which can 
be considered a short circuit, because this is much smaller than the optimal load 
resistance, as can be seen later. As shown in Fig. 6.7, the open circuit resonance 
frequenciy woc = foc - 2zt and the short circuit resonance frequency «xc = fse - 27t 
are the frequencies at which the corresponding voltages are maximum. The squared 
effective coupling factor k2. can now be calculated according to (2.28): 
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Fig. 6.8 Measured output voltage of the Mide V22B harvester in open circuit (Rz = 10 MQ) and 
short (Ry, = 1 KQ) circuit, versus the operating frequency 


Fig. 6.9 Half-power band- V, 
with method to determine the 
quality factor of the mechani- 
cal resonator (schematic view) 





kore = E. (6.2) 


As shown in Fig. 6.8, the open circuit and the short circuit resonance frequency 
are determined as Woe = 174.0 Hz - 2z and oy. = 172.3 Hz - 27 for the Mide 
V22B harvester, and from that, the squared effective coupling factor is calculated as 
LS 0000. 

From the angular short circuit resonance frequency, the stiffness of the mechanical 
resonator including the piezoelectric layer, k, is determined using (2.29): 


For the Mide V22B harvester, the stiffness is k = 1,447.8 N/m. 

The open circuit voltage curve can be used to determine the quality factor of the 
mechanical resonator. According to Fig. 6.9, the half-power bandwidth œ — @, = 
27 fo — 27x fı is the bandwith at the point where the voltage drops to Vp 34B = 
VP max: 1/ Jin VP max ‘0.707 compared to the maximum voltage, which represents 
a drop to approximately —3 dB. Due to P = V?/R, the power drops by 1/2 at that 
point compared to the maximum power. The quality factor can now be calculated 
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according to 
Woc 


C=. (6.3) 


o — @1 


The quality factor is linked with the damping, the seismic mass and the short- 
circuit resonance frequency [8, 10] by the equation 


In Osce 


Q 


For the Mide V22B harvester, the quality factor is determined as Q = 46.8, and so 
the damping can be calculated as d = 28.6 x 107? kg/s. 

When the short circuit current amplitude ioe and the open circuit voltage ampli- 
tude VP oc at their respective resonance frequency are known, the piezoelectric capac- 
itance can be found by 


d= 





(6.4) 


IP sc 


Cp = (6.5) 


D , 
Wsc VP.oc 


under the assumption that the imaginary part of the piezoelectric harvester impedance 
dominates the real part, i.e. the internal impedance of the piezoelectric harvester is 
purely capacitive [6]. From the measured short circuit current amplitude [T E 
0.290 mA and open circuit current amplitude VP oc = 13.93 V, the electrical piezo- 
electric capacitance is calculated as Cp = 19.0 nF. Finally, the remaining parameter 
I' can be found by means of (2.30): 


I' = kett/ KC p, 


which leads to I? = 0.745 x 107? As/m from the values calculated above. The 
simulated output power using the resistive load fits better to the measured value 
when the piezoelectric capacitance value is slightly increased to Cp = 19.5 nF, 
while I” is kept at the calculated value. Table 6.1 summarizes all the values of the 
Mide V22B harvester derived above. 

As a final step, the lumped element parameters of the piezoelectric model shown 
in Fig. 2.12 can now be calculated using (2.35) as 


Ving = 16.6V-sin(@t) Lyc = 2.22KH Rye = 51.5kK22 Cyc = 0.383 nF. 
(6.6) 
In order to validate the behavior of the simulation model, some characteristic 
simulation curves should be performed and compared to the measured ones. The 
open circuit and short circuit frequency sweeps are important to check whether the 
open circuit and short circuit resonance frequency match. Moreover, the maximum 
open circuit voltage amplitude and short circuit current amplitude should be equal 
as well. Due to non-linearities of the real devices (e.g. a non-linear spring constant) 
which are not implemented in the model, it is possible that the voltages/currents do 
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Table 6.1 Parameters of the 


Parameter (unit) Value 

Mide V22B piezoelectric — ; 

harvester with additional 1 g Acceleration amplitude @ (m/s^) 10 

tip mass Total effective mass m (kg) 1.236 x 107? 
Piezoelectric capacitance Cp (F) 19.5 x 107? 
Coupling coefficient J” (As/m) 0.745 x 107? 
Stiffness k (N/m) 1,447.8 
Damping d (Ns/m) 28.6 x 107? 
Quality factor Q () 46.8 


Squared effective coupling coefficient ae O| 0.020 
Maximum extractable power Piim (UW) 667.4 


not perfectly match for frequencies different from the resonance frequency. As a final 
verification indicator, it is important that the root mean square power extracted with 
the optimal load resistor in measurement is in good agreement with the simulation. 

In case some of the parameters mentioned in the previous paragraph differ strongly 
between simulation and measurement, adjustments are necessary. If the simulated 
maximum voltage amplitude and the extracted power are both too small compared 
to the measured values, the effective mass should be increased. In case the voltage 
matches but the extracted power is too low in simulation, the piezoelectric capacitance 
must be increased. After each adjustment step, it is mandatory to recalculate all 
parameters that depend on the changed parameter. 


6.2.4 Evaluation of the Model 


In this section, itis discussed to what extent the parameters extracted in Sect. 6.2.3 and 
summarized in (6.6) represent the real Mide V22B harvester correctly. As indicated in 
Table 6.1, an excitation acceleration of à = 10 m/s? is used throughout this section. 

Due to the electromechanical feedback, the resonance frequency changes with 
alternating load resistance, so in order to get perfect results, the resonance frequency 
could be searched for each load resistance. Since this method represents an immense 
measuring effort which is error prone, it is common practise to set the excitation 
frequency to either the short circuit resonance frequency [8] or the open circuit 
resonance frequency [11], accepting small deviations from theory. Here, the corre- 
sponding open circuit resonance frequency is chosen as an operating point for the 
simulation, the measurement and the calculation. 


6.2.4.1 Resonance Characteristics 


In order to validate the behavior of the simulation model, first of all the open circuit 
and the short circuit voltage amplitude of the piezoelectric harvester are simulated as 
a function of the excitation frequency and compared to the measured characteristics. 
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Fig. 6.10 Simulated and measured piezoelectric voltage amplitude as a function of the excitation 
frequency. a Open circuit (Ry, = 10 MQ). b Short circuit (Ry, = 1 KQ) 


Therefore, the coupled harvester model as shown e.g. in Fig. 6.12a and as discussed 
in Sect. 2.4.3 has been used, where the values of the piezoelectric model parameters 
are taken from Table 6.1, in order to represent the Mide V22B harvester. For the open 
circuit case, a load resistance Ry, = 10 MQ is attached to the piezoelectric terminals 
in order to simulate the oscilloscope probe. The 1 kQ resistor represents the short 
circuit condition, as written in Sect. 6.2.3. The measurements have been performed 
using the setup described in Sect. 6.1 with a shaker acceleration amplitude set to 
10 m/s?. The frequency has been swept from 160 to 190Hz using the Vibration View 
controller, maintaining the constant vibration amplitude for each frequency. 

The resulting comparison plots are demonstrated in Fig. 6.10. It 1s obvious that 
the corresponding resonant frequencies in open circuit, 174 Hz, and in short circuit, 
172.3 Hz, fit to the measured characteristics very accurately. Whereas the simulated 
curves exhibit a very symmetric shape around the resonant frequencies, the mea- 
sured ones look somehow distorted. The reason for this behavior is the occurrence of 
nonlinearities in the real harvester, probably the mechanical and/or the piezoelectric 
stiffness may contain nonlinear terms which are not considered in the simulation 
model. Therefore, the simulated open circuit and short circuit voltage amplitudes 
in respective resonance are slightly larger than the measured values. But since the 
optimum output power using the resistive load, which is the more important para- 
meter, matches almost perfectly between simulation and measurement as shown in 
Sect. 6.2.4.2, the deviations in the voltages are acceptable. 

The resonance characteristics using the SECE load is also performed by measure- 
ment and simulation, the according setups are described in Sect. 6.2.4.4. Figure 6.11 
shows the simulated and measured characteristics of the piezoelectric voltage ampli- 
tude for connected SECE load. Similar to the open circuit and short circuit load 
discussed before, the measured voltage is slightly lower than the simulated one, and 
the measured curve looks somehow distorted compared to the nice-looking sym- 
metric simulated one. But whereas the simulated resonance frequency has slightly 
decreased by 0.25-173.75 Hz compared to the open circuit load, the measured reso- 
nance frequency has increased significantly by 1.6—175.6 Hz. Due to the fact that the 
SECE technique drives the piezoelectric harvester most of the time in open circuit, 
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Fig. 6.11 Simulated and measured piezoelectric voltage amplitude as a function of the excitation 
frequency, for SECE load 
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Fig. 6.12 Setups for determination of the output power using the resistive load. a Simulation setup. 
b Measurement setup 


the simulated characteristics meets the expectations quite well, whereas the mea- 
sured behavior cannot be explained by means of the linear model taken as a basis 
for the considerations used in this chapter and throughout the book. For a deeper 
analysis of this phenomenon, a nonlinear model of the piezoelectric harvester should 
be derived, but since the modelling was not the main target of this work, this issue 
is not further pursued here. 


6.2.4.2 Resistive Load 


Due to its simplicity, the resistive load is the interface circuit best suited to evaluate 
the correctness of the piezoelectric model. For the simulation, the configuration 
shown in Fig.6.12a is used. The coupled model of the piezoelectric harvester as 
discussed in Sect. 2.4.3 incorporating the values from Table 6.1 is taken in order to 
consider the electromechanical feedback on the harvester. Figure 6.12b shows the 
measurement setup which is used to determine the output power. An oscilloscope is 
taken to measure the RMS power dissipated in the load resistance. The probe with a 
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Fig. 6.13 Comparison of the analytical, simulated and measured output power using the resistive 
load versus the load resistor at resonance (f = 174 Hz) 


high input impedance of 10 MQ does not affect the actual load resistor unless its value 
is close to the probe impedance. By means of a resistor decade, the load resistance 
is swept in the range 1 kQ-10 MQ, resulting in an effective range 1 k$22-5 MQ due 
to the probe impedance. 

Figure 6.13 compares the analytical, the simulated and the measured output power 
as a function of the load resistance. The maximum analytical power can be determined 
by means of (3.5), whereas Ry, = Ryes opt is assumed which is first calculated by 
means of (3.6). From these equations, the maximum output power using the resistive 
load is 533 uW at an optimal load resistance of 64k. Note that this value is lower 
than Piim = 667 pW as stated in Table 6.1. This is because Pj can only be achieved 
at a single value of the squared coupling factor, p as given in (3.7). Obviously, 
the actual coupling factor of the Mide harvester is different from this optimal value. 

As shown in Fig. 6.13, all the three curves match almost perfectly. The simulated 
curve exhibits an optimum power of 530 uW at an optimal load resistance of 65 kQ 
for the simulation. For the measurement, the same optimum power is achieved for 
an optimum load resistance of 50 kQ. This is a very small deviation which can be 
justified by slight imperfections in the measurement setup. 


6.2.4.3 Full-Wave Rectifier with Capacitor Load 


In order to evaluate the piezoelectric model using the full-wave rectifier, first of all it 
has to be defined which type of full-wave rectifier is best suited. Therefore, the active 
rectifier, the NVC rectifier and the diode rectifier are simulated using the setup shown 
in Fig. 6.14a. The active rectifier uses the active NVC (see Fig. 5.10a) in combination 
with a PMOS switch in order to block the reverse current. All involved MOSFETs 
are actively driven by a VerilogA control circuitry which consumes no power, in 
order to make the corresponding performances comparable to each other. The NVC 
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Fig. 6.14 Setups for determination of the output power using the full-wave rectifier load (f = 
174 Hz). a Simulation setup. b Measurement setup 


rectifier is composed of the passive NVC (see Fig. 5.10b) and a Schottky diode. The 
multiplier of the Schottky diode is set to 150, in order to achieve a voltage drop of 
200 mV which corresponds to the measured voltage drop across a BAT 46 Schottky 
diode. The same number of parallel Schottky diodes is used for the diode bridge 
rectifier. 

Figure 6.15 depicts the simulated output power curves for the three different rec- 
tifiers, for operation at resonance and at far off-resonance. For the operation at reso- 
nance, the optimum output power using the active rectifier is 393 uW, representing 
the highest value compared to the other rectifier topologies. By means of the NVC 
rectifier and the diode rectifier, 98.7 and 94.4 96 of the highest value can be harvested. 
The NVC rectifier performs better than the diode rectifier, because it suffers from 
only one diode voltage drop instead of two. For the operation at far off-resonance, 
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Fig. 6.15 Comparison of the simulated output power for the active rectifier, the NVC rectifier and 
the diode rectifier. a Resonance. b Far off-resonance 


the optimum output power using the active rectifier is 8.8 uW, again representing 
the highest value compared to the other rectifier topologies. By means of the NVC 
rectifier and the diode rectifier, 89.7 and 77.3 96 of the highest value can be harvested, 
which is less compared to the operation at resonance. This discrepancy is due to the 
fact that the power dissipated by the Schottky diode is playing a more dominant role 
if less power is harvested. 

From the simulation results described above, the active rectifier represents the full- 
wave rectifier having the best performance. However, the implementation is much 
more difficult in reality compared to the other rectifier types. Hence, the second best 
architecture which is the NVC rectifier is chosen to serve as a full-wave rectifier for 
the measurements. 

Figure 6.14b illustrates the measurement setup of the NVC rectifier, which 1s 
established simply by connecting one BAT 46 Schottky diode having a low voltage 
drop of 200 mV to the PSCE chip at the Vz, terminal. A buffer capacitor with 
Couf = 22 pF smoothes the voltage ripple so that a standard digital multimeter can 
be used to measure the steady state buffer voltage Vbuf ss. The buffer voltage settles 
to a stable value shortly after the load resistance has been set to a new value, ending 
up in a state where the power removed from the buffer capacitor by the load resistor 
equals the power flowing into the buffer capacitor, 1.e. the harvested power, which 
can be calculated using 


2 


V 
iS (6.7) 
Rr 


Similar to the resistive load shown in Sect. 6.2.4.2, the excitation frequency has 
been set to 174 Hz, corresponding to the open circuit resonance frequency of the 
piezoelectric harvester. 

Figure 6.16 depicts the comparison of the analytical, the simulated and the mea- 
sured output power using the full-wave NVC rectifier. According to (3.20) and (3.18), 
the optimum achievable output power is calculated as 465 uW at a load resistance of 
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74 k&2. The discrepancy between the analytical and the two other curves is due to the 
fact that the analytical equation (3.15) does not consider the losses in the rectifier. 
The measured and the simulated curve match almost perfectly. Their optimum load 
resistance is 100 kQ, and the corresponding output power values are 390 uW for 
measurement and 400 uW for simulation. 


6.2.4.4 SECE Interface Circuit 


In order to evaluate the harvester characteristics using the SECE load without com- 
plex loss mechanisms, a simplified simulation model of the SECE interface circuit 
has been generated, as depicted in Fig.6.17a. The control signals for the switch 
and the NVC MOSFETs are generated by a control block programmed in the hard- 
ware description language VerilogA. Thus, the simulation time is drastically reduced 
compared to a simulation using a transistor level design. In addition, the power con- 
sumption of the control block is omitted, because its influence is not a matter of 
concern here, and any other disturbing influences are ignored. The high level of the 
control signals is set to 30 V in order to minimize the channel resistance of the MOS- 
FETs. The harvested power is determined by averaging the output power product, 
Pharv (£) = Vouf /but (t), over several periods after reaching the steady state. 

For the measurement, the PSCE chip as evaluated in detail in Sect.5.8 is taken, 
leading to the configuration shown in Fig. 6.17b. A Panasonic inductor with an induc- 
tance of 10 mH and a parasitic wiring resistance of 4.3 Q is connected (see Table 6.5), 
leading to the highest efficiency. The same values are used in the simulation model. 
SECE23 is enabled to reduce the conduction losses and thus to come closer to the 
theoretically achievable power. Unfortunately, it is not possible to supply the control 
circuitry with an external power source, so especially the conduction losses in the 
parasitic series resistances of the MOSFET channels, the inductor etc. are high for 
low buffer voltage values, leading to a drop in the harvested power. But also the 


144 6 Performance Analysis of the PSCE Chip 


(a) 


Piezoelectric Harvester SECE Interface 


—————ÁÓÁ————————— ———M——ÓM————————————— 


| C E 10mH 4.3Q 


a eee EREE ee ee ee et 


VerilogA Control 


(b) Off-chip 


Panasonic Inductor Digital 
pe en em Multimeter 





Fig. 6.17 Measurement and simulation setup for determination of the output power using the SECE 
interface circuit. a Simulation setup (all MOSFETs are 1,000 jum/0.35 jum, and the control signal 
high level is 30 V). b Measurement setup 


power consumption of the control circuitry and the losses in the gate drivers of the 
switch MOSFETs are considered, all representing undesired losses for this analysis. 

Similar to the rectifier load, the measured harvested power curve is determined 
varying the load resistance, and calculating each single value by means of (6.7). But 
since it is more meaningful to know the buffer voltage instead of the load resistance 
in case of the SECE interface, the output power is given in terms of the buffer voltage. 
This makes it more appropriate to vary the buffer voltage directly in the simulation, 
instead of varying the load resistance. In order to drive the piezoelectric harvester in 
resonance, different excitation frequencies are set for simulation and measurement, 
due to the shifting of the resonance frequency observed in measurements as discussed 
in Sect. 6.2.4.1. 

Figure 6.18 depicts the calculated, measured and simulated output power curves 
using the SECE interface. According to (3.29), the theoretical output power does 
not consider any losses and is independent of the load, so it represents a straight 
horizontal line at 664 uW. The simulated curve shows a maximum harvested power 
of 615, W at Viuf = 5 V. Due to the conduction losses present in the simulation 
model, the corresponding curve does not achieve the theoretically achievable power. 
These losses are also responsible for the sharp drop for buffer voltages tending 
towards zero, corresponding to the theoretical analysis performed in Sect. 4.5.3.2. 
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Table 6.2 Overview of the output power values (given in pW) at the respective optimum resistance 
of the Mide V22B harvester for different load circuits 


Full-wave rectifier SECE 


Analytical 2323 465 664 





As expected, the measured performance is lower due to the significant chip losses, 
reducing the maximum harvested power to 477 uW at Viuf = 5 V. 


6.2.4.5 Conclusion 


In this section, the power which can be harvested from a piezoelectric harvester 
using the resistive load and the rectifier-capacitor load has been discussed by means 
of analytical calculation, simulation and measurement. For the analytical calculation 
and the simulation, the parameters summarized in Table6.1 have been extracted 
from the Mide V22B harvester using the procedure described in Sect. 6.2.3. Table 6.2 
summarizes the output power values at the respective optimum resistances. Due to 
the absence of additional loss mechanisms, the resistive load is the most appropriate 
interface circuit to evaluate the correctness of the piezoelectric harvester model as 
described in Sect. 2.4. Since the optimum output power values especially for the 
resistive load, but also for the rectifier load, are in very good agreement, it can be 
concluded that the discussed harvester model with the extracted values is well suited 
to model the behavior of the real piezoelectric harvester device—at least for excitation 
at resonance. Due to the nonlinear characteristics of the real device, the excitation 
frequency differs in measurement and simulation for operation at off-resonance. 
Since the design of the SECE load is more complex compared to the others, 
incorporating more complex loss mechanisms, the SECE load cannot serve as an 
indicator for the quality of the harvester model. But nevertheless, the comparison 
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Fig. 6.19 PCB utilized for measurements with the PSCE chip 


of the theoretical/simulated with the measured curves allow a statement about the 
chip losses without summing up each loss component, as performed in Sect. 6.4.3.1. 
Also, due to the nonlinearity of the harvester device, the SECE load does result in a 
frequency shift which cannot be predicted by means of simulation. 


6.3 Demonstration Platform 


An early version of the PSCE chip which implemented the SECE13 technique was 
used in combination with the Eichhorn harvester (see Sect. 6.2.1), in order to set up 
a demonstration platform, showing the usefulness of the PSCE chip in contrast to 
standard diode rectifiers [3]. First, the performance of that PSCE chip is shown, and 
then the demonstrator board using that chip is described. Note: Whenever “PSCE 
chip” is mentioned in this section, the early version is meant. 

For the performance analysis, the shaker vibration amplitude has been set to 
5 m/s”, and the vibration frequency is tuned to 118 Hz corresponding to the resonance 
frequency of the Eichhorn harvester (see Sect. 6.2.1). 

Figure 6.19 shows the PCB used for measurements with the PSCE chip, which is 
placed in a CLCC68 socket to establish conection to the outer world. The two-wire 
piezoelectric harvester output is connected via the corresponding connector pins. A 
2.2 pF aluminum electrolyte capacitor serves as the buffering element, a load resistor 
can be connected in parallel. Four different inductors are mounted which are jumper- 
selectable. A multi-pin connector is used to monitor the digital outputs. For testing 
purpose, external voltages can be applied by several banana plugs, e.g. in order to 
supply the switch drivers externally. 

Table 6.3 summarizes the specifications of the inductors used in combination with 
the PSCE chip. 
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Table 6.3 Inductors used in combination with the PSCE chip 
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Fig. 6.20 Measured transient characteristics using the PSCE chip (Vbuf = 1.9 V, L = 2.2 mH). 
a VL+ over several periods. b Analog and digital signals during the transfer process 


Figure 6.20 shows the NVC voltage/rectified piezoelectric voltage V; ,. in steady 
state over several cycles, assuming a constant buffer voltage Ving = 1.9 V. The 
characteristics of the piezoelectric voltage match the theoretically predicted char- 
acteristics as illustrated in Fig.3.7. The transfer happens at the instants when the 
piezoelectric voltage peaks, resulting in a rapid decrease decay of the piezoelec- 
tric voltage. The rest of the time, the piezoelectric harvester is driven at open 
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Fig. 6.21 Measured output power using the PSCE chip with different inductors, compared to the 
output power using the Schottky diode rectifier (L — 2.2 mH) 


resonance. The amplitude Vp — 4.3 V overrides the open circuit amplitude by 1.2 V 
(see Sect. 6.2.1), illustrating the improvement achieved by the SECE technique. The 
small bumps after the transfer process indicate that the transfer process is initiated 
slightly prematurely, which is due to the vibration frequency which is below the 
frequency the peak detector was designed for (see Sect. 5.6.2). 

According to Fig. 6.20b showing the NVC voltage V;, the inductor current 7/7, 
and the digital control signals during the transfer process, the principle operation 
works properly. The peak detector pulse initiaties phase 1 of the transfer process, 
causing Ds; to adopt high level. Hence, Vz, decreases and 7z to increase, while the 
energy is transferred from the harvester into the inductor. At the same time, the high 
level of Ds» follows Vz+, keeping S» turned off. Shortly before Vz reaches zero, 
the zero crossing event is detected, hence phase 2 starts by shifting Ds; and Ds» low 
and Ds3 high. The subsequent transportation of the energy from the inductor into 
the buffer capacitor is indicated by Zz which is decreasing. The transfer process is 
terminated by Ds» going high and Ds3 going low. 

In order to compare the performance of the PSCE chip to the performance of a 
simple diode rectifier, the output power has been measured, as shown in Fig. 6.21. 
The optimum output power of 40 uW is achieved at Vbuf = 1.5 V. Referring to the 
optimal output power achieved by the diode rectifier, this means an improvement of 
30 96. At Vyyt = 2.1 V, the PSCE chip extracts 38 yW, whereas the diode rectifier 
harvests 22.5 iW. This corresponds to an improvement of 69 96, and even more can 
be realized for higher buffer voltages. The decreasing curve indicates the increasing 
influence of switching losses. 

Figure 6.22a shows a photograph of the demonstrator board, incorporating the 
PSCE chip and the full bridge rectifier composed of BAT 46 diodes. The switch can 
be used to toggle manually between both interface circuits. The load is composed 
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Fig. 6.22 Photograph of the demonstrator board a showing its components and b mounted to the 
Eichhorn harvester 
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Fig. 6.23 Schematics of the demonstrator board 


of a voltage detector and an light-emitting diode (LED). If mounted to the Eichhorn 
harvester as shown in Fig. 6.22b, the LED flashes with a frequency corresponding to 
the power generated by the harvester. Hence, the advantageous effect of the PSCE 
chip is visualized. 

Figure 6.23 illustrates the schematics of the demonstrator board. As stated above, 
the single pole, double throw (SPDT) switch is used to toggle between the Schottky 
diode rectifier and the PSCE chip. A 100 F capacitor buffers the energy extracted 
from the piezoelectric harvester. If Vpuf exceeds a fixed upper threshold of 2.16 V, 
the voltage detector closes a switch from the buffer capacitor to the LED. Hence, 
the LED flashes and removes energy from the buffer capacitor, causing Vpy to 
decay. If Vbuf falls below the fixed lower threshold of 2.09 V, the switch inside 
the voltage detector opens, and so Vyyr can increase again. Thus, Vbuf is showing 
a sawtooth characteristics, as depicted in Fig.6.24. The flashing frequency is an 
indicator for the power generated by the harvester. For the Schottky diode rectifier, 
the flashing frequency is 1.4 Hz, and for the PSCE chip, the flashing frequency is 
2.3 Hz, representing a 64 % improvement. This is almost identical to the improvement 
of 69 96 observed at Vyyr = 2.1 V, as depicted in Fig. 6.21. 
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Fig. 6.24 Measured buffer voltage with connected Schottky diode rectifier and PSCE chip, showing 
the characteristic sawtooth shape (L = 2.2 mH) 


6.4 PSCE Chip Performance 


In this section, the performance of the PSCE chip is evaluated by means of mea- 
surement and simulation. All the measurements shown in this section have been 
performed using the Mide V22B harvester as a power source which has been char- 
acterized in Sect. 6.2.3. 


6.4.1 Measurement Setup and Definitions 


Figure 6.25 shows the PCB used for measurements with the PSCE chip, which is 
attached in a CLCCAS socket in order to establish the connection between the package 
and the outer world. The two-wire piezoelectric harvester output is connected via the 
corresponding connector pins. Arbitrary (electrolyte) capacitors and load resistances 
can be plugged into the corresponding connectors. Four different Murata inductors 
are mounted which can be selected via a jumper. Due to their larger size, one single 
Panasonic inductor can be plugged into the corresponding connector. A multi-pin 
connector is used to monitor the digital outputs. A jumper can be used to toggle 
between SECEI3 and SECE1323 switching techniques. In combination with small 
shunt resistors, several INA110 instrumentation amplifiers with selectable gain factor 
allow for measuring the inductor current, the buffer current, and the piezoelectric 
current. For testing purpose, external voltages can be applied by several banana 
plugs, e.g. in order to supply the switch drivers externally. 
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Fig. 6.25 PCB utilized for measurements with the PSCE chip 


Table 6.4 Definition of the excitation levels with corresponding measured characteristic voltage 
and power values 


Excitation level f Vp „OC Vp ,PSCE Pres,max PPSCE,in min 
(Hz) (V) (V) (uW) (uW) (m/s? 


Open circuit resonance 10.0 


The acceleration of the shaker has been chosen such that the switching technique 
SECE23 can be enabled over the whole tolerated buffer voltage range 1.4V < 
Vout < 5 V. Assuming that SECE23 is active for 


Vado e. (6.8) 


according to (5.6), the piezoelectric voltage amplitude must be higher than 2.2.5 V — 
11 V. According to Fig.6.8, a shaker acceleration amplitude of @ = 10 m/s? 
sufficient to meet this criterion, assuming that the harvester is driven at resonance. 
As the SECEI3 technique can be activated for any Vp > 0, as given in (4.26b), this 
technique has not to be considered for the proper choice of acceleration. 
Maintaining the acceleration amplitude of @ = 10 m/s?, smaller piezoelectric 
output voltages can be achieved by driving the harvester out of resonance. It is 
important to keep the acceleration constant since the model parameters of the piezo- 
electric harvester as summarized in Table6.1 are only valid for this acceleration, 
due to the nonlinear behavior of the harvester. For characterizing the PSCE chip, 
different excitation levels are defined as listed in Table 6.4 in order to have equal 


152 6 Performance Analysis of the PSCE Chip 


comparable conditions. ! The lower limit is only applied for the power loss analysis 
at Vour = 2.5 V, representing the only buffer voltage value for which the efficiency 
is positive (see Sect. 6.4.3.2). If the excitation frequency is increased to 166.7 Hz, the 
power provided by the piezoelectric harvester is sufficient to overcome the chip losses 
over the whole tolerated buffer voltage range, even for the inductors resulting in the 
worst performance. The corresponding excitation level is called far off-resonance. 
Since the PSCE voltage amplitude is as low as 2.9 V, there exists no buffer voltage 
within the tolerated voltage range for which (6.8) is true. Thus, SECEI3 is always 
active. The excitation level called off-resonance is defined such that both switching 
techniques SECEI3 and SECE23 are active within the tolerated buffer voltage range. 
The excitation level called PSCE resonance defines the state with the highest possible 
output power, leading to the fact that SECE23 is enabled over the whole tolerated 
buffer voltage range. 

For each excitation level, the following measurement values are given: The vibra- 
tion frequency f required to achieve the corresponding excitation level, the open 
circuit voltage amplitude VP. oc; the voltage amplitude with connected PSCE chip, 
Vp. PscE. the optimum RMS value of the piezoelectric output power dissipated in 
a resistive load, Pres max, and the piezoelectric power flowing into the PSCE chip, 
Pesce, in- In addition, Gmin gives the minimum shaker acceleration which is needed 
to achieve the indicated harvester voltage and power values, assuming respective 
operation at resonance. 2 The values of VP oc and Vp PscE are provided in order to 
observe the electromechanical feedback influence on the improvement achievable 
by the PSCE circuit (see Sect. 6.4.4.2). Whereas the value of Pres max has been mea- 
sured using load sweeps, the values of Ppscg,in are calculated based on Vp.PscE. 
according to 


PPSCE,in = T = fCp VP PSCE> (6.9) 


where T = 1/f denotes the period of the external excitation, and Epp = 0.5Cp 
Y cst indicates the energy stored on the piezoelectric capacitance Cp which is 
removed twice per period during a transfer process. 

In order to determine the harvested power using the PSCE chip, the setup shown 
in Fig.6.17b is used. Table6.5 summarizes the measured DC resistances of the 
employed inductors. The volumes of the Panasonic and the Murata inductors are 
0.63 and 0.35 cm? , respectively. One has to connect a sufficiently large buffer capac- 
itor in the microfarad range in order to achieve a small ripple in the buffer voltage. 


! As discussed in Sect. 6.2.4.1, due to the different shapes of the frequency curves as a result from 
the nonlinearity of the piezoelectric harvester, it is evident that different frequencies are required in 
measurement and simulation to achieve the respective excitation levels. For the off-resonant and far 
off-resonant excitation level, the frequency has been set such that the open circuit voltage amplitude 
matches between simulation and measurement. 

? The minimum shaker acceleration values are indicated only for information, in order to give a 
feeling of the values necessary in real applications—as stated before, the different excitation levels 
are achieved varying the excitation frequency, maintaining a constant vibration amplitude. 
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Table 6.5 Measured 


E Murata 2200R 
parasitic DC resistance of the 
inductors (given in Ohms) 6.4 
used in combination with the 10.9 
PSCE chip 397 





In combination with a load resistor Rz connected in parallel, an invariant buffer 
voltage Vbuf,ss will establish after some time which can be measured by a digital 
multimeter. This so-called steady state is achieved when the power dissipated in the 
load resistor equals the extracted power Ppscg, which can be calculated according 
to (6.7): 





PPSCE, out = Phary = (6.10) 


From (6.9) and (6.10), the efficiency of the PSCE chip can be calculated using 


) 
mM PPSCE,out u Vout.ss l 6.11 
chip = pun R, o fC.02 ue 
PSCE,in L | fCpVg 


The definition of the chip efficiency assumes that the piezoelectric capacitance is 
completely discharged during a transfer process. Due to imperfections in the control 
circuitry, this is not the case in reality, so the actual piezoelectric output power is 
smaller and thus the chip efficiency is higher than calculated. 

The chip efficiency takes into account all the losses incorporated in the chip, 
i.e. the conduction losses, the switching losses, and the power consumption of the 
control circuitry. A detailed comparison of these loss components is discussed in 
Sect. 6.4.3. Besides the chip efficiency, it is also interesting to know how much power 
can be extracted in relation to the theoretical maximum power extractable from the 
piezoelectric harvester. This so-called harvesting efficiency is determined using 


P ys 8d 
hary = “PSCE,out ^ _buf,ss , —: (6.12) 
Piim RL (ma) 








where Piim is given in (2.38). Since Pj is larger than PPsCE,out, the harvesting 
efficiency is smaller than the chip efficiency. 

The inductor current and the buffer current are measured using 1 $2 shunt resistors, 
as depicted in Fig. 6.26, in order to not modify the circuit characteristics significantly. 
The voltage drop at the shunt resistors is amplified with a factor of 100 by an INA 110 
instrumentation amplifier which has very highinput impedance inputs (7:107). 
From the single ended voltages representing the amplified voltage drop, the actual 
current can be calculated by 


Vir /100 V 100 
pos HEP Lnd dispu BURMA. 


19 —* io 
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Fig. 6.27 Measured steady state transient characteristics of the piezoelectric voltage and the buffer 
voltage over several periods for operation at resonance (Cyyr = 1 WF, Ry, = 14 kQ) 


6.4.2 Transient Characteristics 


All the measurements shown in this subsection are performed using the 10 mH Pana- 
sonic inductor. 

First of all, the transient characteristics of the PSCE chip are evaluated using the 
setup shown in Fig. 6.17b. Figure 6.27 shows the piezoelectric voltage and the buffer 
voltage in steady state over several cycles. The characteristics of the piezoelectric 
voltage match the theoretically predicted characteristics as illustrated in Fig. 3.7. 
The transfer happens at the instants when the piezoelectric voltage peaks, resulting 
in a rapid decay of the piezoelectric voltage. The rest of the time, the piezoelectric 
harvester is driven at open circuit resonance. For demonstration purposes, a relatively 
small buffer capacitance has been chosen for the plot in Fig. 6.27, hence a significant 
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Fig. 6.28 Measured steady state transient characteristics of the piezoelectric voltage, the SECE 
control signal and the NVC control signal over several periods for operation at different excitation 
levels (Vbuf = 2.5 V). a Far off-resonance. b Resonance 


ripple of the buffer voltage can be observed. The sawtooth shape of the buffer voltage 
indicates that the buffer capacitor is charged during a very small time interval only, 
and during the rest of the time, the load resistor discharges the buffer capacitor. The 
small buffer capacitance is chosen only to visualize the proper operating principle 
here. In reality, a small ripple is desired, so a much larger buffer capacitance would 
be chosen. From Fig. 6.27, the current harvested power can be calculated according 
to (6.10). The mean value of the buffer voltage is Vbuf,ss = 2.45 V, and using a load 
resistance Rz = 14 kQ, the harvested power is calculated as 


2 
Vout „SS 


L 





Fhar = = 429 uW. 


Figure 6.28 shows the piezoelectric voltage, the NVC voltage, and the digital con- 
trol signals of the NVC and the SECE selector, for operation at far off-resonance 
and at resonance. According to Sect. 5.4, there are two NVC control signals DNvca 
and Dwvcp. Since Dwvcp equals Dynvyca, only one signal is shown in the plots 
here. From Fig.6.28a, it becomes evident that the peak is detected slightly before 
the actual maximum of the piezoelectric voltage, i.e. slightly before the piezoelectric 
beam has reached the point of maximum deflection (see Fig. 3.7). This means that the 
beam moves towards the maximum deflection after the transfer process has occured, 
slightly increasing the voltage in the same direction. This is the reason for the small 
bumps appearing shortly after the transfer process. In contrast, for the resonant exci- 
tation shown in Fig.6.28b, the peak is detected exactly at the point of maximum 
beam deflection, leading to a smooth increase of the NVC voltage after the transfer 
process. The reason for the premature peak detection for low voltage amplitudes 
has been discussed in Sect. 5.6.2. An example shows the loss incorporated due to 
a premature peak detection. Assuming that the peak is detected at 3 V instead of 
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3.2 V, a power of 166.7 Hz - 19.5 nF - (3 Vv)? = 29.34W is harvested instead of 
166.7 Hz - 19.5nF - (3.2 V)? = 33.3 uW according to (6.9), representing a loss of 
12 %. Hence, since the output power is a function of the squared voltage amplitude, 
special care has to be taken for the design of the peak detector. 

The control signal of the SECE selector, Dsgce, indicates which switching tech- 
nique is currently active (see Sect. 5.2). According to (6.8), the piezoelectric ampli- 
tude must exceed 2.2 Vpuf for SECE23 to be enabled. Here, the buffer voltage is set to 
2.5 V, so the critical point is 2.2 - 2.5 V — 5.5 V. For the far off-resonant excitation, 
Dsgcg is permanently high, meaning that SECEI3 is active due to the small VL+ 
amplitude. In contrast, for the excitation at resonance, Dsgcg goes low when VL+ 
rises above 5.5 V, indicating that SECE23 is active. 

Figure 6.29 illustrates significant analog and digital control signals during the 
transfer process for different representative excitation levels. For all the digital control 
signals except Ds», the low level is GND and the high level is Vbuf, so these signals 
are represented as real digital signals being high or low. As described in Sect. 5.1.3, 
the high level of the control signal Ds» is given by max (Vyyr, Vr, .), meaning that 
the high level equals V; if Vr- > Vbur holds. Thus, Ds» is displayed as an analog 
signal. The inductor and the buffer current are measured as described in Sect. 6.4.1. 

As discussed in Chap. 4, the NVC voltage VL+ = |Vp| and the inductor current 
Ij, both define the characteristics of the switching techniques. The transfer process 
is always initiated by the peak detector which provides a short pulse on Dpp when 
Vr, has reached a maximum. The characteristics of SECEI3 can be observed by 
means of Fig. 6.29a and b (see also Sect. 4.4.2 for theoretical treatment of SECE13). 
During the first phase of the transfer process, the switch MOSFET Ms; is on, and 
the ZCD is activated in order to be ready to determine the instant when Vz+ crosses 
zero. Since all the current provided by the harvester is flowing towards GND, the 
buffer current is zero. During the second phase which is activated when Vz has 
reached zero and /z has reached a maximum, indicated by a short pulse of Dzcp, 
Ms» and Ms; are triggered. So, Iz and wf are (ideally) equal. In addition, the RCD 
is enabled, as described in Sect. 5.8.2, in order to ready to detect the instant when 7z 
has reached zero. 

The main difference of SECE23 compared to SECE13 is the fact that 77 is flowing 
directly into the buffer instead to GND during the first phase of the transfer process 
(see also Sect. 4.4.3 for theoretical treatment of SECE23). This means that switch 
Ms; is not needed at all, so Ds; is always low, as can be seen in Fig. 6.29c. Hence, the 
buffer current Jpuf is equal to 77, for the whole transfer process. Another consequence 
is the fact that Iz peaks before VL has reached zero. The direct comparison between 
Fig. 6.29b and c shows that /7 is lower compared to SECE13, leading to a lower power 
loss in the parasitic series resistances, as discussed theoretically in Sect. 4.5.3.2. 

According to the theoretical discussion presented in Sects. 4.4.2 and 4.4.3, the 
duration of the first and the second phase of the transfer process is given by 


JU 
Ti = t3 —to = y LCp- 3 and 75 = fgi4 — t3 =V LCp x (6.13) 
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Fig. 6.29 Measured transient characteristics of significant analog and digital signals during the 
transfer process for operation at different excitation levels (Vbuf = 2.5 V). a Far off-resonance. 
b Resonance (SECE13). c Resonance (SECE23) 
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Table 6.6 Comparison of the duration of the first and second phase of the transfer process for 
calculation and measurement 


Ti (ws) Tə (us) 
Calculation Measurement Calculation Measurement 


for SECE13 and 


—] | 2 
Ti = to3 — tg = y LC p - arccos (=) and 75 = fgo23 — h3 = / LCp :x,/| l — — 
X X 


(6.14) 


for SECE23, where x = Vp/ Vout. For the given case, the parameters in (6.13) 
and (6.14) are L = 10 mH, Cp = 19.5nF, Vy; = 2.5 V, Vp = 2.9V for the 
excitation at far off-resonance, and Vp = 12.8 V for excitation at resonance. The 
corresponding durations observed in measurement are extracted from the plots as 
indicated in Fig.6.29. Table6.6 summarizes the calculated and measured values. 
Due to the power loss which is not considered in theory, the measured time intervals 
are always shorter than the calculated ones. For T1, the values match better between 
simulation and calculation than for 75. The reason for this behavior is the offset 
incorporated in the RCD in order to avoid reverse currents, as described in Sect. 5.6.4, 
causing a premature end of the transfer process. 


6.4.3 Losses 


This section discusses the different loss mechanisms incorporated in the PSCE chip, 
reducing the harvested power. At first, the total power loss will be investigated in 
Sect. 6.4.3.1 by means of simulations and measurements. In order to understand the 
loss mechanisms which make up the total power loss, a detailed breakdown of the 
loss components is given in Sect. 6.4.3.2. This analysis is performed solely based on 
transistor-level simulations using the coupled piezoelectric generator model shown in 
Fig. 2.10b, considering electromechanical feedback. Hence, a realistic and consistent 
insight into at least the relative relations between the loss components is possible. 
In order to capture each loss component in measurement, i.e. to determine absolute 
losses, ultra-wideband and very precise shunt amplifiers would have been necessary 
due to the occurence of fast switching transitions during energy transfer. To the 
author's opinion, the benefit of purchasing such equipment is out of all proportion 
to the costs. 
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Fig. 6.30 Comparison of chip input and output power at far off-resonance. a Chip input and output 
power. b Power loss [difference of the curves from (a)] 


All the results shown in this section are gained using the 10 mH Panasonic inductor 
with a series resistance of 4.3 Q. 


6.4.3.1 Total Power Loss 


In order to get a first insight into the power losses, the power provided by the piezo- 
electric harvester which represents the chip input power, and the power being actually 
transferred into the buffer which represents the chip output power, are compared with 
respect to each other. The difference between both indicates the power consumed by 
and lost in the chip. The input power for measurement and simulation is determined 
according to (6.9) which makes it necessary to know the excitation frequency, the 
electric piezoelectric harvester capacitance and the piezoelectric voltage amplitude 
with connected PSCE chip. Since (6.9) assumes that the energy stored in the piezo- 
electric capacitance is extracted completely during a transfer process, this equation 
provides an upper estimate for the input power. In order to guarantee a robust trans- 
fer process, the zero crossing detector terminates the energy extraction out of the 
harvester at a voltage which is slightly larger than zero (see Sect. 5.6.3), so a little 
amount of energy could remain in the piezoelectric capacitance. The output power in 
measurement and simulation is determined by means of (6.10). In the following, the 
comparison of the input and the output power is performed at a very low and a very 
high input power, i.e. for excitation at far off-resonance and at resonance, according 


to Table 6.4. 
Far Off-Resonance 


Figure 6.30a compares the measured and simulated input and output power for exci- 
tation at far off-resonance. As expected, the measured input power is quasi-constant 
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over the buffer voltage, meaning that the load does not affect the chip input power. 
Due to the chip losses, the output power is below the input power, and has a maximum 
at Vyyr = 2.5 V. The simulated input power curve shows a slight dependency on the 
buffer voltage. Moreover, the simulated input and output power curves exhibit higher 
values compared to the measured ones. These two observations indicate that the sim- 
ulation model does not perfectly fit to the real harvester device in the off-resonant 
regime. The difference between the input and output power curves representing the 
chip losses is plotted in Fig. 6.30b. The measurement and simulation curves show a 
very good matching, at least for buffer voltages below 3.5 V. For higher buffer volt- 
ages, the measured losses exceed the simulation ones, indicating that some losses are 
not included in the simulation. The power loss curves suggest that there is an optimal 
operating point at Vbuf = 2.5 V where the losses are minimum, corresponding to the 
point of maximum output power as noticed before. A detailed insight into the power 
loss mechanisms is given in Sect. 6.4.3.2. 


Resonance 


At resonance, there are two possible switching techniques: SECE23 and SECEI3. 
Thus, the input and output power comparison is performed using those two techniques 
for measurement and simulation in the following. 

Figure 6.31a and b shows the measured and simulated input and output power 
for SECE23 and SECE13. Due to the excitation at resonance, the power values are 
much higher compared to the excitation at far off-resonance. The input power is 
quasi-constant with regard to the buffer voltage, and is the same for SECE23 and 
SECEI3. The output power is lowest for Vbuf = 1.5 V and increases for increasing 
buffer voltage, indicating that the chip losses decrease. Comparing the graphs for 
SECE23 and SECE13, it is evident that SECE23 allows to harvest more power. 
The power loss curves shown in Fig. 6.31c representing the difference between the 
input and output power directly illustrate the advantageous influence of SECE23. 
Using SECE23, the highest power loss reduction of 35 % is possible at Vbuf = 5 V 
in measurement. Generally, the power loss niveau is quite high for the excitation 
at resonance, especially for low buffer voltages. As discussed in Sect. 6.4.3.2, the 
conduction losses play a major role. 

Similar to the excitation at off-resonance, the simulated input and output power 
curves are higher than the measured ones. In constrast, their shape is the same, in the 
sense that lowest and highest power is harvested at the borders of the buffer voltage 
range. The power loss comparison confirms that SECE23 reduces the power losses, 
but the achievable improvement is lower compared to the measurement. This indi- 
cates that the real PSCE chip incorporates higher losses compared to the simulation 
due to the presence of vias, wires, etc., resulting in more power which can be saved 
by using SECE23. 
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Fig. 6.34 Comparison of the chip input power and the output power at resonance. a SECE23. 
b SECEI3. c Power loss comparison for SECE23/SECEI3 


6.4.3.2 Breakdown of the Loss Components 


In the following, the total losses discussed in Sect. 6.4.3.1 are broken down, clari- 
fying which loss components contribute how much to the total power loss. The loss 
components are discussed as a function of the buffer voltage for the excitation at far 
off-resonance and resonance. Finally, a comparison is performed at a fixed buffer 
voltage, including a loss breakdown at the lower limit of acceptable input power for 
which the PSCE chip is operating. The whole analysis is based on simulations. 

One of the main loss components is the conduction loss (Peona) due to parasitic 
series resistances of the negative voltage converter ( PNyc), the MOSFET switches 
(Ps1, Ps2, Ps3), and the parasitic series resistance of the inductor ( Pac) and the buffer 
capacitor (not considered here). The position of these resistances is shown in Fig. 4.3. 
The total conduction losses can be calculated as 


Peona = PNvc + Psi + Ps2 + Ps3 + Pac. (6.15) 


According to (4.9), in general, the average power loss dissipated in a series resistance 
Rs can be calculated as 
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tE 
Pw = 58s | i dt, (6.16) 
10 
where fo and tg denote the start and the end of the current flow indicated by Iz. As 
discussed in Sect. 4.4, the integral over the squared current which is substituted by 
the symbol p can be evaluated for the two phases of the transfer process, separated 


into SECEI3 and SECE23. The corresponding expressions are displayed again for 
convenience: 


Pits = Epo 


P313 = Epo 





P2,23 = Epo 


P3,23 = Epo 
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3 x 
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where x = V and Epo = +Cp Ve. Note that the p coefficients are used for 
explanations only; the power losses shown in the plots are determined directly from 
simulations. 

According to (5.3), the resistance of the conducting channel of a MOSFET can 
be calculated as 


“19 ae » JJ? 


1 
Po o 6.18 
" æ (Vas — Vr) ) 


where oa is a factor incorporating process parameters and the MOSFET dimensions, 
Vgs is the gate-source voltage, and V7 the threshold voltage. Combining (6.16) and 
(6.18), the power dissipated in the passive part and the active part of the NVC can 
be expressed as 


[pe l E E . j) pi for SECE13 


Vit—Vrp)  @(VL+—Vr,N 
PNVC,pass = (6. 19) 
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PNVC.act = "reme Vr Nb) + a Vrn) | 500 
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Note that the current extracted from the piezoelectric harvester is flowing through 
two NVC MOSFETs, so there are two terms for the power loss (see Sect. 5.4). The 
threshold voltage of the n-channel and p-channel MOSFETs is denoted as Vr y and 
Vr p. The threshold voltage of the upper n-channel MOSFETs in the active NVC is 
denoted by Vr np, indicating bulk effect. Whereas the transistors of the active NVC 
are actively driven by comparators, the piezoelectric voltage itself turns on and off the 
transistors of the passive NVC. So, the gate-source voltage is time-variant during the 
first phase of the transfer process, and the average value V+ is taken as gate-source 
voltage for the passive NVC MOSFETs. 

The power loss incorporated in the switch MOSFETs can be written as 


1 
-m CENE for SECE13 
Ps; = f o (Vout — Vr,N) P1,13 TOF (6.21) 
0 for SECE23 
1 
An Vow Veo) Pole for SECE13 
puc a (Vout — Vr, P) MEN (6.22) 
o (Vi. —Vr.P) p223 + o (Vout —Vr.p) D323. FOL 
1 
a ee E for SECE13 
Pss = f a(Vout Vra] PS UT (6.23) 
am - Vi) P3905 for SECE23 
Finally, the power loss in the inductor can be expressed as 
1.13 + p3,13). for SECEI3 
Pac = f Rac (p E ) ; (6.24) 
(p2,23 + pa23) for SECE23 


Additionally, there are the parasitic series resistances of the chip wiring, the bond 
wires, the tracks on the PCB etc., which play a role only in the measurements. Thus, 
they are not considered here. 

But there are also other loss components which are important. The active control 
circuitry consumes static (Pact,stat) and dynamic (Pact,dyn) power in order to drive 
the MOSFET switches. Static power is consumed by the components which have to 
be enabled permanently, and dynamic power is dissipated by the components which 
are enabled during the transfer process only: 


P. act,stat — ‘stat Vout (6.25) 


tE 
Pact,dyn = Vout f Tayn dt. (6.26) 


to 


Due to the fact that a supply independent biasing circuitry is used as described in 
Sect. 5.7, the static and the dynamic supply current Istat and Tayn are ideally constant 
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with respect to the buffer voltage, so Pact,stat and Pact,dyn linearly depend on the 
buffer voltage. The dynamic supply current is zero beyond the transfer process and 
rises to values of several ten microamperes during the transfer process. 

Last but not least, there are the cross current losses (Peross) of the digital circuitry, 
mainly caused by the switch MOSFET drivers (see Sect. 5.1.4). Non-inifinitely sharp 
edges in the control signal cause a shoot through current [gross in the CMOS inverters, 
which can be pretty large for high buffer voltages: 


tE 
Poross = Vout f J Toros dt. 
to 


The common switching losses due to charging and discharging of the switch MOS- 
FET gates play an inferior role due to the low switching frequency which is given 
by the external excitation frequency. 

In the following, the loss components are discussed for the excitation at far off- 
resonance and at resonance. Since the conduction losses play a major role for the 
effectiveness of SECE23, these losses are highlighted more detailed than the other 
loss mechanisms. 


Far Off-Resonance 
Figure 6.32 shows the conduction losses, the cross current losses and the power con- 


sumption of the active control circuitry, for excitation at far off-resonance. Note that 
SECE13 is active for the whole buffer voltage range. The conduction losses are domi- 
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Fig. 6.32 Components of the power loss at far off-resonance, determined from simulations: Con- 
duction losses, cross current losses and power consumption of the active control circuitry 
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Fig. 6.33 Components of the conduction loss at far-offresonance, determined from simulations 


nating for small buffer voltages and loosing importance for increasing buffer voltages. 
This corresponds to the characteristics of the energy loss depicted in Fig. 4.16 and 
theoretically discussed in Sect. 4.5.3.2. In the theoretical discussion, constant para- 
sitic series resistances were assumed. In the chip implementation, the control circuitry 
is directly supplied by the buffer voltage, so the buffer voltage determines the gate- 
source voltage of the MOSFET switches. For larger buffer voltages, the gate-source 
voltage is also larger, leading to reduced on-resistance according to (6.18) und thus 
to lower conduction losses. Hence, the influence of the varying channel resistances 
also contribute to the decay of the conduction losses for increasing buffer voltages. 

The active and static power consumption of the control circuitry show an upward 
characteristic for increasing buffer voltages. According to (6.25) and (6.26), the 
curves should exhibit a linear dependency of the buffer voltage. Due to the fact that 
the supply current is not perfectly constant with regard to the buffer voltage, as shown 
in Fig. 5.28, the linearity of the simulated curves is not perfectly given. 

As described before, the cross current losses are negligible at low buffer voltages 
but grow significantly for higher buffer voltages. 

A detailed breakdown of the conduction losses is given in Fig. 6.33. In the follow- 
ing, the characteristic of the NVC is clarified. According to (6.19), one would expect 
the passive NVC to consume a constant power with regard to the buffer voltage, 
because the characteristics of the rectified piezoelectric voltage, V+, and pi, 13 both 
are not a function of the buffer voltage. But as discussed in Sect. 5.4, the active NVC 
takes over more and more current as the buffer voltage increases, reducing the current 
flowing through the passive NVC. Hence, the square integral p1 13 decreases, and so 
does the power loss in the passive NVC. For buffer voltages up to 3.5 V, the power 
loss characteristic of the active NVC corresponds to (6.20), because the on-resistance 
of the NVC MOSFETS decreases and so does the power loss. But for higher buffer 
voltages, the active NVC is conducting more and more current, increasing its square 
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integral p1.13, compensating for the reduction of the on-resistance. Thus, the power 
loss of the active NVC is almost constant for buffer voltages between 3.5 and 5 V. It 
is important to mention that the effect of sharing currents between the passive and the 
active NVC is not considered in the theoretical expression of p1,13 given in (6.17). 

According to (6.21) and (6.23), the power loss curve for the switches S1 and S3 
is decreasing for increasing buffer voltage, because the on-resistance is decreasing. 
Since p.13 is not a function of Vyyr whereas p3.15 is inversely proportional to Vor, 
the power loss curve of S3 is expected to drop faster than S1. This effect is hardly 
visible, but it can be verified by means of the graph. After the power loss curve of S2 
decreases up to Vbuf = 2.5 V, itis increasing while the buffer voltage increases. Since 
(6.22) suggests a similar behavior for So compared to S5, this behavior cannot be 
explained by means of the given equation. A closer look into the transient simulation 
results manifests a timing problem in the reverse current detector, causing a premature 
end of the transfer process. Hence, the remaining current increases the voltage at the 
negative inductor terminal, V; —, so that the channel of S2 can conduct the remaining 
current. Since the timing mismatch increases when the buffer voltage increases, this 
leads to an increasing power loss at the switch S». 

According to (6.24), the power loss within the inductor should decrease for 
increasing buffer voltage, because p3 13 is inversely proportional to Vpyr. But the 
graph suggests a constant characteristic. This can be explained by the fact that the 
input energy per transfer process, E po, is increasing for increasing buffer voltage, 
as shown in Fig. 6.30b. Since the term E po appears in the expression of p3.15, this 
effect causes the inductor power loss to remain constant. 


Resonance 


Figure 6.34 depicts the loss components at resonance for SECE13 and SECE23. From 
Fig. 6.34a and b, it becomes directly evident that the conduction losses dominate the 
other losses by far over the whole buffer voltage range. This is because the conduction 
losses directly depend on the piezoelectric voltage amplitude Vp which is part of the 
p-expressions summarized in (6.17), whereas the other loss components are more 
or less independent on Vp. Similar to the far off-resonant excitation, the conduction 
losses decrease for growing buffer voltages according to the theoretical analysis and 
due to the decreasing on-resistance of the switch MOSFETs. While the other loss 
mechanisms are hardly affected when switching between SECEI3 and SECE23, 
the reduction of the conduction losses due to SECE23 is clearly visible. While the 
conduction losses at Vbuf = 1.5 V can be reduced by 20 %, the reduction factor at 
Vour = 5 V is 100 2o. 

It is worth to have a closer look at the cross current losses and the power consump- 
tion of the control circuitry, as shown in Fig. 6.34c and d. As expected, the selected 
technique, i.e. SECE13 or SECE23, has no significant impact on the control circuitry. 
Comparing the curves with the far off-resonance excitation shown in Fig. 6.32, it is 
noticeable that the static power consumption at resonance is about twice the power 
consumption at far off-resonance. This becomes directly clear having a closer look at 
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Fig. 6.34 Components of the power loss at resonance, determined from simulations: Conduction 
losses, cross current losses and power consumption of the active control circuitry. a SECEI3. 
b SECE23. c SECEI3 (zoomed). d SECE23 (zoomed) 


the SECE selector topology depicted in Fig. 5.7. The resistor divider is connected to 
the output of the NVC, i.e. the rectified output voltage of the harvester, which means 
that the current flowing through the resistor divider is depending on the excitation of 
the harvester. Since the output voltage of the harvester is much higher when driven 
at resonance, the corresponding power loss within the resistor divider goes the same 
direction. 

Figure 6.35 shows a breakdown of the components which make up the conduc- 
tion loss, for SECEI3 and SECE23. Similar to the off-resonant excitation, the power 
loss of the passive NVC is decreasing with increasing buffer voltage since the active 
NVC takes over more and more of the current flowing from the harvester into the 
buffer capacitor. In concordance with the respective equations, all the other loss 
components exhibit a decaying characteristic with respect to the buffer voltage, 
except for the power loss at switch S; which is permanently zero for SECE23. 
The latter observation comes from the fact that S; 1s not used at all in SECE23. It 
is noticeable that at resonance, the power loss in S5 is dominating over the other 
loss mechanisms, whereas this 1s not the case for excitation at far off-resonance. The 


explanation for this observation is given in the next paragraph about the percentaged 
breakdown. 


168 6 Performance Analysis of the PSCE Chip 








(a) 70 (b) 70 
A —~ 60 
a E 50 
2 | *—* PNVC.pass 2 40 
«——4 P T 

cs NVC, E -1 30 
5 | 8 

2 z 20 
z | ô 

am 3 am 10 

1.5 2.0 2.5 3.0 35 40 4.5 5.0 15 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
Buffer Voltage Vbuf (V) Buffer Voltage Vyur (V) 


Fig. 6.35 Components of conduction loss at resonance, determined from simulations. a SECE13. 
b SECE23 


Table 6.7 Breakdown of the chip losses as percentage of total losses, based on simulation using 

the coupled piezoelectric model (Vbuf = 2.5 V) 
Loss mechanism Low power | Far off- Off- 

resonance |resonance | PSCE resonance 


SECE23 | SECE13 
Total absolute loss 3.6 uW 7.6 uW 19.7 uW 64.7uW |92.0 uW 
Control circuitry (dynamic) 1.4 
Control circuitry (static) 349 | 243 [i51 | 80 | 5.6 


Cross currents .0 1.7 1.9 
all 


limit 





Msı channel resistance 


Ms? channel resistance 


Mss; channel resistance 


LN | 
32.4 25.5 
Passive NVC channel resistance 


Active NVC channel resistance |33.8 6.0 


Comparing the curves in Fig. 6.35a and b, it becomes evident that each loss compo- 
nent is reduced by the use of SECE23. Because of that, the total amount of conduction 
loss can be diminished as shown in Fig. 6.34, leading to an improved efficiency. 


Percentaged Breakdown at a Fixed Buffer Voltage 


Table 6.7 summarizes the loss components mentioned before as percentage of total 
losses at a fixed buffer voltage of 2.5 V. A column exists for each excitation level 
listed in Table 6.4. The respective switching technique is automatically determined 
by the SECE selector: SECE13 for the low power limit and the excitation at far off- 
resonance, and SECE23 for the excitations at off-resonance and resonance. In order 
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to have a direct comparison between SECEI3 and SECE23 under equal conditions, 
SECE13 has been additionally forced for the resonance case. The zero power con- 
sumption of S; is always an indicator that SECE23 is enabled. The total absolute 
losses which are used as a reference for the percentaged loss components are given 
in the second row. While the transferred power increases from the low power limit 
to the excitation at resonance, the total losses also increase, mainly caused by the 
conduction losses. 

In the following, the contributions of the main loss components to the total power 
loss are discussed for the two extremal excitation levels. At the low power limit, the 
control circuitry (static and dynamic) contributes with 44.4 %, the conduction losses 
(switches, NVC, inductor) with 45.1 96, and the cross current losses with 10.5 %. At 
resonance (SECE13), 8.0 % of the losses are caused by the control circuitry, 91.2 96 
by the conduction losses, and 1.9 96 by cross currents. Hence, the cross current 
losses always play a minor role. If the total losses are very low, the control circuitry 
contributes the same amount of losses to the total power loss as the conduction 
losses. For excitation at resonance, the conduction losses totally dominate the other 
loss components. 

Highlighting the conduction losses in detail for SECE13, it is evident that at 
low input power, the NVC loss represents a very large fraction compared to the 
total loss (almost 43 96 at the low power limit). For high input power, its impact 
reduces by about one half, ending up at 21.9 %. In contrast, the contribution of the 
switches S» and S5 and the inductor increases to significant fractions (25.5, 15.3 and 
20.7 96), whereas the fraction of S, remains below 8 96. The reason for the different 
characteristics can be found by having a closer look at the p-expressions shown in 


(6.17). From x — we and Epo = 3Cp Vou: we can conclude that p1 13 and p2,23 


are a function of V2 whereas p3,13 and p3 23 are a function of V3. assuming that all 
the other parameters remain constant. Thus, the losses in S; and the NVC which are 
a function of pı 13 and p2 23 show a weaker dependency on Vp than the losses in S2, 
$3 and the inductor which are a function of p3,13 and pa »3. In addition, Sı and S3 
are NMOS transistors, whereas S5 is a PMOS with a higher channel on-resistance. 
So, the cubic dependency on the piezoelectric voltage amplitude has an even greater 
impact on the loss in S5 than on the other loss mechanisms. 

From the total absolute losses for the excitation at resonance, it becomes evident 
that SECE23 can save lots of power compared to SECE13. Comparing the loss mech- 
anisms, the percentage of S;, S3 and the inductor reduces for SECE23, indicating a 
very strong saving potential. Even the conduction loss mechanisms which increase 
their percentages (1.e. 5» and NVC) actually are decreasing their absolute power loss, 
but not as strong as the other ones. 
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6.4.4 Efficiency and Harvested Power 


In this section, the harvested power and the efficiency are discussed as a function of 
the buffer voltage, for different inductors and under varying excitation conditions. 
Efficiency plots as a function of the vibration frequency and the piezoelectric voltage 
amplitude are shown in Sect. 6.4.5. 


6.4.4.1 Efficiency 


Figure 6.36a shows the chip efficiency using the 10 mH Panasonic inductor for dif- 
ferent excitation levels from resonance to far off-resonance. According to (6.11), 
the chip efficiency indicates the ratio of the harvested power to the output power 
provided by the harvester, i.e. the ratio of the chip output power to the chip input 
power. Due to the chip losses, this ratio is always smaller than 100 96. Since the load 
does not affect the piezoelectric harvester characteristic, the input power is constant 
with regard to the buffer voltage, according to (6.9). Thus, the shape of the efficiency 
curves is similar to the shapes of the curves showing the harvested power discussed 
in Sect. 6.4.3.1. The highest efficiency of 85 96 is achieved for the excitation at res- 
onance at Vbuf = 5 V. The off-resonant curve shows a sharp drop at Vbuf = 3.2 V, 
indicating the automatic transition between SECE23 which is active for Vbuf < 3.2 V 
and SECEI3 which is active for Vpy¢ > 3.2 V. The vertical height of the transition 
indicates the improvement factor of SECE23 compared to SECE13. From the plot, 
this height is determined to be 8 %. 

The far off-resonant excitation suffers from the lowest efficiency values due the 
low input power and the SECEI3 technique. It is noticeable that the SECE23 effi- 
ciency of the off-resonant excitation (for Vbuf < 3.2 V) is higher than the efficiency 
of the resonant excitation (which is SECE23 for all buffer voltages). This is because 
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Fig. 6.36 Measured efficiency curves for different excitation levels using the 10 mH Panasonic 
inductor (active switching technique indicated in parentheses). a Chip efficiency. b Harvesting 
efficiency 
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of the conduction losses which are higher for the resonant excitation. Obviously, for 
the off-resonant excitation, the best loss-to-input-power ratio can be observed. 

Besides the chip efficiency, it is also important to put into relation the harvested 
power to the absolute maximum power which can be theoretically extracted. This 
so-called harvesting efficiency is calculated according to (6.12) and presented in 
Fig. 6.36b for different excitation levels. Compared to the chip efficiency curves, 
the resulting harvesting efficiency curves show less efficiency, because the absolute 
maximum extractable power is of course higher than the chip input power. Whereas 
the harvesting efficiency is only about 12 % below the chip efficiency, the curves 
representing the off-resonant and the far off-resonant excitation suggest that 20 % 
or even 2 % of the absolute maximum extractable power is harvested. This is clear 
because the maximum power as defined in (2.38) is independent of the excitation 
frequency. So, in order to come closest to the maximum power, the harvester has to 
be driven at resonance. 


Now, the efficiency improvement between SECE23 and SECEI3 is highlighted, 


as theoretically discussed in Sect. 4.5.3.2. In the optimum, i.e. when Vbuf = us 


holds, the theoretical analysis predicted an efficiency improvement of 5 % using the 
Panasonic inductor with Rac = 4.3 Q, and 10% using the Murata inductor with 
Rac = 23.7 Q, both with L = 10 mH. In the measurement, the optimum buffer 
voltage is 128v — 6.4 V (see Table 6.4). The highest tolerated buffer voltage for 
the PSCE chip is 5 V, so the efficiency improvement is supposed to be close to the 
theoretical value. Figure 6.37 shows the corresponding plots for different inductors 
and inductances. First of all, the general shape of the measured curves correspond 
roughly to the theoretical curves, in the sense that the efficiency improvement is 
increasing for increasing buffer voltage. The optimum improvement of 8 % for the 
Panasonic inductor and 13 % for the Murata inductor is achieved at Vyyg = 5 V, 
as illustrated in Fig.6.37a. Hence, the efficiency improvement is higher than the 
theoretical values. This can be explained by the fact that for the theoretical analysis, 
parasitic series resistances were determined by means of simulations. In reality, these 
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Fig. 6.37 Measured chip efficiency improvement for operation at resonance, for different inductors. 
a L = Il0mH. b L = 4.7 mH 
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resistances are always higher due to the additional wiring inside the chip, the bond 
wires, the tracks on the PCB etc. Hence, the achievable improvement is also larger. 

As can be seen in Fig. 6.37b, a smaller inductance leads to a higherefficiencyim- 
provement. This is because more current is flowing from the harvester to the buffer 
capacitor during a shorter time interval compared to an inductor with a higher induc- 
tance, which causes higher losses. Hence, a greater improvement is possible using 
inductors with smaller inductances. 


6.4.4.2 Power 


This section highlights the performance of the PSCE chip regarding the output power 
and the efficiency using different inductors. For measuring the output power, the setup 
described in Sect. 6.4.1 is used. 

Figure. 6.38 compares the output power using the Panasonic and Murata inductor 
having 10 mH each. As explained in Sect. 6.4.1, the excitation at resonance allows 
the SECE23 technique to be active over the whole buffer voltage range. Thus, a 
comparison of SECE23 and SECE13 over the whole buffer voltage range is possi- 
ble, whereas the SECEI3 can be forced externally as described in Sect. 5.2. All the 
different curves in Fig. 6.38a suggest that with higher buffer voltage, more power can 
be harvested, which is due to the lower conduction losses, as discussed in Sect. 6.4.3. 
As expected, the performance of the Panasonic inductor is higher compared to the 
Murata inductor. This can be explained easily by the wiring resistance which is 
lower for the Panasonic inductor (see Table 6.5), resulting in a smaller conduction 
loss. The graph clearly shows that SECE23 overrides the performance of SECE13 
over the whole buffer voltage range. Interestingly, the smaller-sized Murata inductor 
with SECE23 allows to harvest the same or even more power than the Panasonic 
inductor with the standard SECEI3 technique. Hence, SECE23 helps reducing the 
package size and thus cost of the harvesting system. 
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Fig. 6.38 Measurement of the extracted power at resonance, showing a the absolute extracted 
power b the extracted power improvement (L = 10 mH) 
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Fig. 6.39 Measurement of the extracted power at resonance, showing a the absolute extracted 
power b the relative improvement (L = 4.7 mH) 


Figure 6.38b represents the improvement achieved by SECE23 compared to 
SECE13, given as a percentage of the harvested power attainable with SECE13. 
Whereas the improvement using the Panasonic inductor is around 10%, varying 
only slightly over the buffer voltage range, the improvement achieved by the Murata 
inductor is minimum 14 % at Vbuf = 2 V and rising up to 20 96 at Vous = 5 V. So, if 
the losses are higher due to a larger series resistance, the obtainable improvement is 
also higher. 

A similar analysis has been performed taking the 4.7 mH Panasonic and Murata 
inductor. Due to the smaller inductance, more current is flowing from the harvester 
to the buffer capacitor during a shorter time interval. This causes higher losses, 
reducing the harvested power, as shown in Fig. 6.39a. But at the same time, Fig. 6.39b 
illustrates that the achievable improvement by applying SECE23 is higher. Compared 
to the 10 mH inductance, the improvement can get as high as 25 % using the Murata 
inductor and 15 96 using the Panasonic inductor. 

In order to classify the output power performance of the PSCE chip, it will be 
compared to the output power attained by the NVC rectifier using one BAT 46 Schot- 
tky diode (see Sect. 6.2.4.3). For the curves shown in Fig. 6.40, the PSCE resonance 
frequency has been set for the measurement using the PSCE chip, and the open circuit 
resonance frequency has been set for the measurement using the NVC rectifier (see 
Table 6.4). So, all the curves represent the highest output power achievable with their 
respective interface circuit. The highest power can be harvested using the Panasonic 
inductor with the highest inductance, 10 mH, as shown in Fig. 6.40a. The lower the 
inductance, the higher the losses due to the higher current flow from the harvester to 
the buffer capacitor during a shorter time interval, and so the lower the harvestable 
power. An important performance indicator is the gain which can be achieved com- 
pared to the NVC rectifier. Considering the individual optimal output power values, 
the gain is 80 uW or 20 % in relation to the NVC rectifier output power. Due to the 
fact that the PSCE curves do not decay as much as the rectifier curve for decreasing 
buffer voltage values, the gain is dependent on the operating point. The gain is high- 
est at Vbuf = 1.5 V, there the PSCE chip can extract 187 yW or 97 % more power 
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Fig. 6.40 Measured extracted power using the PSCE chip and the rectifiers for different inductors 
at resonance. a Panasonic inductors. b Murata inductors 
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Fig. 6.41 Measured extracted power using the PSCE chip and the rectifiers for different inductors 
at far off-resonance. a Panasonic inductors. b Murata inductors 


compared to the NVC rectifier. For the Murata inductors, the overall output power 
using the PSCE chip is lower compared to the Panasonic inductors. Considering the 
individual optimal output power values, the gain is 61 yW or 17 % in relation to the 
NVC rectifier output power, and at Vbuf = 1.5 V, the gain is 187 uW or 76 96. 

In contrast to the excitation at resonance, a more interesting situation arises for 
the excitation at far off-resonance, as shown in Fig. 6.41. In the way that the smaller 
inductances lead to lower output power values, the characteristics are similar to 
the excitation at resonance described before. But since the piezoelectric harvester is 
driven far away from its resonance, the overall power level has decreased significantly. 
On one hand, this results in a more dominant influence of the power losses, causing 
the respective curves to show a maximum at Vbuf = 2.5 V and a strong decay 
for buffer voltages above 3 V due to the cross current losses in the switch drivers, as 
depicted in Fig. 6.32. But on the other hand, the NVC rectifier is able to harvest power 
only in a range 0 < Vp < 1.4 V. Remember that for the excitation at resonance, 
the NVC rectifier is able to harvest power over the whole range 0 € Vpurp < SV. 
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So, for excitation at far off-resonance, the PSCE chip is able to harvest power in a 
range where the NVC rectifier would harvest nothing. Considering the individual 
optimum power values, the PSCE chip using the 10mH Panasonic inductor can 
harvest 10 uW or 127 % more power compared to the NVC rectifier, as shown in 
Fig.6.41a. For the 10 mH Murata inductor, the respective gain numbers are 8.5 uy W 
or 108 96. 

Comparing the excitation at resonance and far off-resonance, it is evident that 
the power gain considering the individual optimum power is greatly increased from 
20 to 127 %. Another indicator for this characteristic is the ratio of the open circuit 
voltage amplitude Jose and the amplitude with connected PSCE chip, Vp psc. As 
depicted in Table 6.4, there is an almost doubling for the excitation at the low power 
limit, but for the excitation at resonance, there is only a very slight improvement. 
This verifies the theoretical discussion in Sect. 3.4, predicting a rising gain factor 
for an excitation frequency which moves away from the resonance frequency. The 
theoretical maximum is 400 96, but the losses of the PSCE chip reduce the practically 
achievable gain factor. 


6.4.5 Operation Limits 


In real applications, the excitation level available for energy harvesting may vary in 
terms of vibration amplitude and/or frequency. Hence, it is essential for the inter- 
face circuit to be operable over a wide range of voltage amplitudes, frequencies and 
inductances. In addition, the buffer voltage may not be fixed to a certain voltage, so a 
flexibility with regard to these parameters 1s advantageous as well. In the following, 
the operability of the PSCE chip with regard to these parameters will be discussed, 
especially highlighting the limits. Whereas the analysis of the piezoelectric voltage 
amplitude, the frequency, and the buffer voltage is performed by means of measure- 
ments, the influence of the inductance 1s discussed based on simulations. For the 
sake of comparability, all the measurements have been made using the 10 mH Pana- 
sonic inductor, and the efficiency values have been determined by means of (6.11). 
Generally, whenever one parameter is discussed, all the other parameters are kept 
constant. 


6.4.5.1 Frequency 


In order to discuss the chip operation with regard to the excitation frequency, the 
piezoeletric harvester is replaced by a sine wave generator modelling the vibrating 
beam. Therefore, the thevenin equivalent of the uncoupled piezoelectric model shown 
in Fig.2.10 is used, where a 33 nF capacitor models the piezoelectric capacitance. 
Since the PSCE chip needs a differential input, the 180? phase shifted signals of two 
signal generators with an amplitude of 2.25 V is fed to the chip. The resulting open 
circuit amplitude is 4.5 V, so that a PSCE amplitude of 9 V can be observed due 
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Fig. 6.43 Measured transient characteristics of the NVC output voltage for different frequencies 
(Vout = 3 V) 


to the missing electromechanical feedback. Figure 6.42 depicts the corresponding 
measurement setup. This way, it is easy to sweep through a broad band of frequency. 
Using real harvesters, the resonance frequency of the harvester would have to be 
readjusted for each excitation frequency, causing a huge effort without gaining more 
information. 

Figure 6.43 shows the NVC output voltage for different frequencies at a fixed 
buffer voltage. It is directly noticeable that the voltage amplitudes differ between 
consecutive peaks. This is probably due to an inaccurate phase shift which slightly 
deviates from 180°. For the characteristics at different frequencies, this misbalance 


6.4 PSCE Chip Performance 177 


Chip Efficiency 





10? 10? 
Frequency (Hz) 


Fig. 6.44 Measured chip efficiency as a function of the frequency (Vp =9V, Vinee = 3V) 


is irrelevant. The first important observation 1s the fact that the chip is automatically 
adapting its timing for excitation frequencies down to 50 Hz and up to 2,000 Hz, 
without external interference. However, the value of the series capacitance used in 
the peak detector described in Sect.5.6.2 has been optimized for a frequency in the 
range of several hundred hertz. The transient signal for 150 Hz suggests that at this 
frequency, the peak is detected exactly at its maximum voltage. But 500 Hz also 
results in acceptable peak detection events, although the peak is detected slightly 
late. 

For frequencies lower than the optimum frequency 150 Hz, the peak is detected 
prematurely, and for higher frequencies, the peak is detected late. For frequencies 
lower than the lower limit of 50 Hz, the current flowing through the series capacitance 
will not exceed the fixed bias current, thus no peak can be detected. For frequencies 
above 2,000 Hz, the peak detection will still work, but at some point, the excitation 
period comes into the range of the transfer process duration, and in this case no benefit 
can be achieved compared to a simple rectifier. For frequencies beyond the given 
tolerated frequency range, it may also happen that the transfer process is initiated at 
a piezoelectric voltage as low so that transferred energy is below the chip losses. In 
such cases, the chip will stop operation as well. 

The chip efficiency curve given in Fig. 6.44 shows a plateau for frequencies in a 
range 120...600 Hz for which the highest efficiency is achieved. Due to the non- 
accurately timed peak detection, power is lost at frequencies beyond this range, 
resulting in reduced efficiency. Note that the calculated chip input power may not 
be precise because the piezoelectric peak voltages are not constant over time, as 
described before. So, the absolute efficiency numbers may be inaccurate, but the 
general characteristics should have acceptable quality. 
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Fig. 6.45 Measured transient characteristics of the piezoelectric voltage for different piezoelectric 
voltage amplitudes (Vbuf = 2.5 V) 


6.4.5.2 Piezoelectric Voltage Amplitude 


The piezoelectric voltage amplitude may easily change when the vibration conditions 
change with respect to vibration amplitude or frequency. In addition, piezoelectric 
material properties may vary over time, leading to different output voltage levels for 
the same acceleration. So, it is important to maintain proper timing over a wide range 
of voltage amplitudes. 

The maximum tolerated voltage amplitude is 20 V due to the used transistors. 
Once the PSCE chip has performed startup properly, the piezoelectric voltage ampli- 
tude can be decreased to 1.3 V in order to achieve a stable Vbuf = 2.5 V with a 
load resistance of 10 MQ. According to Table 6.4, the corresponding input power 
is 5.7 pW, so this state is defining the excitation level called low power limit. For 
even lower piezoelectric voltage amplitudes, the chip losses override the power pro- 
vided by the harvester, causing the chip to stop operation. If the buffer voltage is set 
externally by a separate voltage supply, the voltage amplitude can be further reduced 
to 1.1 V, showing the absolute tolerated minimum. For amplitudes lower than that, 
the peak detector is not operable because the current flowing through the series capac- 
itance does not exceed the biasing current. For the used harvester, this value has no 
practical meaning, because the PSCE chip is not able to boost the buffer voltage to 
a minimum required value for proper startup. 
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Fig. 6.46 Measured chip efficiency as a function of the piezoelectric voltage amplitude (Vbuf = 
2.5V) 


Figure 6.45 compares the piezoelectric voltage characteristics for amplitudes from 
1.3 to 20 V. Except for the 1.3 V amplitude, the peak instants are detected very 
accurately at the respective maximum voltage. The circuit is still operable for Vp = 
1.3 V, but the peaks are detected slightly prematurely, leading to a small power drop. 

The chip efficiency as a function of the piezoelectric voltage amplitude is provided 
in Fig. 6.46. At the low power limit for Vp = 1.3V, the chip losses are only slightly 
larger than the input power, leading to a bad efficiency. The optimum efficiency 
of 70 96 is achieved at Vp = 7.5 V, and while increasing the piezoelectric voltage 
amplitude further towards 20 V, the efficiency decreases only slightly by 5 %. This 
indicates that the chip losses and the input power increase in a similar way as the 
voltage amplitude increases. 


6.4.5.3 Buffer Voltage 


The SECE principle predicts that the load does not affect the input characteristics, 
as discussed in Sect. 3.3. Hence, a varying buffer voltage should not alter the piezo- 
electric voltage amplitude. For the excitation at resonance, transient plots showing 
the piezoelectric voltage at three different buffer voltage levels are given in Fig. 6.47. 
Visually, there is no difference in the amplitude. The graphs showing the input power 
for resonance and far off-resonance in Sect. 6.4.3.1 exhibit a marginal dependency 
on the buffer voltage. This is at least true for the measurements —-the simulation 
shows a deviation from that behavior, probably caused by the piezoelectric model. 
Since the input power is calculated by means of the piezoelectric voltage amplitude, 
a constant input power is an indicator for a constant piezoelectric voltage amplitude. 

Figure 6.48 depicts the efficiency curve already shown in Fig.6.36a. Ideally, i.e. 
without any chip losses, the curve should be congruent with the input power curve. 
Any deviation from that are caused by the chip losses, as discussed in detail in 
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Fig. 6.47 Measured transient characteristics of the piezoelectric voltage at resonance for different 
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Sect. 6.4.3. As written in Sect. 6.4.6, the minimum trigger voltage is 1.4 V, so this 
voltage defines the low buffer voltage limit. The upper limit is given by the absolute 
maximum voltage being tolerated by the 3.3 V transistors used in the switch control 
circuitry, which is 5 V (see second paragraph in Sect. 5.8). 


6.4.5.4 Inductance of the External Inductor 


For the analysis of the operation limits concerning the inductance, the parasitic series 
resistance should not affect the operation of the PSCE chip too much. The parasitic 
series resistance of real inductors strongly depends on its inductance and its size, 
which makes it very difficult to find or design inductors covering a broad range of 
inductances, and at the same time keeping its series resistance constantly low. In 
addition, it is difficult to find inductors with inductances higher than 100 mH at all. 
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Fig. 6.49 Simulated transient characteristics of the piezoelectric voltage at resonance for different 
inductances (Vgyg = 2.5 V) 


Therefore, the analysis of the operation limits concerning the inductance is performed 
by means of simulations instead of measurements, setting the series resistance to 
4.3 Q. 

Figure 6.49 shows the transient piezoelectric voltage amplitude for different induc- 
tances, whereas the buffer voltage is set to 2.5 V. As expected from the theoretical 
analysis, the duration of the transfer process is affected by the varying inductances. 
For the 100 pH inductance, the oscillation frequency of the tank circuit composed of 
the piezoelectric capacitance C p and the inductance L is very high. Thus, the energy 
transfer process is much shorter compared to the vibration period, indicated by the 
piezoelectric voltage which looks like a vertical line during the transfer process. This 
means that the energy is transferred instantaneously to the load. In contrast, for the 
] H inductance, the decay of the piezoelectric voltage is clearly visible, meaning that 
the transfer process takes much longer due to the low oscillation frequency of the 
Cp — L tank circuit. The instant of the peak detection itself is not affected by the 
varying inductance. 

Figure 6.50 illustrates the chip efficiency as a function of the inductance. As 
theoretically discussed in Sect. 4.5.3.2, the efficiency is increasing for increasing 
inductance, because the current flowing through the chip during the transfer process 
is decreasing, causing less losses in the parasitic series resistances. But on the other 
hand, the transfer process takes longer. Hence, since the power consumption of the 
PSCE chip during the transfer process is quite high compared to the off period, the 
chip consumes more average power for increasing inductance, reducing the efficiency 
for high inductances. For the given piezoelectric harvester, the optimum efficiency is 
achieved for the 100 mH inductance. Note that in reality, due to the series resistance, 
the optimum is supposed to be at lower inductances. 
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Fig. 6.50 Simulated chip efficiency as a function of the inductance at resonance (Vbuf = 2.5 V) 


For inductances below 100 WH, the piezoelectric voltage is falling too fast, causing 
the zero crossing detector (Fig.5.20) to fail. When the inductances are increasing 
above 1H, the piezoelectric voltage shows oscillations after the transfer process, 
because the oscillation cancellation cannot remove the residual energy on the inductor 
fast enough. These oscillations cause the PSCE chip to fail, as described in Sect. 5.3. 


6.4.6 Startup 


The startup functionality has been evaluated by varying the output power of the 
harvester, starting from an uncharged buffer capacitor. Therefore, the excitation fre- 
quency of the shaker has been increased with a slew rate of 100 Hz/min, starting 
from 165 Hz which means far off-resonance. At an excitation frequency of 173 Hz 
which is close to the resonance frequency, the shaker holds the frequency for 2 s, and 
afterwards the frequency is reduced back to 165 Hz. In order to protect the PSCE 
chip from overvoltages, a load resistance of 80 kQ is connected in parallel to the 
buffer capacitor. Its 1 F capacitance is chosen to be quite low, in order to produce 
a significant ripple to demonstrate the buffer voltage characteristics for passive and 
active charging. Figure 6.51a shows the corresponding characteristics of the buffer 
voltage and the switch control signals. After the shaker has been started, the buffer 
voltage is slowly increasing as the buffer capacitor is charged passively via the CMOS 
diode inside the bypass. At the same time, the power supply of the control circuitry is 
decoupled from the buffer capacitor, so the PSCE chip is turned off. This is indicated 
by the switch control signals which do not show any activity. The low ripple of the 
buffer voltage shortly before the chip starts up 1s typical for the passive charging via 
a diode. 

Figure 6.51b illustrates a zoomed view into the transition from the off state to 
the on state. Immediately after startup has occurred at the startup trigger voltage 
which is 1.4 V, the buffer voltage rises much faster, illustrating the advantageous 
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Fig. 6.51 Transient measurement of the startup (RL = 80kQ, Chur = 1 WF), using the follow- 
ing shaker profile: 165 Hz — 173 Hz (2 shold) — 165 Hz. a Complete power-up/-down cycle. 
b Zoomed view into the transition from off to on state. c Zoomed view into the transition from on 
to off state 


effect of the PSCE chip. The sawtooth shape of Vbuf indicates that the PSCE chip 
is operating properly. Since the harvester voltage amplitude starts with a large value 
immediately after startup and afterwards is decaying towards its steady state, and the 
buffer voltage is increasing at the same time, there is a transition from SECE»; to 
SECE;5 which is visualized by means of the control signal for Ms. After the buffer 
voltage has reached and maintained 4.5 V for two seconds, the harvester power is 
decreased by decreasing the shaker frequency back to 165 Hz, resulting in a drop of 
Vout due to the connected Ri. At the shutdown trigger voltage, the PSCE chip stops 
operation and R, draws the buffer capacitor empty. 

As can be observed in Fig.6.51c, the PSCE chip shuts down at 1.0 V which is 
below the startup trigger voltage. This is due to the hysteresis present in the startup 
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trigger, as described in Sect.5.5.1. As Ds» is no longer showing activity after the 
shutdown, the SECE principle is no longer active. The activity of Ds; after the 
shutdown is due to the fact that some components of the digital logic are directly 
connected to the buffer capacitor instead of the control circuitry supply. When the 
buffer voltage falls below a certain point, the behavior of the logic components is 
not predictable and produces uncorrelated pulses for a short period of time. Since 
the whole startup/shutdown procedure has been recorded by an oscilloscope, some 
pulses of Ds3 were not captured since they are very short. This is also the reason 
why the high level of Ds; seems to vary although it should follow the level of Vor. 

Figure 6.52 compares the startup characteristics using the PSCE chip with SECE13 
and SECE23 and the NVC rectifier (see Sect. 6.2.4.3), while driving the shaker at 
PSCE resonance (see Table 6.4). Since a quite large 100 jLF-buffer capacitance is 
used here, there is no overshoot in the buffer voltage as observed in Fig. 6.51a. Dur- 
ing the passive charging, there is no visible difference in the curves for SECE13 and 
SECE23. But as the buffer voltage has approached steady state, SECE23 achieves 


a 0.2 V higher buffer voltage compared to SECE13. The corresponding harvested 


2 2 
power can be calculated as ae = 421 uW for SECEI3 and Us “2 = 462 uW 


for SECE23, so the percentaged improvement is almost 10 % which corresponds to 
the result shown in Fig. 6.38. 

At the beginning of the charging, the buffer capacitor is charged by the bypass, so 
itis obvious that the CMOS diode inside the bypass has a performance which is worse 
than the Schottky diode inside the NVC rectifier. But immediately after the PSCE 
chip has started up, the corresponding buffer voltage curve exceeds the curve of the 
NVC rectifier. From the plot, the PSCE chip with activated SECE23 outperforms 
the NVC rectifier by 237 uW, corresponding to a 105 % improvement. This value is 
much higher than the value given in Sect. 6.4.4.2, for two reaons: First, the chosen 
load resistance of 40 kQ is far away from the optimum resistance (see Sect. 6.2.4.3), 
resulting in non-optimized output power using the NVC rectifier. Second, for the 
analysis shown in this section, the PSCE resonance frequency has been set for both 
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the PSCE circuit and the rectifier. In contrast, for the evaluation shown in Sect. 6.4.4.2, 
the shaker has been tuned to the respective resonance frequency of both circuits. 
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Chapter 7 
Conclusions and Outlook 


A fully autonomous integrated interface circuit for piezoelectric harvesters has been 
presented in this book. By the use of a modified Synchronous Electric Charge Extrac- 
tion (SECE) technique, the conduction losses within the circuit can be significantly 
reduced, increasing the total efficiency. Corresponding Pulsed Synchronous Charge 
Extractor (PSCE) chips demonstrating the superiority compared to standard full- 
wave rectifier circuits have been been fabricated in a 0.35 jum CMOS process. Their 
transistor level design is presented. The chips have been successfully tested under 
laboratory conditions, whereas real piezoelectric generators have been used as a 
power supply. A fully autonomous demonstrator board has been built which allows 
switching between a commonly used full-wave diode rectifier and the PSCE chip, 
whereas a flashing light emitting diode (LED) indicates that the PSCE chip is capa- 
ble of extracting more power compared to the full-wave rectifier. The demonstrator 
board works properly in combination with a custom made piezoelectric generator 
from Christoph Eichhorn which is excited with a 5 m/s? acceleration amplitude at a 
frequency of 114 Hz, resulting in 30 yW output power. 

Extensive measurements have been performed using the PSCE chip in combi- 
nation with a commercial piezoelectric generator from Midé Technologies. All the 
measurements have been driven at an acceleration amplitude of 10 m/s?. By varying 
the vibration frequency, a wide range of output power is achieved. The performance 
in terms of harvested power, efficiency, operation limits and startup behavior has 
been investigated. Detailed breakdown of the different loss mechanisms has been 
performed based on simulations, using a model which corresponds well to the real 
characteristics. The comparison of the measured performance and the simulated 
losses yielded good concordance with the theoretically predicted characteristics. 

The major achievements are already summarized in Sect. 1.4. In short, the PSCE 
chip outperforms standard full-wave rectifiers, offers proper operation over a wide 
range of frequencies, piezoelectric voltages and buffer voltages, consumes low power, 
and achieves high efficiency. Due to the fact that all required signals are generated 
on-chip, and a startup from an uncharged buffer capacitor is possible, the PSCE 
chip enables fully autonomous operation. The improved SECE technique reduces 
the conduction losses, hence smaller inductors can be used while maintaining a 
© Springer Science+Business Media Dordrecht 2015 187 
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high efficiency. This helps reducing the size and thus cost of the whole harvester 
package, and simplifies the integration in real world embedded systems. To conclude, 
the presented PSCE chip overcomes the problems of formerly presented interface 
circuits and is more efficient, contributing to the goal of replacing batteries in small 
embedded sensor systems. 

Several issues in the design and the analysis of the PSCE chip still need further 
research. Some topics will be discussed now. 


e The buffer capacitor voltage level is unregulated and depends on the output per- 
formance of the piezoelectric harvester and the connected load resistance. In order 
to supply electronics, an additional DC/DC converter providing a stable output 
voltage is required. 

e As observed in Sect. 6.2.3, the employed piezoelectric generator exhibits a non- 
linear behavior with regard to the output voltage at frequencies different from 
the resonance frequency. In contrast, the used simulation model is linear. Conse- 
quently, in order to achieve agreement in the piezoelectric output power between 
simulation and measurement, different frequencies have to be set in measurement 
and simulation. Moreover, for the excitation at far off-resonance, the simulated out- 
put power shows deviations from the measured one, as discussed in Sect. 6.4.3.1. 
Hence, a more sophisticated piezoelectric model can be set up, incorporating non- 
linear effects, in order to achieve a better matching between measurement and 
simulation. 

e Theoretical analysis of the switching technique SECEI2 predicts an efficiency 
gain compared to SECEI3 (see Sect. 4.5.2.2), in case the piezoelectric voltage 
is not high enough in order to enable SECE23. However, SECEI2 has not been 
implemented due to the problematic detection of the transition between the first 
and the second phase of the transfer process, as discussed in Sect. 4.6. In future, 
a corresponding detector can be implemented on discrete level in order to verify 
functionality. Afterwards, an integrated implementation can be tackled, including 
an estimation of the corresponding power loss. 

e Currently, the control circuitry is supplied directly by the buffer capacitor which 
has a variable voltage. As discussed in Sect. 6.4.3.2, this results in an increasing 
power consumption for increasing buffer voltage. In addition, a supply independent 
biasing is necessary in order to guarantee a constant switch control timing. If the 
control circuitry would be supplied by a constant voltage, on one hand, the losses 
incorporated in the supply bias independent biasing circuitry and the dependency 
of the power consumption could be avoided. On the other hand, power loss would 
occur in an additional DC/DC converter required to provide a constant supply 
voltage. This matter can be invesigated in the future. 

e The high-voltage power MOSFETs used as switches exhibit a relatively high 
threshold voltage, leading to high conduction losses especially for low buffer 
voltages (see Sect. 6.4.3.2). This issue can be tackled twofold. As a first solution, 
one could think of a DC/DC converter or a charge pump, increasing the gate- 
source voltage of the power MOSFETs in order to reduce the channel resistance. 
As mentioned in the previous point, the additional DC/DC converter losses have to 
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be taken into consideration. Another possibility is to reduce the threshold voltage 
of the power MOSFETS, e.g. by applying the floating gates method [2, 3], which 
could drastically reduce the conduction losses and thus increase efficiency. 
Throughout the book, the vibrations and hence the voltages produced by the piezo- 
electric generator are assumed to be purely sinusoidal. In practical applications 
however, the input vibration source most probably will not be sinusoidal, and will 
contain a spectrum of many different frequencies. It can be worth investigating if 
the SECE technique and the PSCE chip are as performant compared to full-wave 
rectifiers as with sinusoidal voltages. 

As analyzed in [1, 4, 5] and briefly discussed in Sects. 1.4 and 3.4, the SECE 
technique is superior to other power extraction techniques only under certain con- 
ditions, the low coupling factor of piezoelectric generators being one of these. As 
mentioned in Sect. 1.3.2, many different topologies which can be used to harvest 
power from piezoelectric generators have been published up to now. Hence, a 
combination of two or more different techniques in one common interface chip 
can harvest the respective optimum power for the given conditions. 
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Appendix A 
Mathematical Calculations 


A.1 Solution of the Linear Differential Equation Systems 


The initial value problem of the switch configurations 1 and 2 corresponding to 
(4.14), 


o Vp =I 
Ot Cp V V 
= (5) (Os) = zo1 :— te ; (A.1) 
olL Vp L L01 
ðt L 
and (4.17), 
o Vp -— 
Ot Cp V V 
= E (0s) = zoo :— (vr) l (A.2) 
ofr Vp — Vout L LO 
ot L 


can be written as a matrix 


OZ 


(z) (t) = Az(t)+b z(0s) 2zo:— (5 


A.3 
i (A.3) 
with 


being the time vector and 
0Hz —— 
A -( e cr) (A.4) 
I 0 Hz 
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being the state matrix. The b matrices for switch configurations 1 and 2 are 


Zp — (0V5 _,  ._ (Vo 
b= bi = He and Zo = Z0] :— (ure) 


and 


b = b2 -( o ) and Zo = Z2 := p. 


1 
— z Vout Iro2 


The initial value problem (A.3) is a inhomogeneous, linear, ordinary differential 
equation (ODE) with constant coefficients. By means of a coordinate transformation! 


z(t) = x(t) — A^! b, (A.5) 
the initial value problem (A.3) can be written as a homogeneous ODE 
(0;x) (t) = Ax(t) (A.6a) 
with the initial condition 
X(0s) = xo :— zo + A | b. (A.6b) 
The solution of the homogeneous initial value problem (A.6b) is written as 
x(t) = e^'xo, (A.7) 


whereas the matrix exponential function is defined by the power series: 


The solution (A.7) can be inverse transformed through (A.5), leading to the solution 
of the inhomogeneous initial value problem (A.3): 


z(t) = e^! (zo n A~'b) — An". 


Hence, the solutions of the initial value problems (A.1) and (A.2) for switch 
configurations | and 2 are 


' The coordinate origin of x is on the quiescent point of the ODE (A.3), which is indicated by 
! 
(22)(r) = Az(t) +b =0. 
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Zı (t) = g^ (201 + A~'b1) — A^ lbi 


acad l^) (A.8) 
Iroi 


and 


a(t) = e^ (zo j A~'by) = A= bs 
|. At ( VPo2 — UB Vout 
- ( E )« (ox). (A.9) 


whereas the inverse matrix is given by 


—1 


1 1 
es (MAY cas Qr) (s 1) 
E 0 Hz ms 0 Hz —Cp Os 


Rather than using the infinite power series, the matrix exponential function can 
be practically calculated through the matrix diagonalization. The matrix A € R**? 
can be diagonalized if two eigenvalues 4; 2 and two eigenvectors v1, 2 can be found, 
and the matrix exponential function can be written as 


Ait 
e 0 = 
e^! — v-( j PI, i (A.10) 
whereas the columns of the matrix 
V := (vi v2) 


are composed of the eigenvectors of the matrix A. The two eigenvectors vı 2 are 
defined as 


AV1.2 = A4,2V1.2 
& (A — AiaE) v12 =0, (A.11) 


whereas E is the unit matrix in R?*?. If the matrix A is replaced by (A.4), (A.11) 
can be expressed as 


Sjo enc. 
( n Cp )u: = 0. (A.12) 
p —2 


The linear equation system (A.12) has non-trivial solutions v4.» Æ 0 for 
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From that, the two eigenvalues 





Ado = jo with wo:= 


/ LCp 


can be calculated. The corresponding eigenvectors are the non-trivial solutions of 
(A.12). The two rows of this matrix are linearly dependent, because the multiplication 
by Xj /Cp/L of the upper row leads to the lower row. Thus, e.g. the lower row can 
be removed, which results in the following solution set for (A.12): 


1 
v2 =k- with ke R. 
oe an 


Since there are no more demands for the necessary eigenvectors v1 2, k = lis chosen, 
resulting in 


Using these two eigenvectors, the matrices V and V ^! can be expressed as 


V := (vi v2) 


| (ve va) 


and 
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All these preparations lead to the following matrix exponential function: 


Ait 
e 0 zl 


1 ( 1 1 ) (55 0 ) 1 ] e 
5^ : [Cp ; [Cp f —jaot | ` . 
2 Vy EWES KOE e) -ié 
1 ( el@ot e Jeot | 1 j A 
— P 
95 : | Cp jaot ; | CP. .—joot J 
2 —] p eem ] eo J@0 1 =j e 


E z E 

1 eJoot -+ e Jæot J Cs (020 — e Jeot ) 

2 —j E (elo! = e Jof) el?0t NS e Joot 
"EP 

At COS (Wot) J c, sin (Wot) 


" Er sin (cot) COS (Wor) 


Finally, the solutions of the initial value problems for switch configurations 1 and 2 
can be derived from (A.8) and (A.9): 


L å 
Vpo1 cos (oot) — Iroi4/ c, Sin (Wot) 


Z(t) = (A.13) 
Vpo14/ Er sin (Wot) + [91 cos (wot) 
and 
Vout + (Veo — Vout) cos (wot) — Iro, c sin (oot) 
z5(1) = 7 (A.14) 
(Vpo — Vout) y 7 sin (wot) + Izo cos (wot) 
A.2 Flux Property 
It is assumed that the initial value problem 
Ou 
(=) (t) = f(t u(t) (A.15a) 
u(0) — uo (A.15b) 


has the unique solution 
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Fig. A.1 Illustration of the 


Uy .— GQ(At», u 1) 
flux property 






U] := P (At, uo) 


Pu: R > R“ 
tt> y(t). 


The transformation 


p: RUN > RY 
(t, uo) > B(t, Uo) :— Pug (t) 
is called flux of the ordinary differential equation (ODE) (A.15a). If At; € IR, At» € 


R are two arbitrary integration times and u} € R“, u» € IR^ the points as defined 
in Fig. A.1, the following semi-group property holds: 


uz = (An, (An uo)) = (Ati + An, uo) . (A.16) 
Using the specific values 
At := ti —to and Af:=t-t 
and the labels from Fig. A.1, the following relation can be derived: 
m = (Apu) = O(1—11,u1) =O(1 1, (19). (A17 
From (A.17) and (A.16), the desired flux property can be written as 


9 (r =o. i) UTER. (A.18) 
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A.3 Trigonometric Relations 


In the calculations of the switching techniques, the terms sin 
(arccos u) and sin (2 arccos u) often appear. Using complex representation of the 
trigonometric functions, 


l7. | 
sin u = — (cl — ei) 
2j 
arccos u = —jln (u +jv1-— u?) , 


these terms can be rearranged as follows for all u € C: 


- [amen m) u etn) mit n € [1,2] 
J 


sin (n arccos u) — 


1 n 1 
=> (u +jv1-u?) EE | (A.19) 


With n = 1, (A.19) can be written as: 


sin (arccos u) — zl +jv1—u2 — | 
— — u? 


1 MEN UR 
ab^ uds dcus P 
J 


u? + 1 — u? 
Denominator—1 


sin (arccos u) = V1 — u?. (A.20) 


With n — 2, (A.19) can be written as: 


1 1 
sin @areeos W) = y [a - 1 + jan T= ~ 
= - 


1 2u? —1—j2uv1 — u2 
=> rin tae c 
Qu? — 1)" + (j2uvT =u?) 


——M—— 
Denominator=1 


sin (2arccos u) = 2uv 1 — u?. (AAT) 


198 Appendix A: Mathematical Calculations 


A.4 Numerical Calculation 


The exact efficiency can be determined by the energy loss which has to be calculated 
using (4.9) directly from the transient coil current Zz (t) of the respective switching 
technique. The trajectories of the switching techniques SECEI3 and SECE23 can 
be calculated analytically or numerically also for the case that losses in parasitic 
resistances are considered, because their transition condition Vp = OV does not 
depend on the energy loss. For the switching technique SECE12, it is not possible 
to calculate a transition condition in the lossy case, because the transition condition 
(4.41) holds only for the loss-free case (4.39). However, the transition instance t12 
can be calculated numerically as follows. 
The continuity condition 


! 
Fa (Ta, zo) = Fp (Os, Zap) (A.22) 


holds for all switching techniques ab, whereas T, denotes the duration of switching 
configuration a. The solutions F, of the switch configurations n € (1, 2, 3) can be 
calculated numerically or analytically in the lossy case. For the switch configuration 
b, the flux property 

Fp (0s, Zap) = Fp(— Tp, ZE) (A.23) 


holds, whereas Tp indicated the duration of the switch configuration b. Hence, (A.22) 
and (A.23) can be combined: 


Fa (Ta, zo) = Fp (—Tp, zz). 


For the lossy switching configuration ab = 12, this vector equation can be solved 
numerically only for T :— (Ta, Tp) by searching the zero of the function 


Fa(Ta, zo) — Fp (— Tp, ZE) =: Nap (T) 


using Newton’s method. This can be achieved by means of the command fsolve 
which exists in many numerical libraries. The initial values T? := a”, T) are 
chosen to be the known times from Sect. 4.4 for the loss-free case: 


0. 0. 
Ij := tap — to and ly :— te,ab — lab. 


Using (4.45) and (4.46), the zero for the switching technique 12 can be determined 
as follows: 


?) 
I :— y LC p arccos (5) and p ‘= V LC p arccos (: — =) 


Ny (hi, 22) 2= FiGi; 2) — Bo(— 15, Zz). 
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Note 
Newton's method 


pta ATA 


ab (T^). AT" = T" 
( aT ) ): = —Nab(T ) 





requires the derivative (Jacobian matrix) of the function N55 (T), which must be 
either known or approximated numerically by means of finite differences. For the 
given problem, the Jacobian matrix 


0 Nap "ww ONab ONap M 
( aT jar y= (FF mer" 
E oF, T^ 2) T^ 
= (F) az. w. (S azz) 


is known, because the solutions F4 and F, meet their respective differential equations 














OF a oF, 
( 31 ) (t, zo) = fa (F;(t,zo)) | and (=) (t, zo) = fp (Fp (t, zp)). 





So, using the matrix representation of the vector functions 
{,(z) = Agz+b, and f,(z) = Apz+ bp, 


the Jacobian matrix can be written as: 





Nop 
( 2T ) (T") = (Auk. Zo) + Da, ApF p(T, ZE) + bp) 
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